THE PROCEEDINGS OF 
THE PHYSICAL SOCIETY 


Section A 


Vo. 65, Part 10 1 October 1952 No. 394A 


COMPLEX STRUCTURE AND NUCLEAR STRUCTURE 


The following four contributions were read at the Society’s meeting at 
Southampton on 19th December 1951. 


A Theoretical Attempt to Predict the Excited States of Nuclei 
in the Neighbourhood of 7*°Pb 


Pyowio tL. PRYCE 
Clarendon Laboratory, Oxford 


MS. received 21st March1952 ; read before the Society at Southampton on 19th December 1951 


ABSTRACT. **®Pb has a double closed shell of 82 protons and 126 neutrons. The 
classification of the lowest excited states of *°’Pb, °°’ Tl, ®°Pb, ®°°Bi, which differ by one 
nucleon, should thus follow the single-particle model closely. The attempt is made to 
classify the lowest excited states of 2°°Pb, ?°°T1, 2°°Hg, 21°Pb, 7!°Bi and 7!°Po, which differ 
by two particles, in terms of the one-particle model, taking the interaction of the two particles 
as a perturbation. Considerable uncertainties in the order of the energy levels remain, but 
some definite conclusions are possible. ‘The existence of an isomeric state of 7!°Bi (RaE) 
with very high spin (J=8 ?) is predicted. ‘This fits in with known facts about the alpha- 
particle emitting isomeric state. 

- The method can also be applied to 7°°T1, ?°8Bi and the excited states of ?°8Pb which arise 
by lifting one nucleon into the next shell, leaving a hole. It meets with considerable success 
for #°°T]. Discrepancies with known facts in the case of 7°°Pb suggest a more complex 
picture, according to which the first excited state would have J=0 and even parity. Sucha 
state has never been observed, but its existence is not completely excluded by observation. 


§1. INTRODUCTION 

HE shell model proposed by Haxel, Jensen and Suess (1949), and by 

| Mayer (1949, 1950 a) has been so successful in accounting for the shell 
properties and ground state spins of nuclei that it is tempting to subject 

it to further tests. ‘That it throws useful light on the spin of excited states is 
well known from the study of isomeric transitions, but any attempt at a detailed 
prediction of excited states in this domain is in general frustrated, in view of the 
approximate nature of the individual particle description, by the large number 
of levels arising from the various nuclear configurations involved, the uncertainty 
concerning nuclear interactions, and the intricacy of the calculations. ‘This 
complexity is much reduced, at least for the lower lying excited levels, in a nucleus 
which is close to a double closed shell of protons and neutrons. For a nucleus 
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composed of a double closed shell plus one proton, or one neutron, the lowest 


levels are, within the limitations of the model, those of a single particle moving _ 


in the equivalent potential of the core, and constitute no further theoretical 
problem. ‘The same is true, with slight modifications, for nuclei consisting of — 
‘a double closed shell less one proton or neutron. Nuclei differing by two nucleons 


from a closed shell offer a testing ground for more detailed applications of the | 


model, while not yet leading to such a degree of complexity that comparison with 
experiment is impossible. Each extra nucleon may be ascribed a definite energy 
state, and in zero approximation all states arising from this configuration will 
have the same energy. When the interaction between the two nucleons is taken 
into account the degeneracy will be removed, the different resulting energy levels 
being classified by the resultant angular momentum (nuclear spin). 

In order to carry out this programme certain assumptions concerning the 
nuclear forces are necessary. The results will be more sensitive to these if the 
range of the nuclear forces is comparable with the nuclear size (as in the lighter 
nuclei) than if the nucleus is large compared with the range. In the latter case 
it may be expected that the most important parameters are (i) the ‘strength’ 
of the central singlet interaction, and (ii) the ‘strength’ of the central triplet 
interaction, regardless of the exchange nature of the force, while the shape and ~ 
exchange character of the central force, and any tensor or spin—orbit forces, may 
be expected to play only a minor part. The present calculations take into account 
only these two parameters. 

For this reason the paper is confined to nuclei in the neighbourhood of 7°°Pb 
(82 protons and 126 neutrons), as this is the only double closed shell nucleus 
which is (i) sufficiently large in relation to the range of nuclear forces, and 
(ii) in or near the region of nuclear stability. ‘The nuclei of interest are therefore 

(i) 2°°Pb and ?°°Bi, containing one extra nucleon; 

(ii) 74°Pb (RaD), 71°Bi (RaE) and ?!°Po (RaF), with two extra nucleons ; 

(iii) 2°7Pb (AcD) and ?°’Tl (AcC”), deficient by one nucleon; 

(iv) 2°6Pb (RaG), 2°°T'l and 2°%%H¢, deficient by two nucleons; 

(v) ST (‘ThC”), ?°8Bi, with one nucleon and one hole, and the excited states 

or. VeRD itselt. 
Of these, 2°°Pb, ?°’Pb, 2°8Pb and 2°°Bi are stable, ?!°Po emits alpha-particles, 
208Bi probably decays by electron capture, and the remainder are beta-active. 
The nuclear diameter in this region is about 17x 10-" cm (circumference 
5 x 10-" cm), which is to be compared with about 2 x 10° cm for the range of 
the force. 

The biggest stable, or nearly stable, double closed shell before 2°°Pb is 160, 
whose diameter is only 7 x 10- cm, which is not really big enough to justify the 
approximations. ‘I'he energy levels which fill at nucleon numbers 20 and 28 
do not fall into sufficiently separated groups for “Ca, 48Ca or ®®Ni to count as 
double closed shells for our purpose. 


S20) ee P DEAN DE 2 Bt 
Little is known experimentally about the excited levels of these nuclei. The 
spin of the ground state of #°Bi is 9/2 and the magnetic moment indicates that 
the odd (83rd) proton is in an hg), state. The shell model gives hgjp, fy, fs, Pays, 
Piz and iy3/2 as low lying states. Analogy with earlier closed shells would suggest 
the order hyp, frj2, Ps. ++, With is). probably interpolated before psp, or even 
before f,),. The order of the levels in 2°7Pb, a nucleus with one hole in the 
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corresponding neutron shell, discussed in §4, suggests that possibly the order 
is Ngjs, iygja, fyj9.... In lighter nuclei where this neutron shell is starting to fill, 
however, the ground state of nuclei with 83 neutrons is f,j._ There is indirect 
evidence from the following work, and from the levels of #44Bi (AcC), which 
should be similar, that these levels are close together. Information on the excited 
states could in principle be obtained from the inelastic scattering of neutrons 
in bismuth, from the reaction ?°°Pb(dn)?°Bi, and from gamma-radiation following 
electron capture in ?°Po. The energy needed to dissociate 2°Bi into 2°°Pb 
and a proton can be calculated from a closed cycle (Kinsey et al. 1951). It is 
3-79 Mev. 

The ground state spin of 2°°Pb is not known. Information about the excited 
states is available from the reaction ?°°Pb(dp)?°°Pb (Harvey 1951a). ‘There are 
proton groups corresponding to excited states at 0-75 Mev, 1:56Mev, 2-03 Mev, 
and 2:-54Mmev. ‘The dissociation energy into 2°°Pb and a neutron is 3-87 Mev. 
The resolution of this experiment is only 0-5 Mev, and some of these levels may 
be unresolved multiplets. According to the shell model the lowest group of 
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Fig. 1. Assumed energy levels of ?°°Pb and 2°Bi. 


states 18 iy1J, Zoie, rie, As/2, Asie, Sia aNd jysj- ‘The spin-orbit coupling, which 
is presumably large for the high angular momenta, may bring jj5). well down 
towards the bottom of the group, and may push ij,). above go. ‘The ft-value 
of the beta decay of ?°°Pb is compatible with either go). or 1,;/, for the ground 
state, but not with j,;).. The low intensity of the proton group corresponding 
to the first excited state at 0-75 Mev (14.Mev deuterons, 15 Mev protons) might 
be taken as an indication of high angular momentum transfer, and suggests that 
this is the jj5j. level. 

These considerations suggest that the ground state is go (or iyy)/9), the 
0-75 Mev state is js, the 1-56 Mev state is i,4). (or gg.) and the 2-03 Mev state is 
dsj. OF g;jo, but the evidence is very inconclusive. It is summarized in tentative 
form in fig. 1. Further information could be obtained from the study of the 
so far undetected gamma-radiation following the decay of ?°'T'l. 


§ 3. 2°Bi (RaE) 

We shall treat ?!°Bi as a system consisting of a neutron and a proton moving 
in the equivalent field of the closed shell. Concerning this field we shall assume 
that it is spherically symmetrical and gives rise to the one-particle states discussed 
above. This means, of course, that it contains strong spin—orbit interaction to 
give the necessary doublet splittings. We shall treat the perturbation arising 
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from the proton—neutron interaction as small compared with the spin-orbit 
splitting (in other words, we adopt jj-coupling). This is probably a tolerabel 
approximation, but a quantitative estimate indicates that it may break down when 
the proton and neutron wave functions overlap strongly. 

We shall further assume that the perturbation in the energy is proportional 
to the probability p (per unit volume) that the neutron and proton are in each 
other’s immediate neighbourhood. If we were dealing with spinless particles 
whose wave function is ‘V’(r,, r,) this would be 


pe PTE e rs. on ee ae eee (3.1) 


The effect of spin is considered below. 

The justification for this assumption is not very firmly based. If the 
proton-neutron potential V(|r,—r,|) were weak, we could, ignoring spin as 
above, and now also any exchange character of V, write the first-order 
perturbation as 


Ex =f (Pee, PV Gr yer, ee ee eee (3.2) 


and for sufficiently short range of force we could replace V(r)Ad%(r), giving Ap. 
Actually the potential is so strong that (3.2) is a poor approximation. We shall ~ 
nevertheless retain our assumption that the perturbation energy is Ap, for lack 
of a more precise method of calculation, but the ‘strength’ will not be assumed 
to have so direct a connection with the interaction potential as above. 

It is clear that if V contains a coordinate exchange operator (as in the 
Heisenberg or Majorana potential) this will make no difference to (3.2) in the 
limit of zero range. 

If now we take the spin into account, we can split the joint wave function into 
a singlet part (S=0) and a triplet part (S=1). We now have to introduce two 
‘strengths’ A, and A, for the singlet and triplet interactions. If we denote by 
p, and p, the probability of the neutron and proton being in each other’s 
immediate vicinity in the singlet and triplet states respectively, then we can write 
the mutual energy as 

Ep Peet hie. yee 2 ee (3.3) 


To evaluate p, and p, we proceed in two steps. ‘The unperturbed wave function 
is in terms of jj-coupling, and this has first to be translated into LS-coupling. 
Then for each resultant L which occurs, the appropriate p must be calculated, 
and these summed over the relevant values. Cross-terms involving different L’s 
do not occur, since for the central forces we are assuming Z and S are constants 
of the motion. 

The neutron will be assumed to be in a state of angular momentum j, 
compounded of an orbital momentum / and spin $ (7 =/+ 4); and similarly for 
the proton we assume j’ and /’. ‘The resultant angular momentum of the system 
is I. We denote the state of the system by |J; j, j’);; to emphasize its jj-coupling 
nature. It can be expressed as a superposition of states of definite total orbital 
momentum JL and definite total spin S which we denote by |J;L, S)z, to 
emphasize their L,S-coupling character. Since S only takes on the values 0 and 1, 
the sum consists of only four terms: 


ALS 9,9 Dig = VoL, 0) gt B31, 1)z5 
+ 4/y|E; 141, 1)754+ VOU; 7-1, 1)z5. die 6 ste . (3.4) 
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The calculation of the probability amplitudes ./« etc. is fairly straightforward. 
For the case j=/+ 4, j’=l'+4, their squares, « etc. which are what we shail 
need, are 


_ (tlt l'+2)1441-D ] 
7 2(21+ 1)(2/' + 1) : 
(J—V)PT+ 14+ +2)14+ 41-2 
2(22+ 1)(21’ + 1)1(1 +1) : 
— 7+ 1tl-lyi+14l-Dl4+-Dls+l4+1-D 
a 2(22+ 1)(27' + 1)(2+1)(22 + 1) ‘ 
ae (f+ l- EYL + DT ++ + T+ bl 42) 
2(21-+ 1)(20’ + 1)1(27 + 1) 
The expression for the other combinations of j=/+4, j’=l' +4 can be obtained 
from (3.5) by the simple pansion mato’ of pe hane lL by —(/+1) if j=l-4 
and /’ by —(/'+1) ifj’=l’—-}. 

The wave function of a single neutron can be written as the product of a 
radial function %,(7,) and the appropriate linear combination of products of 
spherical harmonics and spin functions, and similarly for the proton. Each 
of the terms in (3.4) will therefore be represented by a wave function with the 
radial dependence ¥%,(7,)/,(7)), which is multiplied by the appropriate angular 
and spin wave function. ‘The wave function of specified LZ (and m, = M), ignoring 
the spin part, is of the form 

Pal'ap(%>) Chat m YO, dn) Y,"' (9, dy). 
Its contribution to p is 

JF |S CHV", $)¥v"'(8,b) [2 sin 8 dO def Whe2(r)hy2(0) dre ee. (3.6) 
Let us call the angular integral in (3.6) 47P, and the radial integral Q/47. 
O is the same for all the terms arising from the same configuration, depending 
only on the radial wave functions, while P depends only on L,/ andl’. ‘The 
factor 47 has been included to normalize P(L;/,/') in such a way that it is unity 
if the angular correlation of the two nucleons is purely random. P is therefore 


a measure of their correlation in solid angle. 
It is given by (cf. Carlson and Rushbrooke 1950, eqn. 17) 


(21+ 1)(2U' +1)(2g—2L)! (2g — 21) ! (2g — 21’) !(¢!)2 
2g+I)iig—L)l(g-)Me—L)! 


B= 


P(L;1,1) = 


if 2g=L+/+l' is even and 
Men iat Veet). ot 2 tere Vin of ih Aiedie os (3.7) 
if L+1+I' is odd. 
The interaction energy is therefore given by 


E, =%P(;1,1')A,O + [BP(U31,U)+yP1+1;h,1')+8PU-1;1, 0) AQ, 


in which the quantities P, «, 8, y and 6 are directly calculable. For convenience 
we write 
=A,Q, e,=A,Q, 
and ee eels) een) de® Pull Pele A iro a iubtyr- esta (3.9) 
a,=BP(I;1U)+yP04+1;h1)+8PT7-13;1,0), 
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so that the expression for E,,, becomes 


Exp = Gear Gey 0 EES (3.10) 


The quantities «, and «, have the dimensions of energy. They are determined 


by the degree of overlap of the radial wave functions of the neutron and proton, 
and by the strength of the proton—neutron potential. ‘The numerical quantities 
a, and a, depend, for a given configuration (i.e. j, j’, J, l’ given), only on the 
resultant angular momentum J. 

At this stage it is instructive to make a rough estimate of the order of 
magnitude of the <’s. If we consider those single-particle states in which the 


J-value appears for the first time in the enumeration of the states in order of — 


energy, the radial wave functions have no nodes, and their maxima are fairly 
sharp and overlap considerably (e.g. the configuration (j15), hgj)). Such 


configurations will have the largest value of Q. Assuming that the nuclear — 
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Fig. 2. Dependence of the interaction energy Enp in *!°Bi on nuclear spin J for the nine configura- 
tions considered. Energy scale is arbitrary. (3a¢+2as drawn to same scale for each con- 
figuration, corresponding to the assumption et/«s=1°5). 


radius R is about 8-4 x 10-%cm, and that the effective ‘spread’ S of the radial 
wave function is about 2 x 10- cm, one estimates 
(47Q)-? ~s(R— 4s)? ~110 x 10-9 cm?. 

Further, if we suppose that V is adequately represented by a square well 
35 Mev deep and of range 2x 10cm, and for A we take 47fr?V dr, then 
A/4r=--93 x 10-*°Mevcm®. So «e=AQ=—0-85 Mev. We shall take this figure 
as a guide to the order of magnitude only. For configurations in which there is 
little radial overlap it may be much smaller. 

‘Table 1 shows the numerical parameters a, and a, for the nine configurations 
formed from jp, jisj2. lisja for the neutron, and hgj, ij3/2, fj. for the proton. 
in order to show the relative behaviour inside each configuration in graphical 
form, the definite value 1-5 has been chosen for the ratio ¢,/e,=A,/A,. This 
value corresponds roughly to the relative strengths of Yukawa potentials of 
equal range which will account for the low-energy proton-neutron and 
proton—proton scattering data. Having fixed this ratio, E,, now depends only 
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d scale factor. Figure 2 shows E,, 


nfigurations, the energy scale being arbitrary. 
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The energy of the complete system is the sum of the energy of the closed shell, 


the energy of the neutron, the energy of the proton, and the interaction energy : 
E=Eybhgt lgthe: o-) 1 werner (3.11) 


Owing to the uncertainty discussed in §2 concerning E, and E, only limited — 


conclusions can be drawn, but certain points emerge clearly. 

It will be seen from fig. 2 that in each configuration (i) the lowest energy is 
attained for the lowest possible value of J, (ii) even and odd values of J behave 
differently, (iii) proceeding from the lowest /-value in steps of 2, the energy 
increases monotonically, (iv) descending from the highest J-value in steps of 2, 


the energy first increases to a maximum and (with one exception) decreases again. — 


From this it can be inferred that the ground state of RaE corresponds to the 
lowest /-value of some configuration. 

Another conclusion of great interest is that there must be an isomeric state 
of very high spin (at least 8). This follows from the fact that in each configuration 
the state of maximum J is lower in energy than all states whose J differs by less 
than 7 units, and it is a fair inference that the maximum J-state which has the 
lowest energy among the various configurations is below all other states which 
do not differ in J by less than 6 or 7 units. Such a state will obviously have an 
astronomical lifetime against gamma-ray emission, and against beta-emission to 
any plausible energetically allowed state of 74°Po. Its most probable mode of 
decay will be by alpha-particle emission to ?°°Tl. An alpha-emitting isomer of 
RaE is in fact known (Neumann, Howland and Perlman 1950), whose energy 
above the ground state is estimated to be 0-18+0-ll1Mev. Of the nine 
configurations which have been studied all have a maximum J/-value of at least J =8. 
Consideration of the single-particle energies (£,,, E,) and the overlap of the radial 
wave functions make the following combinations possible competitors for the 
isomeric State: (Boj0, hojo)r-9, (Bola tusie)r—11» (Borer frje)r—a» (Jrsiar Moje) r—r2» 
(Jas/2> l1s/2)7—14- 

A certain amount is known experimentally concerning the energy levels of 
RaE. It is formed in the natural radioactive sequence by the beta-decay of RaD, 
leading to an excited state. The transition energy is very low, probably about 
18 kev (Bannerman and Curran 1952). No measurable decay to the ground state 
takes place. The lifetime of RaD (taken with the assumption J/=0, even, for 
RaD) is compatible with J=0 or 1 for the resultant state of RaE. There are at 
least six soft gamma-rays following the beta-decay, of which the most energetic 
has an energy of 46-7kev (Tsien 1946, Feather 1949). This indicates at least 
three intermediate states. In addition there is the isomeric state already 
mentioned. Also the work of Harvey (1951a) with the reaction ?°Bi(dp)?!°Bi 
indicates levels at 1-65, 2:2 and 2:°3mev. The present analysis would suggest 
that these ‘levels’ are the unresolved configurations (or rather the densest part) 
formed from the proton in its ground (hg) state and the neutron in the states 
corresponding to the 1-56, 2-03 and 2-54 Mev levels in 2°Pb. (The absence of a 
group corresponding to the 0-75 level can be ascribed to a combination of the 
high spin and the large splitting of the configuration). The work of Kinsey, 
Bartholomew and Walker (1951) on the neutron-capture gamma-rays of bismuth 
also throws light on the question, and is discussed more fully below. 

The spin and parity of the ground state of RaE have been deduced from the 
shape of the beta-spectrum (RaE-> RaF) to be J=2 (even). This assignment, 
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though plausible, is subject to some doubt as it rests on assumptions concerning 
the form of the fundamental beta-decay interaction. If it is correct, the ground 


State 18 (So/2, ligi2)7—2- 
The state in which RaE is formed in the beta-decay of RaD is one of the 


following : (8g)2, Moiz)r—o» (Jusias tisi2)r—15 (Boles frlz)r—19 (ines Dore)r—1» (inie» tusie)r—1- 

The assignment of the intermediate states is narrowed down by this analysis, 
but much ambiguity remains. 

In the configurations we have studied there are several states with the same 
parity and same value of J. If any of these fall close together in energy on a 
first-order calculation, then second-order effects may be important (configura- 
tional interaction), with the result that the levels are pushed apart and their 
character mixed. This further increases the ambiguity of the interpretation. 

The gamma-ray spectrum resulting from the capture of thermal neutrons 
by 7°°Bi (leading to RaE) consists of a broad line of energy 4:17 Mev (maximum) 
(Kinsey, Bartholomew and Walker 1951). A comparison of this energy with the 
binding energy of a neutron in RaE obtained from a closed cycle shows a 
discrepancy of 0-37 Mev (Kinsey, Bartholomew and Walker 1951), the binding 
energy deduced being 4:54mev. A comparison of this value with the value 
3-87 Mev in *°°Pb gives some indication of the magnitude of the proton—neutron 
interaction (4:54—3-87=0-67 Mev, not all of which is necessarily due to E,, 
since the field of the core is altered slightly). On the basis of the present 
calculations one would not expect the capture radiation to lead to the ground 
state, but rather to a group of excited states, for the spin of the capture state 
(Syj2, Hygjz) is either 4 or 5 (odd), so that gamma-emission is most likely to states of 
spin J =2 (odd), 3, 4, 5, 6, (7). The observed width of the line (0-07 Mev more 
than the instrumental width) and the energy discrepancy of 0-37 Mev can be 
fitted into the scheme which we have obtained above, as may be seen from fig. 3, 
which gives a highly speculative picture of the levels. ‘The subsequent 
de-excitation of these levels is presumably by a cascade of soft gamma-rays (which 
would not have been detected with a pair spectrometer), and may lead either to 
the ground state or to the isomeric state (in conformity with its known mode of 
production). 

To show that the various data can be reconciled with the ideas outlined above, 
fig. 3 shows the level scheme which results from a definite choice of numerical 
values for the parameters. For convenience levels with different /’s have been 
separated. ‘These are far from being uniquely determined by the data, and the 
figure should not be taken as more than extremely speculative. In particular, 
no significance should be given to the precise order of the levels. It will be seen 
that the known facts fit satisfactorily into the scheme. In the 7°Bi(dp)?!°Bi 
reaction the groups of levels reached are those marked a and d ; in the 
capture radiation of bismuth it will be mainly a selection of a levels in the spin 
range J=3-6. ‘The parameters have been chosen so as to bring five levels into 
the lowest 50 kev in a way which might be fitted to the known gamma-rays. * 
This coincidence of so many levels in a narrow energy range appears to be the 
result of chance, and this is a rather unsatisfactory feature of the proposed scheme. 
It will be noticed that three of these five levels have spin 1, the other two having - 


' * The level scheme proposed by Bannerman and Curran involves six levels. More experimental 
work is still needed to decide the level scheme, but if a sixth level is required, (i,;/2, hg/.)7—4 might 
supply it. 


782 IM dio Mery ce 


I=0 and 2. The existence of possible isomeric levels with an energy around 
0:2 Mev, however, follows naturally from the parameters chosen to give the 
correct discrepancy in the capture radiation energy, and speaks in favour of the 
scheme. The values of «, chosen lie in the range 0-1—0-5 Mev, in satisfactory — 
agreement with the calculated order of magnitude. . 
The situation is further complicated by the possibility of low levels arising 
by raising a nucleon from the core and pairing it off with one of the odd nucleons 
in a state of high /, as discussed at the end of §7. : 
Our result, that in any configuration the state of minimum spin has the lowest 
energy, contradicts a rule enunciated by Nordheim (1950). This states that in ~ 
odd—odd nuclei the lowest level is that of minimum or maximum J according to 
which makes the spins of the neutron and proton most nearly parallel. 
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Fig. 3. Very tentative diagram of energy levels of 7°Bi arising from the nine configurations 
considered. The figure is to be interpreted only as suggesting in a general way the predicted 
complexity of the level scheme. Details are not meant to be correct. Key: A=(gojp, i132); 
B=(jisja5 hgj2), C= (jrs/2, £72), D=(in,2, 113/2)>. a=(go2, hg)2), b=(g9)2, fr )2), c=(ji5)2, ly3/2)> 

=(iy1)2, hyo), e=(iyy/2, f7/2). Capital letters indicate even states, lower case odd states. 


Empirically it is fairly reliable for light and medium nuclei, though it does not 
hold without exception. Our formulae for the interaction energy show that for 
reasonably low j-values the state of maximum spin is indeed the lowest level in 
the cases where Nordheim’s rule demands it, provided a suitably high value 
of e,/e, is taken, but that for high j-values the value of ¢,/e, required to make the 
rule hold is impossibly high. It is quite possible that the effect of the range of 
the force, which is more important in lighter nuclei, favours the high spins. 
In any case, in view of the known exceptions, the present results should be taken 
rather as elucidating its range of validity than as a contradiction to it.* 


* Actually, even for lighter nuclei the present results, although not directly applicable to nuclei 
far from closed ‘shells, seem to be in better agreement with experiment than Nordheim’s. In 
‘Al, for example, the present calculations applied to (d5/)? give [=0 lowest, in agreement_with the 
superallowed character of the beta-decay, as contrasted to Nordheim’s J=5. In ®8Co and ®Co 
(f7/2, Ps/2), the present calculations indicate that I=2, J=5 are very close (crossing over for 


«t/€s=1-:59), in agreement with known isomeric states of probable spins 2 and 5, whose order. is 
opposite in the two isotopes. tas 
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§ 4. Pb (RaD) AND Po (RaF) 


These nuclei will be regarded as two neutrons, or two protons, moving in the 
field of the core. The identity of the two particles introduces a small modification 
to the calculations of §3. The wave functions must be made antisymmetrical, 
and this means that the two nucleons can only be close together in the singlet 
state, so that only the singlet interaction is effective. The mutual energy is 


2uP(I;1, l')e, mabe nara} 


(except when j =j’, /=/', when the factor 2 is absent). A general feature of (4.1) 
is that the interaction energy is zero for alternate values of J. 

If we consider RaD, and use the assignment of neutron levels discussed in § 2, 
it follows that the ground state probably arises from (ggj.)?;_9. ‘The remaining 
states, J=2, 4, 6, 8, of this configuration (odd values are forbidden by the Pauli 
principle) are higher and follow in that order. The gap between 0 and 2 is 
considerably greater than between the remainder (the numerical values suggested 
by the fitting of the RaE levels would indicate that J=0 lies about 0-5 Mev 
below J=2). There will be other levels corresponding to (29), jy5)2.) of which 
those with J =4, 6, 8, 10, 12 will have an energy about 0-75 Mev above the ‘head’ 
of (g9/2)", while those with 7=3, 5, 7, 9, 11 will be lowered by neutron—neutron 
interaction (I=3 lowest). Similarly there will be levels arising from (j,5)2)?, 
whose head will be about 1-5 Mev above the head of (gg).)?.__ Owing to the high 
j-value and the concentrated nature of the radial wave function, (j,5/2)?;-9 may 
be strongly depressed and may possibly even compete with (g9j.)?;-9 for the 
ground state. 

RaD is formed in a highly excited state in the beta decay of RaC’(?!°T'l). 
Gamma-rays of energy 3-1, 4-0 and 4-9 mev are known. If we apply the present 
ideas to RaC’, regarding it as a proton hole (in an sj). state) and three neutrons, 
and assume the three neutrons in their ground state to act as one neutron in a 
Zoo State, then we conclude that the ground state is [=5 (even). Since the 
beta-decay appears to be allowed, RaD must be formed ina state of spin at least 4. 
The gamma-rays therefore presumably do not go direct to the ground state. 
In the alpha-decay of RaC’, RaD is formed in the ground state, and since its 
formation from RaC” is a rare mode of decay, nothing seems to be known 
experimentally concerning its lower excited levels. It is therefore not profitable 
at present to work out anything in great detail. 

Similar results apply to ?!°Po, with two protons. Since, however, there is 
some reason to believe that the hgj., i,3/. and f,). states are closer together than the 
corresponding neutron states, there will be more overlap of the configurations. 
In particular, it is not clear which will be the lower of the two states (hgjp)?;~o, 
(i}3/2)27-9, and they may in fact be quite close together (though configurational 
interaction will tend to push them apart). The claim has been made by 
Zavelsky, Umarov and Matushevsky (1949) that the beta decay of Rak is 
complex, leading most of the time to an excited state of 7!°Po 85 kev above ground. 
Although their interpretation of their experiments does not seem very convincing, 
such an excited state could be fitted into the theoretical scheme. It would not 
emit gamma-rays. Otherwise the only known gamma-ray emitted by ?!Po is 
that following electron capture in 7!°At, of energy about 1 Mey. 
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$5. ’Pb anv 2°?T1 (AcC”) 


These nuclei can be treated as holes in the closed shell. The energy required 
to remove a particle, leaving a hole in a definite state, will roughly speaking be 
the one-particle energy of that state, plus the ‘pairing’ energy, which in the 
short-range approximation is the same as the interaction energy (j+4)e, of two 
similar particles in the same orbital with resultant J=0 (Mayer 1950b). The 
low-lying energy levels of these nuclei will therefore follow approximately inversely 
the order of the one-particle states. The available low levels for ?°’Pb_on the 
shell model are p4j., Psy, f52, fr/2, Hog and iygja, and for *°’T'l syjo, dgjo, dsj, Brie 
and hij. For 2°’Pb there is enough experimental material to fix all these levels 
except hg. The level diagram is shown in fig. 4 (the reasons for this assignment 
of levels are given in the Appendix). The hgj level is probably somewhere about 
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Energies are in key. 


2-3mev. ‘The levels shown above 2 Mev (which are formed by electron capture 
in °°7Bi) are more complex in structure, involving excitation of other nucleons 
in the closed shell. 

*°7'T'l is known to have an excited state 0-35 Mev above ground (long-range 
alphas of AcC), but experimental evidence beyond that is meagre. The ground 
state and first excited state are presumably s,). and dg). respectively. 


§ 6. AD9 2) os 206'T"] AND AN 8 is3 


‘The methods already developed for RaD, RaE and RaF apply to these nuclei. 
‘The fact that we are dealing with two holes instead of two particles makes only 
minor differences. ‘The energy of a single hole must be taken to be minus the 
energy of the particle in the corresponding state. ‘The interaction energy, on 
the other hand, has the same sign for two holes as for two particles, and is given 
by the same formula. It is also readily verified that, more generally, in 
jj-coupling, the interaction of a hole with a particle (as in 2°°T] or 2°8Bi, or the 
excited states of ?°8Pb) has the same sign, so that no distinction need be made. 
between holes and particles in discussing the effects of mutual interaction. 
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We shall discuss only ?°*Pb in detail. Its levels are based on those of 2°7Pb 
(one hole), the lowest five of which are well established. Fifteen gamma-rays of 
*°6Pb are known (Alburger and Friedlander 1951). There is therefore some 
hope of correlating theory and experiment. 

Let us first consider the energies of the heads of the several configurations. 
Were it not for the pairing energy, these would in zero approximation (neglecting 
configurational interaction and the effects of the difference in the effective field of 
the core) be given by adding the two corresponding single-hole energies. This is 
shown in the second column of table 2. The third column shows the coefficients 
a, of «,, in the expression (4.1) which gives the effect of mutual interaction. 


Table 2. Levels of 2°*Pb 
Energy of head 


Configuration (uncorrected) <_______-___—- a, > 
(Mev) a0) 1 2 3 4 5 6 a 
{P2)* 0:0 1 
- 6 
{Pres f5/2) 0-565 Z 0 
a37 
(Pij2, Psj2) 0-87 0 : 
24 2 
- _— _ 
(£52) ids 3 = = 
8 
(P12) £72) lctal 0 ; 
NP 8 
(f5/2, Ps/2) 1-43 0 ae 0 - 
i 14 
(Pi2; 13/2) 1-62 0) ae 
8 40 200 
(£5/2, £72) LO:7, 0 = 0 a 0 ne 
2 
(P3y2)” 1-74 2) : 


The correction for the pairing effect in the levels of 2°’Pb is the sum of the 
pairing energies of the two states in the configuration (a positive correction), 
together with a negative correction to allow for the change of zero (conventionally 
the ground state of *°*Pb). Theoretically the pairing energy in an orbit 7 is 
(j+4)e,, being evaluated for the configuration j*, so that in principle the pairing 
correction and the configurational splitting are not independent. 

The following results emerge. ‘The ground state is (pjj.)?;-9. The first 
excited state is either (pyjo, f5j2)7-2 OF (Pil) Psi2)r—2 OF possibly (f;/2)?;-o. All the 
states up to at least 1-5 Mev are even states (the first odd states are (pj), ij3/2) with 
I=6 and 7). 

There is experimental evidence that the first excited state is at 803 key. 
It results from the emission of short-range alpha-particles from ?!°Po, which 
is followed by a single gamma-ray (Alburger and Friedlander 1951, Grace et al. 
1951). The spin is inferred to be J=2 (even) from the alpha~gamma angular 
correlation (de Benedetti, private communication to H. Halban). 'This:evidence 
does not exclude states of parity opposite to that of J being lower, since these 
could not be reached in the alpha-decay of Po; nor does it exclude J=0 being 
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lower, as such a state would not emit gamma-rays (though it is unlikely that the 
resultant orbital electron emission would have been missed). Nevertheless it 
seems reasonable to assume that the first excited state is indeed at 803 kev with 
I=2 (even). Although the gamma-rays have been resolved to 1:2%, and their 
energy measured to $%, this accuracy is not sufficient to fit them unambiguously 
into a term diagram. The present considerations help to some extent, but very 
considerable ambiguity remains. An attempt at a level scheme bearing the 
various factors in mind is shown in fig. 5. No part of this is to be regarded as 
more than very tentative. This is particularly true of the assignment of the first 
excited state of (py, P3j2), which makes some of the selection rules easier than 
(P12. £52) but the indications are very slight. On the other hand, the picture that 
this gives of many levels within the first 2 Mev is probably reliable. 


87. 28 2B VAND LT AE EXCL EDs WATE Ss Ome 2 Pb 


208'T'] (ThC”) can be regarded as made up of a proton”hole and a neutron. 
The relevant one-particle levels are therefore those of *°”T'l for the hole, and 
of 2°°Pb for the neutron. We therefore expect the lowest levels to arise from 
(S1j25 Zg2) and (dsj, 9/2)—which should form a doublet and a quartet respectively. 
The levels of ThC” have been exhaustively studied (Surugue 1937) and the 
experimental evidence is in good agreement with this assignment (though other 
assignments of the neutron state would fit just as well). ‘Table 3 lists the 
experimental levels (line 1), the theoretical spin (line 2) and calculated values of 
a, and a, (lines 3 and 4). All the states have even parity. ‘The intensities of the 
long-range alpha groups of ThC fit in well with the predicted spins. It will be 
noticed that the quartet is reasonably fitted by the values e,= —225kev, 


Table 3. Levels of ?°8TI 


Experimental energy (kev) 0 40 328 472 492 616 
Predicted spin 5 4 3 6 4 5 

a, 1-000 0-444 1:667 0:734 0:495 0-394 
ds 0:000 0-556 0:000 0-629 0-202 # 0-000 


e,~—130kev. ‘The doublet splitting indicates «,—«,= —72kev (the «’s are of 
course different in the two cases). ‘The spacing between the two groups is about 
150 kev more than would be calculated using the 350 kev spacing in 2°7T1. 

208Bi is similarly made up of a neutron hole (levels pj), f5/2, pj etc., as in 
207Pb) and a proton (levels hgj, 143/2, fy. etc., as in 2°°Bi). Since the first three 
proton levels are probably fairly close together, and the first two levels of 2°7Pb 
are only 0-54 Mev apart, there are probably at least 24 levels of 2°°Bi below 1 mev. 
The ground state is probably /=5 (even), 

Only negative evidence exists concerning ?°SBi. It is supposed to be unstable 
against electron capture to °°*Pb, with a very long life. The high predicted spin 
value supports this. 

One would imagine that the lowest excited levels of the closed shell nucleus 
28Pb itself would arise from the excitation of either a proton or a neutron, and 
that they could be treated in terms of a particle and a hole. The high energy 
(2:62mev) of the first observed level would at first sight appear to be 
understandable on this view, but a more detailed comparison with known 
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facts reveals difficulties. For one thing, the difference in the binding energies 
of a neutron in ?°Pb (7-38 Mev) and in 2°Pb (3-87 Mev), and of a proton in 
208Pb (8-05 Mev) and in ?°Bi (3-79 Mev), suggests that about 3:5 Mev and 4-2 Mev 
is required to lift a neutron and a proton respectively to a level in the next shell. 
It is therefore difficult to understand in this way a level as low as 2-62 Mev. 
In addition, the lowest levels would be expected all to have either rather high 
spins, or odd parity, or both. If a neutron is excited, for example, one would 
expect the levels to arise from (py), Z9/2) : 1 =4, 5 (odd) or (f5/0, Zoi2): 1=2—7 (odd). 
ieaeproton, then’ (s),,/ hye): 2=4; 5° (odd), <(83, “ijg) 7 2=6, 7? (even), 
(Syjo, fy): 1=3, 4 (odd), or (dsj, hg): =3—7 (odd), etc. There is evidence, 
however, that the 2-62 Mev level has J=1 or 2 (even). It is clear that none of the 
above combinations can explain this. 

A possible way out of the difficulty is to suppose that the lowest excited levels 
correspond to the excitation of two nucleons, which both go into a state of high 
angular momentum and pair off with zero resultant. The pairing energy gained 
in this way can be very high (e.g. for the j,;,. neutron state it might well amount 
to 5 Mev), so that it may be energetically more economical to excite two nucleons 
than one. The different excited levels would then correspond to the configuration 
of the holes out of which the two nucleons came. The resulting level diagram 
would bear some resemblance to the level schemes of 2°6Pb, 7°°Tl and ?°*Hg, all 
raised up in energy and superposed. On any such view, the theoretically expected 
level diagram is extremely complex. 

This explanation is not without its difficulties, however, for one inevitable 
consequence is that the lowest excited state should have J=0 (even). No 
experimental evidence exists for such a state, which would presumably emit 
electron—positron pairs and orbital electrons. Its energy would presumably be 
about 2 Mev and, if it were excited, its characteristic radiations would be easily 
detected. If it exists, it may be that failure to detect it arises from its not being 
formed to any appreciable extent following the beta decay of ‘ThC’—which 
probably leads in the first instance to a state of fairly high spin, from which it 
would not be fed—and that when the 2:62Mev is reached the overwhelming 
majority of gamma-rays lead to the ground state. 

Although this conjecture may appear fanciful, one must remember that in 
another double closed shell nucleus, 1°O, the first excited state has J=0 (even) 
and emits pairs—presumably for very similar reasons to those outlined above. 

If this explanation is correct, one might expect a similar phenomenon in RaE. 
Raising one neutron from the core to pair off the odd neutron in a J,5/. orbit, for 
instance, may give rise to a low energy state. The RaE problem must therefore 
be considered as still open. Similar levels can also occur in 7°°Bi and 7°°Pb. 


§8. DISCUSSION 

In this paper an elaborate deductive structure has been erected on extremely 
flimsy foundations. Not only has it been assumed that it is safe to calculate 
nuclear energies as though one were dealing with an essentially independent- 
particle model—and the recent successes of the strong spin-orbit model indicate 
that this approach is valid beyond what appears to be superficially reasonable— 
but use has been made of definite assignments of orbital states and energies of the 
individual particles, often based on very little direct evidence. ‘The intention of 
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this work is frankly exploratory, and for that reason the methods of calculation — 
have been as simple as possible. This has been achieved by assuming that the 


forces are of very short range, ignoring non-central forces and assuming rigorous 


jj-coupling. ‘The apparent success of much of the resulting theory in reconciling _ 


available experimental data must not be taken too seriously as support for the 


hypotheses. Other assignments might do just as well or better, and fatal | 


discrepancies may yet be found. Further analysis, both theoretical and 
experimental, is now required. 
On the theoretical side, much can still be learnt from the apes study of 


the available experimental material. The present work was started two years — 


ago with the view of understanding the energy levels of RaE, as a modest beginning 


to a study of nuclear energies generally, but no real progress was made until the — 


interpretation of the data on 7°7Pb was found. This gave sufficient confidence 
in the general predictions of the shell model, and definite enough results, for an 
extrapolation to the next shell to be attempted. Not until further qualitative 
work of this kind has been carried out will it be profitable to refine the methods 
of calculation. 

On the experimental side, more work to test the predictions of the theory 
would amost certainly yield results. In particular, more work on the gamma-ray 
spectrum of ?°*Pb, either with greater precision or using coincidence techniques, 
is to be welcomed. Further study of the gamma-rays of RaE (i.e. following the 
decay of RaD), with the view of establishing the spins and parities of the levels, 
would be useful. A direct determination of the spin of the ground state of Rak 
would be an invaluable check on this work and on beta-decay theory. Much 
higher resolution in (dp) experiments would be capable of giving extremely 
valuable information. ‘The levels revealed with the present resolution are in 
most cases whole groups of levels. The development of a similarly fine technique 
with (dn) reactions would also contribute to the understanding of nuclear 
structure. An attempt to find a possible pair-emitting state in ?°°Pb (and possibly 
in 2°6Pb and #1°Po), perhaps by direct excitation with fast neutrons or gamma-rays, 
would throw light on the difficult problem of the excited states of a closed shell. 
An analysis of ee gamma-rays from ?°°T', ee and RaC” would also be very 
useful. 

The predicted levels of RaE and ?°*Pb, in their large numbers, make it clear 
that, unless the present ideas are very far wrong, heavy nuclei possess, in general, 
many more levels than have been detected experimentally so far. The reason is, 
of course, that the primary excitation is very limited, and the subsequent cascade 
of gamma-rays can only explore a very limited region of the level diagram. 
In *8Pb, for instance, there are probably about a hundred levels between 
2 and 4 Mev of which only some six are known at present. 
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ABPRPENDIX 


THE LEVELS OF *’Pb 

Recent experimental work furnishes an almost certain level scheme for 2°7Pb. 
The crucial pieces of information are: 

(i) The measured spin and magnetic moment of the ground state, 
indicating pj). 

(ii) The short-range alpha-groups of *44Po(AcC’), indicating three excited 
levels at 0-54+0-04, 0-88 + —0-04 and 1-11 +0-06 Mev (Neumann and Perlman 
#951). 

(iii) The discovery by Campbell and Goodrich (1950) of an isomeric transition 
of energy about 1-05 Mev. Its lifetime, 0-9sec, places it unambiguously in the 
M4 class (Goldhaber and Sunyar 1951), in line with the expected is). —f5) 
transition in this neutron shell. 

(iv) The measurement of the gamma-rays following electron capture in 
207Bi (Neumann and Perlman 1951, Grace and Prescott 1951) of energies 
0-565, 1-063, 1-46, 2-05, 2-20, 2-33, 2-49 mev (and a soft one whose energy is 
either 0-137 or 0-064Mev). The 1-063 Mev radiation can be identified with the 
isomeric transition. It is in cascade with the 0-565 Mev gamma-ray. 

(v) The existence of a weak gamma-ray of energy 0-870 Mev in the spectrum 
of 2°7T1(AcC”) +*°7Pb(AcD) (Surugue 1946). 

(vi) The short-range groups in the reactions 7°*Pb(dp)?°’ Pb and 7° Pb(dt)?°7 Pb. 
Both show unresolved groups interpreted by Harvey (1951 a, b) as two levels at 
0-61 and 0-95 Mev, but which are also consistent with three levels at 0:54, 0-87 
and 1-11 Mev, and broad groups corresponding to levels in the region 2:0—2:8 Mev. 
The (dt) reaction indicates a level at 1-61 Mev and another at 2-33 Mev. The 
(dp) reactions reveals further levels at 3-60, 4-42, 4-66 and 5-28 Mev (incompletely 
resolved, and so possibly more numerous). 

The limited energy available in the decay of AcC” suggests that the 
0-870 Mev gamma-ray is a transition to the ground state of ?°’Pb, implying a level 
at 0-870 Mev. This is clearly the 0-88 + 0-04 mev level revealed by the alpha-group. 

The 1-063 Mev isomeric transition must be interpreted as ij3), f5j2, and so 
cannot be to ground. Its coincidence with the 0-565 Mev gamma-ray indicates 
that it is between levels at 1-628 and 0-565 Mev, the latter clearly being the level 
at 0-54 +0-04 Mev in the alpha-groups. The i,3). level is then at 1-63 Mev, in 
excellent agreement with the level observed in the (dt) reaction, and there are four 
experimentally established levels below it (revealed by the four alpha-groups 
of AcC’). Since it is long lived, the hg). level cannot be below it, for the 
M2 transition to it would take place much too quickly. At first sight it would 
also appear that the f,). level cannot be below it, since the E3 transition to that 
would be fairly rapid. But this would leave only the pj), pg and f;). levels to 
be fitted into the four experimentally established places, and the only way to find 
the fourth is by assigning the f,j. to the 1-11 Mev level. ‘This now fixes the 
P3j2 level at 0-870 Mev, and the f;). level at 0-565 Mev. ‘The energy difference 
between i,3/. and f,). is thus 0-52 + 0-06 Mev, and a comparison with the empirical 
results of Goldhaber and Sunyar shows that the partial lifetime could be anywhere 
between 0-01 and 100sec. The observed predominant decay of the i,5). state 
to f;). with a lifetime of 0-9 sec is therefore compatible with the proposed scheme. 
It may be remarked that the alternative hypotheses, that the 0-565 Mey transition, 

PROC. PHYS. SOC. LXV, ILO—-A 52 


790 M. H. L. Pryce 


and not the 1-063 Mev transition, is the lifetime determining one in the isomeric 
state cannot be made to fit the facts. 

This assignment of the 1-063 Mev gamma-ray in the level scheme conflicts. 
with the scheme proposed by Neumann and Perlman, who, however, were 
unaware of the isomeric nature of the transition, and assumed that it corresponded 
to the transition from the 1-11 Mev level to ground. 

The four gamma-rays of energy above 2Mev most probably correspond to 
ground state transitions, for, as pointed out by Neumann and Perlman, the 
estimated energy available in the decay of ®°’Bi is insufficient to lead to levels 
above 2:5 Mev. The 1:46 Mev gamma-ray can be fitted in within the experimental 
error as a transition between the 2-05 and 0-565 Mev levels, or between the 
2:33 and 0-870 Mev levels. 

According to Prescott (unpublished) not more than 10% of the transitions 
from the isomeric state go via the fy). level, which would give rise to a cascade of 
three gamma-rays all with about 0-5 Mev energy. This means that any possible 
confusion of these gamma-rays does not seriously affect the value deduced forthe 
energy of the f;). level, which we shall take to be 0-565 + 0-01 Mev. 

Prescott also finds that in coincidence with the K x-rays he observes a_ 
gamma-ray of energy approximately 0-55 Mev. ‘The absence of coincidences 
with the 1-063 Mev line confirms that this is an isomeric transition, while the 
other coincidences indicate that either direct capture into the 0-565 mevy state is 
taking place, and/or that it is formed by prompt gamma-transitions. Prescott 
also finds indications of a gamma-ray of energy 1-78+0-05 Mev in the total 
spectrum. ‘This could be a transition from the 2-33 to the 0-565 Mev level. 

The position of the hg), level is left undetermined. It is unlikely to be formed 
strongly from the direct decay of ?°’Bi, for considerable internal rearrangement 
is necessary, although the spin of the ?°’Bi is probably also hg. One would ~ 
expect the 7°°Pb(dt)?°’Pb reaction to give rise to the levels corresponding to one 
neutron hole, and the level at 2-33 Mev is therefore possibly the hyp. This cannot 
be the same as the 2:33 Mev level found by Neumann and Perlman, since this 
radiates to the ground state. 

The four levels above 2 Mev most probably arise from configurations in which 
at least one nucleon is excited to the next shell. 
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ABSTRACT. The angular distribution of y-radiation from nuclei aligned at low tempera- 
ture is given as a series of Legendre polynomials, [(0)=1-+-2,Ax(T) Px (cos 0) the Ag being 
temperature dependent. The formulation is applicable to cascade processes of arbitrary 
multipole orders and nuclear spins. An approximation valid for low degrees of alignment is 
given, as is the angular distribution from completely aligned nuclei. Application is made to 
S°Co; 


§1. INTRODUCTION 

F radioactive nuclei can be partially aligned at low temperatures, that is, 
] caused to occupy states such that the spins are preponderantly in one 

direction, an anisotropy can be expected in the emitted y-radiation. In 
the case of a single emission the theory of this effect has been given by Spiers 
(1949). In the following the angular distribution is given in the form of a series 
of Legendre polynomials, /(?)=1+2,A,;(T)Px(cos@) where @ is the angle 
between the direction of emission and the axis of alignment and T is the absolute 
temperature. The A, are temperature dependent and are such that the angular 
distribution from an arbitrary cascade can readily be given provided the decay 
scheme is such that et each stage the y-radiation is of one multipole order only, 
or if B-emission from the oriented nucleus precedes the y-emission (or emissions) 
the total angular momentum of the 6—neutrino system has one value only. 

The first two terms of the series give the polar diagram, (a) for arbitrary 
nuclear spins where no multipole orders higher than quadrupole are involved, 
(6) for arbitrary multipcle orders where no initial nuclear spins greater than 
5/2 are involved. 

An approximation valid for low degrees of alignment i is given, applicable to 
arbitrary nuclear spins and multipole orders. ‘The maximum possible value 
of anisotropy, for completely aligned nuclei can also be readily obtained. 

Alignment of nuclear spins can be accomplished by allowing a strong 
magnetic field H (~10° gauss) to act on the nuclear magnetic moment yp, at 
low temperatures (~0-01°K). Under such conditions, if the lifetime of the 
unstable nucleus of spin J is sufficiently long, the population of states of varying 
z-component of spin M assumes a Boltzmann distribution. If W(M) is the 
relative population of state WM 


f 
IV(M)=exp pM | x expfM, 
M=-J 


where 8=Hp/RTJ. ‘This method has been described as the ‘brute force” 


method of aligning. nuclei. 
If 7,,(0) is the angular distribution of y-radiation emitted by nuclei with 
s-component of spin M, the angular distribution from all such states is 


1(0) = e WM)Iy(8). med) 
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Bleaney (1951) has suggested a method whereby the magnetic field of the 
-atomic electrons is used to align nuclei. In this case the interaction is such that, 
in general, M is not a ‘good’ quantum number. However in special cases (e.g. 
60Co, Bleaney 1951) the present work is a good approximation. 

Pound (1949) has suggested a still further method of aligning nuclei in which 
the interaction of a crystalline electric field gradient with the nuclear quadrupole 
‘moment Q is utilized. In this case the population function is, in Pound’s 
‘notation 


J 
W(M)~expyM | x expyM? 
M=-J 
eNO 
4 J(2J—1)kT Oz 


assuming axial symmetry. The angular distribution is then given by eqn. (1) 
with this population function. 


where y=- 


§2. SINGLE EMISSION 
Consider a transition in which a nucleus of spin J decays to a nucleus of 
spin J, emitting a single y-ray of multipole order L, which can be represented 
by (FIA) where |J—J,| <L<J+J,. Spiers (1949) has shown that the 
distribution from a single initial state / can be written 


Ty() ==, |Ch, 4 M— m, Be (9) bia cea «3 (2) 
independently of the character (electric or magnetic) of the radiation. The 
C4, 4—m,~ are Wigner (Clebsch—-Gordon) coefficients which add states of 
angular momentum J, and L to form a single state of angular momentum J 
with z-component M. oe are given by Condon and Shortley (1935). In 
Condon’s notation C% 43m, #% nana M-—m,m|J,,L,J,M). We can write 


14)= & [AMmO)e, 


‘where A?”"(64) represents for c= +1, 0 the components of the magnetic vector 
potential of multipole order L in the form of a spherical wave. Solutions of 
Maxwell’s equations in this form were obtained by Heitler (1936) (see also 
Goertzel 1946). 

The following expansion of [4(@) is derived in Appendix I. 


K(K+1 
Tx(0)=1+ ¥ (2K +1) {1 " mah CEELCL EL Py (cos8) os... (3) 


which can be written 
T(0) = =1 sus Lx a,(L)U x(m, L)P x (cos 0). 
K takes only the even values 2,4,....2Z. ‘The Legendre polynomials Py (cos #) 


‘are normalized in the sense ai P;,(1)=1. The expressions II x(m, L) are 
essentially the coefficients Ema 6% omitting a normalizing factor, 


(—1)(L+m)!(L—m)!K! 
leon SS m—K-+v)! 


nd. a,;(L)= _ K(K+1)) {(2L+1)(2L—K)!K!? 
and e(L)=@K+1) {1 Est | ere 


IIz(m, L) = 


II ,(0, LY. 
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‘The first two IIx(m, L) are 

II,(m, L) =[3m?— L(L +1), 

TI ,(m, L) =3[35m4 —5(6L? + 6L —5)m?+3L(L+1)(L+2)(L—1)]. 
‘The first two a,(L) are , 


3 5 
all)=~ {1~ aEey| GERHOESD 


Ethin {1 = 10 81 
a aa L(L+1)) (2£+5)(2L43)(2L—1)(2L —3)° 
It is not difficult to show that 
ml Ch, tm, HP We(m, L)=Sr(J,J1, Lg MJ) «we (4) 

where Sx is independent of M, for Ch, 44-m, # Ch X% represents a decomposition: 
in which J is decomposed into J, and L, and L is further decomposed into K 
and LZ. An alternative is to decompose J into K and J’ and J’ into J; and L 
obtaining the same final states (this is a generalization of the simple alternatives. 
of LS and jj coupling in atomic spectra), therefore 


Jd LCL KL % L 
Car Moe pO Diy Wy Cie Bot Che, hme 


Multiplying both ae i Cir, —m,% and summing over m we obtain 
En | Cie, B—m, a PCr, Om = Da Wyo, ie 8y, «= WiC, 0, 
whence eqn. (4) follows (W, is determined by this equation). 


The S,(J,J,,L) in all cases of interest assume relatively simple forms. 
The following are derived in Appendix II: 


be _ (2£)1(24—K)! 
ee ge (2L—K)!(2J)!? 
J—L+1 
Sr(J,J —L+1,L)=S;,,J—-L, L) E - PK) 
— (-)A2L)'(2F+1)! 
Sx, +1 1)= Or eKy(aTeK 41’ 
iti 74 REE APR OR POO A 1)|1- ks oe 
Due to the relatively small probability of dipole transitions the following also 
may be needed: eee 3(27 +5)(2F —3) 
OSD a are aha way NED 
Dale! p42) = a 


J(J +1)(2F +3)(2F —1)° 
Thus we have 

Ty (0) =1 4+ Dxax(L)Se(J, J, L) MN g(M, J) Px (cos 8). 
It should be noted that Il,(M,/)=0 if K>2J. Thus the maximum value of K 
is restricted to K <2p where p is integral and satisfies p<J, p< L. Thus if J =0 
or 4 no angular effect at all can be obtained, if J=1 or 3/2 only the term in 
P, (cos @) is required regardless of L, if J =2 or 5/2 only the first two terms are 
required and so forth. 
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The complete polar diagram is given by 


1(0) =DyW(M)Iy(0) =1+E,Ax(T)Px (cos 9), , ries) 
where Al T) Ghd) Set phgpl ene ce 
B,(T)=2yW(M)IR(M, J). a 


If W(M) =exp BM/XyexpBM the summation over M can be carried out 
explicitly since 


J com _ Sinh B(2J + 1)/2 
aes ri sinh 8/2 ; 
a? d® [sinh B(2J + 1)/2 
(wey i ee A a in scat Nee 
aa cars ral sinh B/2 


and so forth, so that 
YyW(M)II,(M, J) =2)(J + 1) — 3[(2J + 1)ab—b?]=B,(T), 
Ly W( MI, M, J) =3 {8J (J + 1)(J +2) — 1) — $(27 + 1)(8J? + 8J — 9)ab 
+ 2(4J? +4 —1)b? —*2 (27 + 1)ab? +1254} =B,(T), 
where Poiheak B(2J +1)/2, b=cothB/2. If more terms are required, or W(M) 


is more complex, it would appear that for calculating purposes the angular _ 
distribution is better left in the form of eqn. (5). 


§3. APPROXIMATIONS* 
(1) For Low Degrees of Alignment 
II,(M, J) involves only even powers of M. This implies, since J,(6) =1_y(8), 
that it is immaterial whether nuclei are aligned in both sense and direction, 
i.e. all spins pointing in the +2 direction, or aligned in direction only, spins 
divided between the +z and —2z directions. Thus, for orientation by a magnetic 


field ebM h 
1(8) = © srep} Fal) = 2 so ag} Ll) = = W (My) 


For smal] values of 8 we can expand cosh8M and take only the term in [? 
obtaining 


1 2 1 2 
W"(M) = We arn! = F3Me—J(d + iyi} — arent! =e -i¢M,J)} 5 
Since Lyllz(M,J)=0, and Lyll,(M,J)0,(M,J)=0 unless K=2 while 
Yylll,(M, J) |? = J (J + 1)(2J —1)(27 + 1)(27 +3) it follows from eqn. (5) that 
we can write 
Ty_, (8) =1 + $B%a(L)h(J,J,)4(3 cos?8—1), sw ee we (7) 
where h(J,J,)=(2J+3)(J+1) if J,=J—L, h(J,J,)=J(2J—1) if J,=J+L, 
o(L)=(3—L(L+1)]/(L+1)(2L+3). The same expressions hold for Pound’s 
method if we replace 6? by 2y. 
An exact criterion for the validity of this approximation would be difficult 
to give but some investigation shows that it would be ‘good’ for 6? <1/J(J +1). 
If « is defined as the anisotropy « =[J(7/2—J(0)]/I(7/2), this can also be 
approximated for small f 
__=PBt( LMI Ss) ap, 
€p>0= T— 4 ?0(L)AV, J) = — GPP a(L)h,J,) wna a es (8) 
for 6? sufficiently small. 


* Some of the results given here have appeared previously (Steenberg 1951). 


Angular Distribution of y-Radiation from Aligned Nuclei 795 


(ii) Completely Oriented Nuclei 
For very large values of f nearly all the nuclei will be in the state with M=J, 
i.e. saturation. This occurs practically for 8 greater than 4. Then e* is much 
greater than e™ for MJ and I;_, ..(9)=1,(8) and since 


| Pee 
II (J, J)= K!(2J—K)! (K even) 
i (0)=14+ Sa,(L)SzVJ,J7 fy cele Bay (cos 6) 
B—> = K fs Tan K!(2J—K)! K ° 
The value of e« at saturation js 
(2J)! (—1)k? K! 
Bee ete TOF RI ORD) =] 


€3-> © 


2J3)! ee el 
ig UO CBS?) MIea a Soe 


A particular simplification occurs if J; =J—L for then 
158) =Eml CF, 3= Im, ml? Lr) = mB mtZn( 2) 

and J, , ,.(0) =1,(0) =I4(6). The I¢(@) are given by Arnold (1950). In Arnold’s 
notation [4(0)=F4(6). The following are taken from Arnold’s results (and 
normalized) : 

I,?(8) = 3(1 —cos* 6), 

I,3(0) ==22 (1 — cos? 6 — cos* @ + cos® 4), 

I,4(0) = 8 (1 —2 cos? 6 + 2 cos* 6 — cos® 6), 

1;5(0) = 422 (1 — 3 cos? 4 + 2 cost @ + 2. cos* @ — 3 cos® @ + cos!” 4). 


§4. MULTIPLE EMISSION 


Consider a decay scheme represented by (J. oy me — 3) where /;, Jo etc. 
are the total angular momenta of emitted particles 1, 2, etc. If, as will subsequently 
be assumed, particles 2 and 3 are photons we write L,, L, in place of j, and J. 
It should be noted that the half lives of the intermediate nuclei of spins J, and 
J. must be sufficiently short that no re-alignment takes place. We consider here 
a three stage cascade but the formulation extends in an obvious way to a cascade 
of any number of particles. 

Assuming particle 3 is a photon, the angular distribution of all photons, 

L, (=j3\, coming from a single initial state M is 

IGX(9) == vat m,| Cis, Mi, e—u, Cut, th, tt, CU, M—ms, mil Tmi(O). ovens (10) 
This result can be siteaiea by a simple extension of the hae used by Spiers 
(1949) to obtain eqn. (2). 

I#:(0) is expanded in the form of eqn. (3), as for single emission, 
T}3(8) =1 + Eg ag(Le) 1 g(ms, Ly) Px (cos 0) 

and the summations indicated in eqn. (10) can all be performed using the properties 
of the II, expressed in eqn. (4). ‘Thus eqn. (10) reduces simply to 

19(0)=1 + Ege Lo) Se(Joy Ios La)Se vd J2)SxIniv I) x MI)Px (008 6). 


K is now restricted to K <2p, where p, is integral and satisfies p; < Ls, p, <J, 
Pe <Ji, P3 <Jp. 
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If particle 2 is a photon as well the same procedure is applied (L. =/2), giving 


2D. i 
1Q(O)=14 E ap (Le)Su(aoIn Le) SaI dave (Mo T)P x (6080), «++ -(12) 
ke 


where p, is integral and satisfies p. <J,J,, Lo. ”. 

If particle 1 is a photon its angular distribution will be calculated as if it 
were a single emission. If it is a B-particle which is not observed eqns. (11) and 
(12) still hold and j, is the total angular momentum of the B-neutrino system. 
Assuming particle 1 is not observed, the observed angular distribution of photons 
2 and 3 is 1(6)=S),W(M)[8_19)(6) + 8,19)(0)] where 5, and 6, are detection 
efficiencies. If W(M)=exp8M/ZexpfM, and it is so desired, the summation 
over M can be performed giving 


1(8) =14+ Xx [8p AP(T) + 85’ AP(T)]P x (cos 9), 
where 5) = 50. o.oo er nen: 


This formula can be approximated for small values of 8 by a method analogous 
to that used for single emission, giving . 


see eee 


a(L) and h are defined in § 3(i), and 
J,(2J,—1) 

Jo \e eee 
Fa 1 2) J(2J,—1) 
UI,Is)=8 Ai) if y= |F—-Si |, 


(J, +1)(2J, +3) 


ed at a Cee Va Rats WV SELL | 


if Sp dy ae 


POE Dis sce ae ate 
= ye eat) een ore 

© Gi 3) Pad l= sto |) ee 
USI) og 3) Sahel yee if? 7,=3 —J, 41. 


For the saturation angular distribution the same considerations apply as in 
the case of single emission, namely J,_, .(@)=1,(@). In particular if 7, =J —J,, 
L,=J,—J,, L,=)n—7),, 

1(8) = 821749) + 331739). 


§5. APPLICATION TO “Co 


This procedure can be illustrated by its application to Co. This nucleus 
has been successfully aligned recently (Daniels, Grace and Robinson 1951) 
by Bleaney’s method (Gorter et al. 1951 also report success with this nucleus). 
It is known that the y-cascade in ®°Ni is (4a OO) both transitions electric 
quadrupole (Brady and Deutsch 1950). Recent work (Deutsch and Scharff- 
Goldhaber 1951) suggests that the ground state of Co has spin 5. We therefore 


“8 1 2(y 2 
assume a transition CRG SSE Ye Now from eqn. (3) 


Ti, (9) =I, (0) =1 — 23m? — 6) P,(cos 0) — 4 (35m* — 155m? + 72)P,(cos 0). 
Since in this case L,=J,—J, and L,=J,—J, the product of the factors Sx in 
ISH(0) and I92(0) are identical and from eqns. (10) and (11) we get 


SOME ce td be 5 1H 12. 11,(M,J)P,(cos6) 
Pxt(8) = 2548) = Lu) =1 — 5 Fey z aay Palos 9) — > Fe Tyo 3 aT) 
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and /(@) =(52+83)2y,W(M)Iy(9); the detection efficiencies appear as a common 
factor and can be disregarded. Finally, with J =5 
1(8) =Xy,W(M)Iy,(0) = 1 —$(40 — 11ab +b) P,(cos 6) 
—zga(128 — 121ab + 516? — 11ab + b*) P,(cos 4) 


a=coth118/2, b=coth$f. 
uA plot of the anisotropy, « =[I(7/2 —1(0)]/I(7/2) derived from this function 
is given in the figure, curve (a). Also included are curves representing the 


a ee ao ee SS LT 
co 19-0 95 475 TxI0°(°k) for(a)and(d) 


Variation of anisotropy, ¢=[J(m/2)—I(0)]/I(7/2), with B. (a) (; ——»4—>2—+> 0), 
6) 6—$— 14 —$_ 2.0), (0) 4 20), @) (6 — 4 2 0). 
The temperature scale refers to curves (a) and (d) only. 


decay schemes (eg aay) curve (bd), (ot curve (c) and 


et 9 >) curve (d). A temperature scale is included for curves (a) and 
(d), J =5, on the assumption that the magnetic moment of ®Co is 3-00 nuclear 
magnetons. ‘The effective value of H is known from the hyperfine splitting in 
59Co (Bleaney 1951) and the magnetic moment of °®Co (Mack 1950). 

The figure of course represents uniaxial alignment. In the experiment 
of Daniels et al. (1951) a crystal was used which gave biaxial alignment, the two 
axes intersecting at an angle of approximately 76°. Intensity measurements 
were made in the meridian plane defined by the bisector of the axis and their 
common perpendicular. ‘This arrangement reduces the anisotropy at saturation 
from 100%, to 38-6°/, (cos*38°) in the case of the (5 >4—>2—>0) transition. 
It should be noted that in this case since Bleaney’s method was used M is not 
strictly a ‘good’ quantum number so that the theoretical curves given here 
require correction which is expected to be slight. 


§6. DISCUSSION 
Certain features of the angular distribution are common to all decay schemes. 
First if J =0 or $ no angular effect can be observed for any emission, if J, =0 or $ 
no angular effect can be observed from a second particle and so on. 
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Secondly, if an angular effect can be observed, from eqn. (8) it can be seen 
that for small values of f, « varies with f? or T-*, the slope of the curve of « 


against 7-* being 

de(T) 1 (Hy ; 

are we nea 289! a Len)Cn 
where C,, is the product of the appropriate constants, and the summation is 
over all observed y-rays. “a only one y-ray is observed, and if it is a dipole, it can 
be seen that since «(1)=;4 and «<0 there will be an excess of radiation along 
the axis of alignment. Fine all higher multipoles «(L)<0 so that «>0 and the 
excess is in the equatorial plane. Where two or more y-rays are observed there 
will always be an excess in the equatorial plane as long as dipoles are not involved. 
If dipoles are involved no general conclusions can be drawn. 

Thirdly, as 8 is increased (T reduced) « approaches its saturation value. 
If dipoles are not ‘involved this is always greater than zero. If furthermore, 
J, ,—J,=J» for all stages of the cascade then J, _, (0) =0 as can be seen from the 
Ti(8) listed in §3, whence ¢;,,.=1. If this condition is not met, (e.g. in such 
decay schemes as (> and 54 the anisotropy 
at saturation will be less than 1. 

If dipoles are involved it may be negative. 

Finally it should be noted that while the above approach has been applied 
only to decay schemes in which pure radiations are involved it would appear 
difficult to include forbidden transitions* and other complex decay schemes 
such as branching transitions and mixed multipole transitions. ‘These would 
then involve nuclear matrix elements. 
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APPENDIX 
THE EXPANSION OF J4(@) 
Al"(6¢) can be written (see Goertzel abe: 


AM(0$) = CE ce °, ; Ur)YE °(¢) 


where Y7°(6) is a normalized spherical harmonic and f;(r). is a function of 7 
which can be disregarded. This can in turn be written 


A$) =L_ oY 799), 


where. L,= “Rb +7L,)[2L(L+1)), Lo=L,/|L(L+1)7? and L,, L,, & 
oe the differential _operators representing angular momentum. "unions 
AM (0b) = —(—1)°L,Y%(0¢). Consider the Ne LA” 1m (0jb;)A(6;4,) 
which reduces to &,| A?(8) |? when 6;,¢, is put equal to 6, 4. Now 
3, Aba) Alm) = —¥,(—1'L QLPV DP A()). 
* Allowed £-transitions under all selection rules involve only one value of the angular momentum 


j of the B-v system. Forbidden transitions may involve as many as three different values of 7. 
See for example Blin-Stoyle and Spiers (1951). 
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But Y%(j)=(— 1yny7(7) and the fone Y,™()Y%(2) can be written 
Yr") Y5)= i ee eRe 7) veeees (A1) 
K=0 


where F is an eigenfunction of (L®)?,(L)? and of the differential operator K? 
which corresponds to an angular momentum K which is the vector sum of 
L‘) and L®. The eigenvalues are L(L®+1), LO(L®+1) and K(K+1) 
respectively. Operating with 

LO.L® 1 
Dial) 2141) 
ad. puttiac & $;= aS we obtain 


On comparing eqn. (a1) with an expansion deduced by Blin-Stoyle (1951), 


: 21 (2L +1) 
A iE Do ORL L ~sh aecisalhce! ava 
Ws (68) ¥'x(00) = ae Co, m, —m Co, 0,0 Vx(4) [47(2K + Die 


2 (~ Wis 1)°L® PAGS a [K2— L@2— 12] 


| . eels | 
we can write FRO\— Cre 0 Y&(9) ee tam : 


oa 21) 
Noting that Ce ae lye Con, 0, m ( + i) 


the final result is, dropping a factor 1/47, 


K(K+1) 


21 
TX(6) = im, a = — 
m(9) ae (9) | 1+ 3 {1 21(E 221) 


haK+ 1) ar Gen Gan Lem COS t). 


APRPENDLXs I] 


The scalar factors S;(J,/J—L, L) and S;(J,J —L +1, L) are readily derived. 
It has been shown that 


x mtCw % M— —m, m me I x(m, L)= SK(S,A1, DL) g( M, J). 
(a) J,=J—L. 


If we put M=J the Wigner coefficients reduce to [C7 y— KP? =6,,,, and 
therefore 
ULL) (2L)'(2J —K)! 
| na ae ae Tis (2L RL 
(6b) J,=J-L+1. 
For M=J, m takes only the values L and L—1 whence 
S,(J,J—L+1, L)=[UgVJ, J) {| CF 5277 FP Wg(L, L) 
+ | C3, 5-757, 7-1 P MeL —1, L)} 
(J —-L+1) 
ds RRR Ke ree 
= S,(J,J—L, L) E £K(K +1) Lanal: 
It is evident how this procedure would extend to S;(J,J—L+2, L) and so 
forth. | 


The remaining cases are considerably more involved, and the derivation will 
be merely sketched. 
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(c) J, =F EL. 
II (mL) is rewritten in the equivalent form 
wi —1)+*?(L+m)!(L—m)!K! 


wep GK +H) 'GK mH) IL — 4K +m—p)!(L— 9K —m tp)! 
Then by writing | Cj, 37%, 2? explicitly, interchanging the order of summation, 
and some manipulation.of factors the following result is obtained : 


Ln | Cir, Mm, a |? Ue(m, L) 

(—1)"*42K!(J'+ M—p)!(J’- M+yp)! o CU, Ith Ui fe 
PEK FRUGK HIRT + MYT ME nage! Cie | 
where J’=J+4K, M’=M-yp, L'=L-iK, A=m-np, and N;; Py are the 
normalizing coefficients of | C4, 4°, 2? and | Cir, 7,7 |? respectively. 

_ (27 +1)! (2L)! p -@I+K+1)\(2L—K)! 
1 (2F4+2L+41)!" at (2J+2L+1)! : 
The summation over A is unity. Putting M=J and rewriting the summation. 
over y in an equivalent form, we are left with the result 
x (-1)(2J +v)!(K—»)!K! 
=H ees WEG = 
Nee a [((K—v)!v!](2J)! 

We now note that II;(J,.7) =(—1)*(2J)!/K!(2J —K)! is the coefficient of «* 
in the binomial expansion of (1+x)~@7-*t», By rewriting (1+2)-~7—** as 
(1+x)~%+(1+4x%)* and expanding the latter expression, the coefficient of 
x* is. also found to be 


=N,P> 


=Sitd JEL, Elen 


K (—1(2J +»)! K! 
0 (Kv) Lv QI) 


ed (—1)4(2L)!(27 +1)! 
Whence -Sa(J.J FL, 1)—(— 1) Mili = open 
(dZ) J,=J+L-1. 
Sz(J,J + L—1, L) is found by a procedure similar to (c) above. 
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I: Selection Rules and Energy Spectra 
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ABSTRACT. The matrix elements of beta-decay theory are expressed in a way which not 
only gives the beta-spectra and angular distributions for arbitrary degree of forbiddenness, 
but also yields very simply the selection rules appropriate to any type of interaction. 
Results for the energy spectra confirm those found by Greuling by inspection ; angular 
distributions and electron—neutrino correlations are considered in a subsequent paper. 


§1. INTRODUCTION 

N the theory of beta-decay an investigation is made of the matrix elements 
] associated with nuclear transitions which result in the emission of an 

electron (or positron) and a neutrino. ‘The properties of these matrix 
elements lead to the classification of beta-emissions by ‘ degrees of forbiddenness’ 
and to the corresponding spin and parity selection rules, energy spectra and 
angular distribution properties. 

The states of the emitted particles are expressed as plane waves or spherical 
waves (Dirac central field solutions) and the selection rules, etc., are then 
usually found by expanding these wave functions in powers of the coordinates 
near the origin (nucleus). From these expressions vectors and tensors are 
constructed whose complexity increases with each degree of forbiddenness 
(Marshak 1942, Konopinski 1943). 

It is the purpose of this paper to give a method of writing the matrix elements 
which avoids this process by using eigenstates of the total angular momentum 


J=0,1,2,.... ofthe system ‘electron plus neutrino’ treated as a whole. Matrix 
elements so written are then expressed as a sum of terms specified by quantities 
L=0,1,2,....and S=0, 1 which play the part respectively of total orbital and 


total spin angular momentum of the system. 

Consideration of these terms (matrix elements for the creation of the system 
electron plus neutrino in a state specified by J, L, S) leads to simple rules for 
the degree of forbiddenness of a beta-emission between nuclear states of given 
spin and parity and to expressions for the energy spectrum and angular 
distribution properties of the emitted particles for arbitrary degree of 
forbiddenness. 

To anticipate, let a particular interaction (Fermi, tensor, etc.) be chosen; 
then for each Dirac operator of that interaction the non-vanishing matrix terms 
can be picked out and classified by degree of forbiddenness as follows : 

(i) If the nuclear spin change is J;>J,, the possible values of J are |J;—//,| 
to J;+ J; in integral steps. 

(ii) For the scalar operators 1, f, y;, Bys, S=0 only; for the vector operators 
a, Ba, 6, Bo, S=1 only (electron and neutrino emitted with spins antiparallel 


and parallel respectively). 
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(iii) For given J, L=|J—S| to J +S in integral steps. 

(iv) The degree of forbiddenness f of a transition creating the system electron 
plus neutrino in the state JS is numerically equal to L for the operators 
1, B, 6, Bo, while f=L+1 for the ‘velocity dependent’ operators a, Ba, ys, Bys 
(0 =allowed, 1 =first forbidden, etc.). 

(v) If there is (is not) a nuclear parity change on emission, only states of the 
system electron plus neutrino with odd (even) f can be created. 

Table 1 is constructed using rules (ii) to (v). Using rule (i) and selecting the 
Dirac operators of the chosen interaction we can enumerate the f-values of all 


Table 1 
Interaction Operators S J=0 1 2 3 4 
I=, il OSA ee. tyes 2 Sik: 3, 4, § 

Scalar : | 

1 DaabP 0 fe0, #) tect Ne eae eg, ae 
Fermi Sa 

ears Ot, Ba 1 f io 2 AG 2; (3) te ae (4) 3: 2 (5) Gar EF (6) ' 
Tensor | 

as 6, bo Lo fC 4 Gl kz) 1, 2543) re Ce 3, ~4; (ap | 
Axial vector oe / 

Been Ys Bys 0 f=1, * oe D2 * ee 3 * ¥ 2 * a 5; * 
Pseudo-scalar ode 


Degrees of forbiddenness f of transitions creating the system electron—neutrino in states 
JSL. Even values refer to parity change ‘No,’ odd values (underlined) to parity 
change ‘Yes.’ * denotes that the term vanishes, ( ) that the term represents a small 
correction to other terms. 


terms contributing to the beta-emission ; the least /-value among these gives the 
degree of forbiddenness F of the emission as a whole.t ‘Thus, for example, 
for a nuclear spin change 0->+0, J =0 only, and table 1 shows that on the Fermi 
interaction, /=0 (allowed) for parity change ‘ No’ and completely forbidden for 
‘Yes’.. For 11, J=0, 1, 2,and_ on the Fermi interaction F—0 (( No: ane 
ice esi): 

§2. FORMAL DEVELOPMENT 


For convenience in considering angular momentum relations we write the 
matrix elements of beta-decay as 


M= | (PFE): Qd)* del 0k te ee (1) 
where ‘’,, ‘I’; are initial and final nuclear wave functions, #, ¢ are wave functions 
of the emitted electron and neutrino, K is a Dirac operator and Q=(KD*)*, 
D being the operator defined by Furry (1937, 1938). 

Spherical wave solutions of Dirac’s equation for a central field have been 
conveniently listed by Rose (1937). Inspection of these solutions shows that the 
wave function of a state of total angular momentum, m may be written in the form 


bima(toP) a Ca BFS aa 21) F449, p). Sree seashe (2) 
Here the C’s are Wigner coefficients used in the theory of addition of 
quantized angular momenta (Racah 1942) and the spherical harmonics are defined 


with the Condon—Shortley (1935) choice of phase factors. ‘The four Dirac 
components of each solution are characterized by the variables 8B = +1 (‘small’ 


t+ Except that for the pure pseudo-scalar interaction By;, F is counted as one unit less owing to 
the absence of any terms with f=0. 
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and ‘large’ components respectively) and o= +4; a takes the values +1, 
corresponding to the two types of solution (Rose’s ‘type a’ and ‘type b’ 
respectively) ; p is the momentum of the electron. For r equal to or greater than 
nuclear radius, the radial functions / become coulomb field functions; the 
relation to Rose's notations #,>t—7/,.,, =e, FS =i, Bret =2_.: 

Note that for aren values of a and #, (2) has the form of the wave function 
of a particle of spin } and orbital angular momentum j + 3af combined to give 
a total angular momentum j with z-component m. 

The neutrino wave functions will have the same form as (2), but with radial 
functions G appropriate to the case of Z=0 and negligible mass, thus: 


for pr<1, 
(Wai) FP +\(Z=0)=G+! = ah : (ae (3) 
(7a) Frh#'(Z=0)=Grh#! = Cae aoa (4) 
where (Qie tts Sh0C 1) = oo rey. (5) 


Here / is electron or neutrino momentum; to avoid confusion the latter will 
subsequently be denoted by gq. 


§3. THE DIRAC OPERATORS 
The operators Q@=(KD*)* will act on the variables o and f in the wave 
functions of the electron and of the neutrino, and must therefore be written as 
outer products of 2 x 2 matrices with matrix elements denoted by Qf. This is 
done in table 2, wheret 


cm): ves G8) (Ne ot 


Y 01 01 ‘ 10 
w= -iv2(5 A w= i(; le w= -iy2(, ab Sieh (7) 
Table 2 
Main term Velocity term 
Interaction pose Stes eS poe ee ae ee 
K a a wee as Oss 
Scalar B —1(6 5) pp(Fy)”” = aaa 
Fermi 1 + (Fy) pp oy)?” a +-1(6,) gpr(2o 4)° 
Tensor fo +1(6 ») ppm)?” Ba —1(8)gar(zw 4) 
Axial vector o — (oy) pp"(w 4)? Ys —1(0,)pB'(Fy)° 
Pseudo-scalar — = By 5 —i(8)ge(o,)°" 


Note that, for any of the interactions listed, we can factorize Q thus 
(Qasr =(Q)ar(Qy, nee (8) 
the first part of the product being independent of A. 


+ We use the notation Az= A,, Ay, A_, for the. components of any vector A, where 
A,=—(AztiAy)/V/2, Av=Az, Aa =(Az—iAy)//2. In this notation the scalar product 


A> Ba 4/3270, 0 AB 4. 


804 J. A. Spiers and R. F. Blin-Stoyle 


§4, STATES OF THE SYSTEM ELECTRON PLUS NEUTRINO 


By the usual theorem of vector addition, we construct states of the system 
electron plus neutrino with total angular momentum quantum numbers//, M, thus 


byujaltoP, r ‘o’B’) = 2% CHP imabiim'a! staat (9) 


where primes are used to distinguish neutrino variables from electron variables, 


and ja stands for jj’aa’. 
The matrix element for creation of the system in such a state is then 


MES SA73CO8 (PR EO emda ee (10) 


. 


where summation over o, o’, 8, 6’, X is implied; the suffix A and the expression 
— +/3Ctil , are to be omitted if K, Q are scalar. 
Now, using (2) and (9), 


bsse( POP, ro’) =D CY CB CY eV VP vt Beep oa ee (11) 
m, mm’ 


Mmm 
where _/, l’ stand for the quantities 7+ a8, 7’+3a’6’; and by a theorem proved 
in the Appendix this can be written as 
= ay (JSLS\CUS ssChL YY =F VG ae Gee (12) 


with L=/+lI' to |/—I'| in integral steps, S=0 or 1, 2=o+o’, and the A’s are 
constants. 
Inserting (12) in (10), we obtain, using the expression for 0% given in table 2, 


M\M i ja 
Myitja = = M. anes) = 2 | PST Re | ieee See (13) 
where RY, ¢57)= pe ee LS) FRGElO, seas (14) 


and for the scalar operators K =1, f, y5, Bys, 


Typ? = 1/27" YT (06)K J 
Ti =0 sisteleiers (15) 


while for the vector operators K=a, Ba, o, Bo 


TIP =0 
TIP = /2sCgt ye oe 
A 


It is evident from the form of T¥/* that it transforms under a spatial rotation 
according to the D’ group; accordingly the matrix elements vanish unless 
J =|J,—J;,| to J;+J; (rule (i) of §1) and can be written 

M TMGAL8) a Cob u M vate) aoe aa ene (17) 
where M)j,47,s) is independent of M, M,, M,. 

The factors r” in T and r~” in R are introduced in order that TJS may be 
equivalent (except for normalizing factors) to the 2 + 1 independent components 
(for M=—J, —J+1,...., +J) of the various Cartesian vectors and tensors 
(linear in K, homogeneous of order L in 7, and transforming according to the 
D’ group) used in the standard treatment. 
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The ‘average over a sphere’ (Greuling 1942) of T¥?5 is 


>| tS a ae pK)? (scalar operator) 
M ee 1% 2a 
_ 2 = a ele A ae Pe, eo ae vital (18) 
p™"4(K)? (vector operator) 


§5. PARITY SELECTION RULE 


From the form of 74% and the invariance of the matrix terms (13) under 
reflection of axes, it follows that if there is (is not) a nuclear parity change on 
emission, then Z must be odd (even) for the operators 1, 6, 6, Bo, and even (odd) 
for a, Pa, y;, By;. Since f=L for the former and L +1 for the latter (see below), 


tule (v) of §1 follows. 


§6. THE ENERGY SPECTRUM AND DEGREE OF FORBIDDENNESS 


From standard time-dependent perturbation theory, the probability per 
unit time of ejection of an electron with enetgy between FE and E+ dE is 


2E 
=k eee 19 
? (19) 


Mi M¢ 
= 2M yarja(ts) 


2 
P(E) dE = 2g" <MyMps>—> ; 
uw guia 27; +1 
gis the Fermi constant. &, denotes summation over the matrix elements for the 
various Dirac operators of the chosen interaction, with the appropriate value 
of S=0, 1. 
In order to obtain relative orders of magnitude (‘degrees of forbiddenness’) 
of transitions between initial and final nuclear states of various spins and parities, 
and to obtain explicit expressions for the energy spectrum P(E), we write 


M Fists) =| i YT id |Rit,s( P, Cat ee Ash eons (20) 
=e ie ua Myrs|Ri13(P aD gad Bae (21) 


where p is the nuclear radius, and the ‘nuclear matrix elements’ M,;, are 
independent of p, g, Mi, M;, M, j, 7’, a, a’; in estimating orders of magnitude, 
they are assumed to have roughly equal values for the operators 1, 8, 6, Po, and 
values of order v/c smaller than these for a, Pa, y;, Py;, the latter being classified 
as one degree of forbiddenness higher than the former. 

Substituting in (19) and using properties of Wigner coefficients, we have 


2g J+8 PBR 
FE eo |S Cota MaRS tee Levee (22) 
T Jia| K L=|J—8| Pp 


For the sake of brevity we ignore cross-terms resulting from different K 
and L values which may occur in transitions such that more than one state of 
the electron—neutrino system exists with the same J and same (lowest) degree of 


forbiddenness f (table 1). 


Substituting (14) and using the form of G for gp <1, we obtain after some 
algebra 


202 a. ome 
Pjps(£) = =| Myr" . 2Epq? <5 ite Fe Pe d 


=) 


pe le Pins a ue 157 pa? Oat il ach eae (23) 
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where* 
fO=LFRAE+ Fr Pepe) > ph nt Sy 
Boe sin Doe tek a ies hd Ep t(2v+1 HN 2v+31Noe 
fO FE a Pe ares aie cee (74) 
and 
ays = 2 (AER THA) PL(2L —2v — 1D 0 Da 

de 


BIS —2(Att Eos ly Aneey ,_s)(2L—2v+1)!!(2L—2v—1)!!}4 +... (25) 
yS— SD (Aner e_,_y)'l(2L—2v+ 1)! 7°. 
a=+1 


The negative sign applies for 1, 6, a, y;, and the positive for B, Bo, Ba, By,. 
Except for this sign the expressions for the energy spectrum depend only onJ, L, S 
(through the A coefficients). In obtaining these expressions we have ignored all 
terms which in the Z=0 approximation would lead to combined powers of pp, 
gp (<1) greater than the lowest (which is 2L, since A#(J LS) =0 unless 1+’ > L. 
It follows that the degree of forbiddenness depends directly on L, rule (iv) of § 1). 

Expansions in powers of «Z of the functions f{?, f, f®, are given by 
Greuling (1942), the relation to Greuling’s functions M,, N,, L, and Fy being 
(Psi M, flO N, f= Lol. 

Using the values for the A coefficients given in the Appendix, we obtain 

1. For the scalar operators (S=0, J = L), 


(L—v)(2v+1)!22-2" — ) 

ee, ) 
(L—v)(2v-+1)!24-% 

Us (lac eae ( 

(v+ 1)(2v)!22—-2" 


I Syd BEM ss 
a, =U, 


eee 


Lt O pe Ne Saad) NES ee 
i Uz (2L—2v+1)!(v!p2° 
_ (LNP24 
Us= OLLI ey ee we eee (27) 


2. For the vector operators (S=1,J=Z+1 or L),t 
obt1,L,1 _ gh+1,L,1_¢ 


(Qi 1 loi 9 oath ae eee (28) 


y= Un OI 4 Tyo 
act Beha ot ay, | 
ght o — BLO: nxn ee ay eet (29) 
yo = (be). | | 


‘These results agree with those found by Greuling (1942) by inspection, except 
for the factors U; and (L+1)/Z. The origin of these factors can be seen if we 
* The — indicates the values in the limit Z=0. 


_t Expressions for J=L—1 are only required in the special case J=0, L=1 with tensor and 
axial vector interaction (see table 1) and so are not considered here. 
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compare the ‘averages over a sphere’ (18) of the 7¥/% with those calculated by 
Greuling for the corresponding Cartesian vectors and tensors, and thus obtain 
the relations 


Kr,r) |? 
87?U;| Myzo ley = ae 
‘ K, r) |? 
87° Uz| Mr+1,2,1lev = fata Sethi (30) 
L+1 Kxr,r) |? 
8772 I U;| Miz |2, = eae 


§7. ALLOWED SPECTRUM 


Putting L=0 (hence J = S=0 or 1 according to the interaction), we obtain, 
pce ge =p, = Uys = U, = 1/47, 


2 
Poys(E) = $5 20 | Myos|*EP9°f0, 


which is the energy spectrum for all allowed transitions and for certain first 
forbidden transitions with L=0. /f%) is the usual correction factor for Z40. 


§8. CONCLUSIONS 
We have shown that by using eigenstates of total angular momentum of the 
system electron—neutrino, the selection rules and energy spectra of beta-emissions 
between nuclear states of arbitrary spin and parity can be derived. 
In a subsequent paper (Spiers and Blin-Stoyle 1952) the method used above 
will be applied to the derivation of angular distributions and electron—neutrino 
correlations for arbitrary degree of forbiddenness. 


APPENDIX 
A COEFFICIENTS 

It is well known from the theory of atomic spectra (Condon and Shortley 1935) 
that a state of two particles of spin 4 characterized in the usual notation by 
quantum numbers J Mj,j, can also be written as a linear superposition of states 
characterized by quantum numbers J MLS. 

Formally, this implies that the following relation exists between Wigner 
coefficients : 


MmyMms~ M449, ~ Mofl2F, 


ee ete ITS) CCI CN: (31) 


where pp =/4, +g, 6 =0, +0, and the constants Bi? (JLS) are to be determined. 
Bearing in mind the following relation between normalized spherical 

harmonics 

(21, + 1)(22,+ 1) 


1/2 
= ch ypyy = [Cees ] CHEY (32) 


Hye 1 
MiP2 


we can then write 
Sy Cec Cizlat Ypyi 


Myf FO, ~ MoflgFn 


Hips ane = 
a = JA OTBRI EH Ee AG os a alla (33) 
ie (21, +1)(244+1) ]¥? Avra, put, 
where Ab(J LS) = Tee (ce TENS) a anion (34) 
i 53-2 
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In the theory of beta-decay, only those values of /, and /, are of interest for 
which J,+/,=L, and it is then easy to calculate the B pation by solving 


eqns. (31). 
The following results are obtained: 
For S=0, J=L, 
Bhh - cide Bre ces. ae 
L+tlh+t (21, as (CARE WO) METS) ae wk h-th+t— Ze | (35) 
eee -lata] Bi? yy,-4=0- 
For S=1,J=L+1, Be4,,,,=1. All other coefficients vanish. 
Fors =1, f=, 
Bile Se iB = - ‘Bhi tl (L414 + 1) 
+4143 [L(21, + 1)(2d, + 1)]¥? 4—$l,+4 L(2l, +1) mall yy: 
(E+ Diep ee ee (36) 
Bayport Br 1,-4=0. 
ali L(21, +1) i j 
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ABSTRACT. The method of a previous paper is used to derive general formulae (valid 
for arbitrary Z) for the angular distributions of beta-emissions. 


§1. INTRODUCTION 

N a previous paper (Spiers and Blin-Stoyle 1952, to be referred to as I), by 
lf using eigenstates of the total angular momentum of the electron—neutrino 

system, expressions are derived for the energy spectra of beta-emissions 
between nuclear states of arbitrary spin and parity. 

The object of the present paper is to apply this method to the angular 
distributions of the beta-particles emitted in transitions between magnetic 
sub-levels of the initial and final nuclear states. These angular distributions. 
are required in the theory of beta~-gamma correlations (Falkoff and Uhlenbeck 
1950) and of emissions from oriented nuclei (Spiers 1949), and hitherto have 
only been calculated for first and second forbidden transitions in the Z=0 
approximation (Falkoff and Uhlenbeck 1950). The formulae obtained here 
apply to an arbitrarily forbidden transition and take account of the coulomb field 
of the nucleus. 

§2. ANGULAR DISTRIBUTIONS 

Consider a transition from an initial nuclear state with angular momentum 
quantum numbers J, ; to a final nuclear state J;M,;, resulting in the emission of 
an electron—neutrino system specified by the total angular momentum quantum 
numbers J, M= M,— M;. 

The state of the electron—neutrino system can then be written in the 
notation of I as 


-) M;M. = ay 
XMjM¢— x eM yiiaya x Cia! Vimal VO P)bjrma(t' oP’). eh eijaaye. 6 (1) 
Jjaj'a’ mm 


The ‘differential angular distribution’ (Falkoff and Uhlenbeck 1950) of the 
electrons with energy between FE and F'+ dE is then given by 


E 2 ‘ 
Ty, (9; E) dQ aE = are XmjMylr—r R? dQ dE, eovvee (2) 


where S represents an integral over the neutrino coordinates and a sum over the 
variables ao’8f’. R is a large but constant distance. 

The normalization of (2) is such that on averaging over initial and summing 
over final nuclear orientations 


1 
25,41 eel Tug >) asd EE}. Seca (Et 
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Substitution of (1) into (2), using (13) and (21) of I and the normalization 
properties of ¢, gives 


2 2 
Laguq(0sB)= SED] E MrasCaibta RHE (bs & )* CH Vin ROB 08) 
ja'm Bo j 
with 2 = Mtn, M=M;—M- UC (4) 
The summation over J, L, S in the square modulus may lead to many cross 
terms involving unknown nuclear matrix elements M,,, if the nuclear transition 
is such that a number of states J, L, S have the same (lowest) degree of 
forbiddenness (see table 1 of I). In order to obtain more definite results we 
consider here only transitions in which all but one matrix element may in fact 
be neglected. In this case 


Tyas (Os 2) = | Ceara Pha Os Ee ee (5) 
cee eaten Crd i [Mask 2 SE RES(. a P)*CHnntbinal ROS 0B) 
an (6) 


F37S(0; E) could be obtained from (6) by a similar procedure to that used for 
the energy spectrum (I, § 6). 

However, for beta-gamma correlation theory (Falkoff and Uhlenbeck 1950) 
and for the theory of emission from oriented nuclei (Spiers 1949) it is sufficient 
(and much easier) to calculate F¥75(0; E) only, since the properties of angular 
distributions under spatial rotation are such that 


FYS(; E) =| DPw(0)PFESO;B), sees (7) 
where the matrix elements D{),,(0) are the irreducible representation of order J 
of the group of real spatial rotations. ‘This reduces, for the special case M =0, to 
fe 
DIT 


Now for large R the electron wave function has the asymptotic form (Rose 
1937, using the notation of I) 


FIS(9; B) = Lar | YM (0) PFSO;E). eee. (8) 


nal ROG; 08)->(“EE)” Chizeay pesos) APUPR Eat fol (a) 
where the phase factors €,, and ¢, are dependent on the variables indicated and 
on the screened coulomb field through which the electron moves. 

In the absence of screening 


‘va — arg T(y+iaeZE/p)+n—yr/2 and C3 =(«ZE/p) In (2pR)—7(1+ B)/2 
where y=[(j+3)?—02Z?]"? and exp (2in)=[a(j+4)+iaZ/p]/[y+iaZE/p]. « is 
the fine structure constant. 


‘To take screening into account these phase factors would have to be computed 
numerically. 


Substituting from (9) into (6) and writing @=0 we obtain, after some algebra, 
2 2 L 
FYS0;8)= “EEEE | Myzalt & {Cady +32 yg @ptr gern 


(BES MAB + esas) PE 5 y2scamygopergra>l 10) 
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which is the required result. The positive or negative sign is to be used if the 
matrix element involves one of the Dirac operators 1, 6, a, y, or B, Bo, Ba, By; 
respectively. The coefficients «7"(M), ... . «78(M) are given in the Appendix.t+ 
The 4(n=1, , 5) are functions ay D, p related to the radial solutions of 
Dirac’s equation fore a stenteal field. ‘They are so defined as to reduce to unity 
in the limit Z=0. 
In the notation of I: 


[(2v+3)! 1] [(2v+3) IP 


go = per f= QprXr+2) ((FyiaP+ [Poe i Ajpaieee 
4 2v+3 aZ re 
=1+ — ae OE + O(a VAs ) 


go = —FQr+ 3) Qv+ I! py iH (2v +3)! v4)! 
‘aa Coa mas aa 7 


2v+3 «ZH 
v+1 2pp? 


ah ey Fa po eet 


[(2v + 1)! !]? 2v+1)!!)? 
¢® = ee P= SS 


p= —[(2v+3)!!? 
a Qpr+2) 


+ O(a2Z2) 
(UF ™ P+ [Fe P)p-? = 1 + O(02Z?) 


1-1 Fl 1 ly) .—2 1 
(hes cos 7,-1.7, FF. 1 Cos 1,)p "FP 


My 2v+3 aZ aan 

=1+ ary Wa CRA 

1E(2v + 3)!!(2v+1)!! . in Rote 

$= ee ee Che! Cosa, al, lt 2 COS Ty, oe ar 
=1+40(a2Z3) - 


where 7,=£,,, 1-431 and mj=7—€,,, 1 +,431- 

Both these factors are zero for Z=0. 

As shown above, the first three ¢’s are related to the f’s used in the expression 
for the energy spectrum in paper I, and do not involve the phase factors €. 


§3. RELATION TO THE ENERGY SPECTRUM 
Using (3) and (7) gives 


1 


P!1S(B) = (a Pe % : Ty,u,(0; 2) 4Q.= ait FYIS(O; BE). .ce0e (12) 


It is easy to show that the sums over M of the coefficients aJS( M), etc. are 
as follows: 


25 +1 | ape 
JLS/ yy) — oS yy JES’ IY) = JIS 
Sega emanate (13) 
2J +1 vet 
J LS = ILS JLS M)=> JLS = 
Pe moar ise SS Ber Onl 


where the «/”* etc. are defined in I. 


+ Each term of the summation over v corresponds to a total angular momentum j= v+#4 of the 
electron. 
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Hence from (10) 
yard 8 On = . 2»Eg?| Myzs|? = age pDpXrt+ Dg »—1) 
a SPP MESES [((2v+3)! 2°” 


27 +1) 92 L—»)—1 
JLS (2) P (y+ \g?k ”) 


» Aine: (3) p2v me 
= Qr43)l (v4 yt! * E “TaD |e 
tars (14) 

From the relation of 4% to f (n=1, 2, 3) set out in (11) we see that this agrees 
with eqn. (23) of I and is therefore a check on the calculations of this paper. 

The expression (14) for the energy spectrum has the advantage over the one 
given previously (I, eqn. (23)) that the ¢ functions have a simple form for small «Z 
and that the numerical coefficients are sufficiently simple to be written explicitly 
for the three cases of importance 


PLACE) = * epkg? U,| Mi 2 
: aes (1) 2(» +1) ,2(L—»—1) 
re: aon ee omperrng 
(Lv) t Qo i Pa 


“(202 be iiCuea E 
(v+1)24 (3) p20 tae 
+ OEE E Rye vale’ o 
4 i 2h P 
L+1, £,1 Ses eee 2 2 (3) 42" ,2(L— v) 
P (Z) = — g*pkq U;,| Mzs3,1] ® forcaepioen q 


IG2o | 
1B 


4 
PUA(E) = ~ gpg? —— U,| Mak 


L QL-1 
AO yay) 

: a; {Ore $y P q 

L(4) Ss OY ald s tabaiicta 

“Y (2L—20 +1 (2y4+3yl 7” E 
sesh Foes (3) p2v .2(L— v) 

 @L—2+Dl(ar+ ivr 1yL+ 2 2 

with the same sign convention as for (10). 


§4..PROPERTIES OF FU#S(9 ; E) 

Inspection of the coefficients «/”5(M) etc. (Appendix) shows that they all 
possess the property «/”*(M)=«/"S(— M) etc., so that F{75(0 ; E) = F/48(0; B), 
and hence F¥7'°(0; £) = F/44(0; E); this implies that FY?S(6; E) is a polynomial 
in even powers of cos 0. 

Further, the way in which @ is involved in (6) and (7) serves to limit the 
maximum power 2n of cos @ thus: is the largest integer such that n < ye a 

As already mentioned, the phase factors 7, and z,’ occurring in dS and 6 
are dependent on the coulomb field through which the beta-particles move, and 
for a screened field would have to be computed numerically. It should be noted, 
however, that in the important case of F¥+1(6;E) the functions 6 and 6? 
are not involved since =8 =8=e=0, Theexpression for the angular distribution 
is then independent of any possible screening effects except in so far as these may 
affect the form of the Dirac radial function at r=p (Reitz 1950). Fit6; FE) 
relates to nuclear transitions in which the parity change is such that |J,;—J-| 
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is greater by one than the degree of forbiddenness (see table 1 of paper I). These 
transitions have the further property that only one state L+1, L, 1 of the 
electron—neutrino system occurs with non-negligible probability; the angular 
distribution and energy spectrum are then uniquely defined, while for all other 
transitions several matrix elements may be non-negligible with unknown ratios 
and cross-terms result. 

§5. PARTICULAR CASES 


The following are tabulated as they constitute a generalization for ZO of 
the angular distribution functions given by Falkoff and Uhlenbeck (1950) for 
use in the study of beta-gamma correlations. 

We use b= +1 for the operators 1, a, 6, y; and b= —1 for B, Ba, Bo, By;. 


(i) L=1 (First or Second Forbidden Transitions) 
Fq(8; E) =3d$ + (49/9 — 2489/9 + 6/9) p? + b(245?/9 — 24(?/9)p2q/E 
+ cos? O[(d{?/3 + 245?/3)p? + 2bp*qd{”/3E] 
Fo(8; E) =$q?/6 + (49/9 + dP/36 + d0?/9)p? — b(269?/9 + 69/9) p?q/E 
+ cos? O[(p{?/12 — $5?/3)p? + b6$"p7q/3E] 
Fy"(8; E) =$§¢?/3 + b2'p?/6 + 3h) p? cos? 0. 


(ii) L=2 (Second or Third Forbidden Transitions) 
FOG; BE) = 2G g*/15 + (46§/27 + 249?/9 — 466/27) p2q? 
+ (246/75 +349)/50 + 246/25) p4 
+ b(8h§?/45 — 46 (?/45)p7q3/E + b(4g?)/45 + 26{?/15)p4q/E 
+ cos? O[(2p{?/9 + 46{)/9) pg? + (26{?/25 — 369/25 — 24g{")/25)p 
+4bpi)p?g?/15E + b(4g?/15 — 84{”/5)p°q/E] 
+ cos* 6[(3¢5?/10 + 66{/5)p* + 2b6)?p4q/E] 
F578; LE) =h§q4/45 + (P§?/27 + p?/54)p2g? + (PP?/150 + gY/150 — 4/75) p* 
; — 2bpip?g? /45E + b(PP/45 — P/45)p*g/E 
+ cos O[ p9'p?q?/18 + (f{?/50 — 6§?/75 + 464?/25)p* 
— b(pP/15 +4h/15)p4q/E] 
+ cos* 6[(f)?/30 — p{?/5)p* + bgp'q/3 E] 
F321(4; L) = 4g'/30 + $Op%42/15 + 664/40 
+ cos? O[26@p2q2/15 — 694/20] + cost O[69p4/8]. 


§6. CONCLUSIONS 


By using eigenstates of the total angular momentum of the system 
electron—neutrino, it has been found possible to derive a general formula (valid for 
arbitrary Z) for the angular distribution of the charged beta-particles in a 
transition between magnetic sub-levels of nuclear states of arbitrary spin and 
parity. We have considered explicitly only cases where all but one nuclear 
matrix element may be neglected, but the method is of course applicable to the 
general case if required. 

A subsequent paper will deal in a similar fashion with electron—neutrino 
correlations for arbitrary degree of forbiddenness. 
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«> GMP Bae N DEX 
We give below the values of the coefficients «7“5(M) . oe <J¥S(M) used in 
the text. 
(i) Scalar Operators (S=0, J =L) 
U,2"-1(L + M)!(L— M)! 
Loy yy) — LE 
ad Wheat OE) All 
1 
« {ET HSPSD IES on Ee wooaeseai} 
2v+3 
LL ) — LLO 
ja (M)= 97-9, 41% (M) 
emaarts U,2#-1  (L+M)!(L—M)! 
7 M)= (5) te (M+ 9, {Ou AOL vale 
(L— M—2Myr+1)? (L+M+2Mv+1) \ 
oon M—y+i)\. @2M-)NL eM 
pum) Ue __ (+ M) (LM) 
DTT (ave se Re EDO ira) 
L-3M-2Mr+1 L+3M+2Mv+li1 } 
WEEMEISS Loco * Ea 
2v+3 
LLO LL0 
ey (M)= 57 — ay 71° (M). 


(it) Vector Operators (S=1, J=L+1 or L) 
ob +11/ 17) = pitt) = OUT LE ME) = eb ti 1) =0 
a 24 (L4+1+M)!(L+1-M)! 
+111 
a ie Li (L+1)(2L+41)(v!)? 


1 1 
;. Pa tt eae eae = MAES 
ee Hee 


ahh M) = aL M), pLe(M)=— —— pt M) 
2v+3 U oe 1 (L+ M)!(L— M)!(2v—L)? 
LIA as LIA L ae 
Uae (8 lone 2L! [(2v+1)(2L—2v4+1)(L+1)v!} 
(L—M-—2Mv+1) Ri (L+ M+2Mv+1) 
(L+M—v)!(L—-M-—v+1)! (L—-M-v)!(L+ M—v+1)! 
symp) = — PAERS susan), eta) = FET casey yy, 
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Studies in Graphite and Related Compounds 
I: Electronic Band Structure in Graphite 
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ABSTRACT. The Bloch (tight-binding) approximation is applied to a study of graphite. 
‘The mean energy per atom and the internuclear distance are calculated in good agreement 
with experiment. Previous accounts are extended by including certain overlap integrals 
and by dealing with a set of parallel layer planes instead of only one plane. The shape 
of the band due to the z-electrons is calculated, and the experimental x-ray emission 
spectrum of graphite is interpreted. It is concluded that the role of c-electrons is greater 
than is often supposed. 


§i. INTRODUCTION 


HE fact that in graphite the individual planes are separated by about 
| 3-35 A, whereas the interatomic distance within any one plane is only 1-42, 
shows conclusively that we regard each hexagonal plane as a large ‘aromatic’ 
molecule, and the binding between the planes as essentially van der Waals in 
character. Current theories of aromatic molecules divide the bonds into two 
types—first there are the o-bonds, localized in the region between two adjacent 
nuclei, and described by electronic wave functions with approximately axial 
symmetry: second there are the z-bonds due to z-electrons. ‘These latter are 
the metallic electrons, since we represent them by orbitals extending over the 
complete plane. It is these latter, and not the former, which give graphite its 
characteristic properties. Discussions of the electronic structure of graphite, 
with the solitary exception of some early calculations by Hund and Mrowka 
(1937), have dealt only with the z-electrons. We may regard these as having 
energies lying within a band, the z-band, and as contributing to the observed 
interatomic distance and heat of sublimation. In two following papers we shall 
see that (a) these z-electrons contribute to the x-ray scattering in a characteristic 
manner, and (0) in the boron nitride crystal, which may be regarded as derived 
from graphite by alternately increasing and decreasing the nuclear charges by 
unity, these electrons exhibit certain differences from graphite. 

A preliminary account of this work has already been given (Coulson 1947). 
Since then a paper by Wallace (1947) has shown how the z-electron distribution 
is related to the electrical conductivity, though he did not attempt to draw the 
band distribution, and a further paper by Barriol and Metzger (1950) has 
reproduced part of the work reported by Coulson (1947). ‘This will enable us 
to abbreviate the first part of the present discussion though, as we shall show, 
there is one serious inaccuracy in this earlier work. In the later part of this paper 
we shall show the influence on the band shape of an overlap factor previously 
neglected, and the way in which the presence of other planes modifies the details 
of this shape. We conclude with a discussion of the relationship between our 
calculated energy—density curve and the observed x-ray emission bands. 


* Now at Mathematical Institute, Oxford. + Now at T'.R.E., Malvern, Worcs. 
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§2. TIGHT-BINDING APPROXIMATION WAVE FUNCTIONS 
FOR AN INFINITE SHEET 
Figure 1 shows a unit cell, which contains two atoms. Let us represent the 
normalized z-orbital of atom (j,k) by ¢,; and that of atom [j,k] by 4, ,. We 
may consider a plane array such as that of fig. 2 in which there are altogether 
2UV atoms. Later both U and V will tend to infinity. It can be seen from 
fig. 2 that the allowed values of and k are 


j= —n, Tele. Ot ee U=2u+1, 
k= 0; —U+td, cee LO ae oe Vatoat. 


If we introduce Born’s cyclic condition (or, alternatively, bend the large sheet 
of fig. 2 into a closed torus) we can use symmetry conditions to give us nearly 


AW yp : A w 
WwW oe v 
Fig. 1. Nomenclature for atoms in a Fig. 2. Nomenclature for the unit cells in 
unit cell of a graphite layer. a graphite layer. 


all that we require about the Bloch wave functions. Let us introduce the auxiliary 


functions 
a pee _{jm kn 
1 Night Sao ng hb wly Ree Sym {= As > 
; j=—uk=—-v U V 
2 +u +0 F im kn 
Ti ges SS boli degpexpem {e+ ts ieee (1) 
j=-uk=—v Cenk 
where —u<m<u, —v<n<v. These auxiliary functions are such that the only 
combining pairs are ®,, , and Y’,, , with the same m and n. This means that 
each molecular wave function (or molecular orbital, m.o.) is of the form 


CE ary ME | Piss, Reet pee we (2) 


mn 
where N,,,, 18 a normalizing constant and c,,,, is a complex quantity, of modulus. 
unity, whose value, for given m and n, has to be found from the secular equations. 
Let £,,,, be the energy associated with (‘),,,,, and Ey the energy of a separate: 
atomic orbital 4,;,; further let 8 be the resonance snkectead across a bond. This. 
means that1f H is the one-electron Hamiltonian, then 


= | bixgllyp;,ar = J diptlibaay pat =etce Ne ane (3) 


We require two other quantities. The first is the overlap integral between: 
adjacent pairs of atomic orbitals, defined by 


S= [bsbiedr=[Ssabsas, p OT = ClC. Pay Mans ences (4), 


as en 


mn 
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The second is closely related to the resonance integral. It is the quantity 
introduced first by Mulliken, Rieke and Brown (1941) and defined by 
y=B-E,S. Bye es 
Resonance and overlap integrals between other than nearest neighbours are 
temporarily neglected. 
We are now in a position to write down the secular equations, which are 
quadratics. A quite straightforward analysis gives 


ze g 
Finn =Eot+ lease? @eccece (6) 


where 


V. 


It may be seen that g? ranges continuously from 0 to 9. The details of the 
corresponding wave function are given by 


2nim Qin 
Cmn = 41+ exp—=— +ex Ft | THY TLS Pe recs et 8 
{ = aay ty UES) (8) 


td ov i(9) 0 and) Nj 2=2UV(1+2S). 0.0. (10) ° 


If we assume, as is often done, that the overlap integral S is zero, these results 
simplify to 


27m 2mn M~ -n 
m, n)|? =3 +2 cos — +2 cos — +2 c0s2 (=-=). aecteiets c 
[g(m, 7)] = z aC (7) 


so that c 


Both pop’ N62 2UV <8 te) (11) 


These last two equations are implicit in the work of Wallace, and of Barriol and 
Metzger. Since, however, the numerical value of S is known (Wheland 1941) 
to be approximately }, it is undesirable to omit it completely as these authors 
did. Furthermore, as Wheland (1941) and Mulliken, Rieke and Brown (1941) 
have shown for aromatic and conjugated molecules, agreement between theory 
and experiment for the heats of sublimation and energies of spectral transitions 
is impossible without its inclusion. We shall find that in our case the shape 
of the complete band of z-electron energies is greatly affected by the presence 
of S. In numerical work we shall put 


SLE Pe Es Pe (12) 


This value is obtained by direct evaluation of the integral (4) when Slater-type 
atomic orbitals are used (Slater 1930), together with the observed value of the 
interatomic distance. The same is true for the inter-layer overlap (0-034, see 
p. 821) between atoms facing each other in adjacent planes. Explicit formulae 
for these integrals may be found in Coulson (1942) and numerical tables in 
Mulliken, Rieke, Orloff and Orloff (1949). 


§3. MEAN ENERGY AND BOND LENGTH 


Equations (6)-(11) enable us, without great difficulty, to estimate the mean 
energy per atom and the internuclear distance, or bond length. Since there is 
exactly one z-electron per carbon atom it is simple to see from (6) and (7) that 
the occupied levels are just those for which g(m,n) is positive. ‘These occupy 
orbitals which, in molecular language, would be called bonding. ‘The others, 
for which g(m, 7) is negative, are antibonding. In the ground state, at absolute 
zero of temperature, they will all be empty. Thus the mean energy per atom is 
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E= 2m. nLmni2UV, where the factor 2 is introduced on account of the two 
spin assignments, and the summation is over all m and m, choosing the positive g 
in eqn. (7). | 
It is simplest, as Coulson and Rushbrooke (1948) have shown for " graphitic j 
strips, to deal with the simpler case of S=0 first. Then, on letting U, V tend 


to oo, we have 
E-E,= || {3+2cosx+2cosy+2cos(x—y)}!?dxdy,... ...... (13) 
where the range of integration is —7<x, y<7. This reduces to 
= 2 
Bone *y EEA + eos dls = een ee (14) 
WT 0 - 


where k? =8cost/(1+2cost)?, and E(k) is the complete elliptic integral of the 
second kind defined by 


7t/ 2 
E(k) = { (1 —R® sin? OY? dO. 
0 


Tables of E(k) are available, but it does not appear possible to complete the ~ 
integration in (14) analytically. Numerical integration gives 


EAE =157665" — J 9, 0s eee (15) 


This is the same as the value given without proof by Wheland (1934) and Coulson 
(1947). It agrees to within 1°% with the value obtained by Barriol and Metzger 
(1950) by considering particular values of U and V. 

If we include the overlap integral S, then the analysis _ is a little more 
complicated. The result is 

BF 1085) 00 gage eee re (16) 
where y is defined by (5). 

It will be noticed that the values in (15) and (16) fit smoothly into the tables 
of binding energies in graphite crystallites of varying finite size as calculated by 
Bradburn, Coulson and Rushbrooke (1948) (see table I on page 340 and table IV 
on page 345). If we take conventional values for 8 and y, then as shown in this 
latter paper, the numerical values in (15) and (16) do fit extremely well with the 
observed heat of sublimation of graphite. 

Since all the bonds are equal and there are 3/2 bonds per atom, we can find 


the z-bond order (Coulson 1939) by dividing E—E, by 38. ‘The result is 
P=0:525.-Fy il Vie Bee eb peat (17) 


As Chirgwin and Coulson (1950) have shown, this value is unaffected by inclusion 
of the overlap integral S. Assuming 1-540, 1-340 and 1-204A respectively for 
the bond lengths of C—C, C=C and C=C the corresponding bond length for 
graphite is 142A. ‘The complete agreement of this value (Bacon 1950) and that 
for the heat of sublimation (Jessup 1938) with the experimental values, provides 
strong support for the general picture that we have been describing. It also 
suggests that the inter-layer energy is quite small, perhaps only of the order of 
5 kcal/g atom. This would be entirely reasonable, even though this is not an 
adequate means of calculating it, on the basis of likely van der Waals forces. 
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§4. DENSITY OF STATES IN THE 7-BAND 
According to (11), with S=0, the energy levels lie in a band stretching from 
Ey +2maxB to Ey— Smax Bo It is dasily verified that the maximum value of g(m, n) 
. when m=n=0, giving 2yax=3. Thus the band stretches from E,+3f to 
E,—3B. So long as we put S=0, it will be a symmetrical band, just as with 
the strips discussed by Coulson and Rushbrooke (1948). ‘The density of states 
within the band is easily found if we notice that as U, V-> 0, the variables 


m/U, n/V become uniformly distributed between +3. It is convenient to 
introduce the function P(g) where P(g) is the proportion of occupied states 
whose g-values lie between g=3 (the bottom of the band) and our chosen value 
of g. Evidently P(3)=0, P(0)=1. It can also be shown that P(1)=2. Other 
values of P can be found, if desired, by numerical integration in the plane defined 


by the continuous variables m/U and n/V. The following asymptotic expressions 


are available: v3 
Near bottom of band P(g) = -— Ser —g)+...=0-276(3 —g)+.... 
Near top occupied region roe —g?/4/3r+...=1—-0-18497 +.... ...... (18) 


0 05 1-0 15 20 25 30 
IJ 
Fig. 3. The shape of the curve P(g), where P(g) is the probability that the function g(x, y) lies 
between g and 3. 


Since P(g) is the proportion of occupied states for which (H— E))/6 >g, it follows 
that —dP/dg is the actual density of energies per carbon atom per unit range of 
(E—E,)/B. Thus at the bottom of the band the density of states is 0-276, and 
near the top of the occupied region, where E— E, and g are small, the density is 
0-368g. The complete density curve N(£) is shown in fig. 4(a). Clearly 
N(£) = —(1/B8)dP/dg. ‘This curve includes both the occupied and unoccupied 
parts (at zero temperature). It is possible, however, to obtain an integral 
representation for P(g) and, which is more important, dP/dg, as follows. If 
we restrict ourselves, for the moment, to the two-dimensional layer with S=0, 
then we know from (11) that E,,,,=£,+gf, where g is given by (7). If we let 
x( =2nm|U) and y(=2nn/V) be current variables appropriate to the infinite 
layer for which U, V— o as in (13), then E—E, =g(x, y)8, where 
Lela) *=3 2 cosx-h2.cosy-+2.cos(¥—4)! Ale esas (19) 


The range of allowed values may be taken to be |x| <z, |y| <7. E-E, is 
a function of x, y represented by the contours shown in fig. 5. Then P(g) is 
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equal to 1/47? times the area between the origin and the contour g. Let us 
first suppose that g>1 so that the contour is single, and of approximately 
elliptical shape. ‘Then on rotating the axes through 45° in an anti-clockwise 
direction we have 


g?={1 +4 cos (¥/2"”) cos (X/242) +4 cos®(X/242)}. ws... (20) 


max 


2 
Then P(g)= a «4 | Y dX, where X ax = 21? cos {(g —1)/2}, and X and 


Y are related by (20). A certain amount of algebra enables this to be put 


in the form ! : 
eat —1-4t dt 
Me) | (o-P2 ea (* at Cae 
From this, by differentiation under the integral sign, it follows that 
dP 4p /t dt 
ae ele [(1 — #){(22 + 1)? —g%} fg? — (22— 1)" 
The substitution £2? = sin? 6+ {4(g—1)}* cos? 6 allows this to be expressed as 


“ dP ee en in do a mI 
dg a \a o [l=ecos* 6}? +B cos? bf eh) 


where «=(3—g)(1+g)/4, B=(3-—g)(1+g)/(g—1)(g+3). Both « and f are 
positive. 
An analogous formula can be obtained for the region 1>g>0. It is 
Epes Lee Se 
dg a Jo [hy cost OP 14 dcas 0p ee ee 


where y =49/(1+g)?, O<y <1; 5=4g/(1—g)(3+g), 5>0. From (21) and (22) it 
is possible, by simple numerical integration, to calculate N(£) for any selected 
value of g. ‘This technique is better than numerical differentiation of the P(g) 
curve. 

There are two improvements to be incorporated in this treatment. The 
first concerns the inclusion of S, the second involves the consideration of more 
than one layer plane. 


§5. INCLUSION OF OVERLAP INTEGRAL 


It is possible, as shown by Coulson and Rushbrooke (1948), to change from 
the case S=0 to the case S40 extremely easily. The new density function, 
which we may call N,(4), is related to the former one N,(£) by the relation 

NAB) = (TP BAN ge). et eh (23) 
N,(£) is the number of electrons that can be accommodated per carbon atom 
per unit range of (E—£))/y. We show, in fig. 4(b), the graph of N,(E). It 
is at once clear that the previous symmetry around the centre EF = Ej is no longer 
present. Just as in the earlier work of Coulson and Rushbrooke, the occupied 
part of the band is much compressed as compared with the unoccupied. But 
the characteristic drop to zero at the top of the occupied part, which formed the 
basis of Wallace’s discussion of conductivity, is unaffected. Indeed, a 


combination of (18) and (23) shows that not even the slope of the N(E) curve 
is changed at H= E). 
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§6. INTRODUCTION OF THE THIRD DIMENSION 


The second improvement in our calculated band shape is the result of 
introducing the whole set of layer planes instead of concentrating solely on 
one of them. Let us first consider a somewhat idealized structure in which each 
carbon atom has other carbon atoms directly above and below it. (In the real 
graphite this is true for only half the atoms.) Then in the conventional Bloch 
treatment we modify the wave functions (1) by introducing a further exponential 
factor to deal with motion across the planes. The energy splits into two parts 
—one part is the two-dimensional energy given in (6) and (7): the other part 
is an energy due to the interplanar momentum. This latter part is, formally, 
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Region 


=———_F = “3B, +2, 


- ess 
Ely 


— fF 


Fig. 4. (a) Graph of N(£) for a single layer with conventional neglect of overlap integrals. 
(6) Graph of N(£) for a single layer taking account of overlap integrals. 
(c) Graph of N(4) for three-dimensional crystal, with neglect of overlap integrals. 
(d) Graph of N(Z) for three-dimensional crystal, taking account of overlap integrals. 


just the same as for a linear chain of atoms, except that, on account of the fact 
that the overlapping lobes of z-orbitals on adjacent planes have opposite algebraic 
signs, the sign of the resonance integral is altered. Thus for any fixed m/U, n/V 
there is nowaband of energies. This band is quite thin, since its width depends on 
the numerical value of the resonance integral. This is not known at all precisely, 
but if, following suggestions of Mulliken, we suppose that the resonance integral 
is proportional to the corresponding overlap integral, we get for the ratio of the 
inter-layer resonance integral to the intra-layer resonance integral a value 
0-034 : 0-25. If we consider, for simplicity, the case where S=0, and denote 
the resonance integrals by 8, and f, respectively, then the total width of the 
two-dimensional part is 6f;,, and of the part due to the third dimension it is 48). 
PROC. PHYS. SOC. LXV, To—A 54 


822 C. A. Coulson and R. Taylor 


This width is less than one-tenth of the other. This means that we can pass 
from the diagrams in figs. 4(a) and 4(5) to corresponding three-dimensional — 
distributions by replacing each separate level by a thin band, centred at the 
same energy value and spreading a small, but constant, distance on either side. 
More accurate calculations than this can be made formally, but they are not 
justified until we know more about the values of 8, and are prepared to include 
B, for second-nearest neighbours as well as nearest ones. Such more detailed 
calculations are, of course, necessary in discussing the shape of the band near 
the point E=£,, where the spreading of each two-dimensional level into a band — 
of constant width is no longer a quite fair description. ‘This would apply to 
conduction phenomena and electronic specific heat calculations, for which the 
shape of the M(Z) curve near E=E, is quite crucial. Since we are not here 
interested in such matters, we shall not enter upon such a discussion, and shall 
be content with the simpler description over the whole range of energy. With 
this proviso the general result of passing from two to three dimensions is easy _ 
to see. Instead of fig. 4(a) we have fig. 4(c); instead of 4(b) we have 4(d). 
The mean energy is scarcely altered by this, nor is the bond length changed. 
We have shown in some unpublished work that the effect of taking real graphite 
instead of our idealized form is substantially the same as above, except that 
each separate two-dimensional level is broadened rather less than our calculations 
above suggest. ‘The infinity given by the two-dimensional model is removed. 
We enlarge on this later. The final density curve will differ very little from that 
shown in fig. 4(d). Wallace (1947) has discussed the small region of this curve 
surrounding the energy H=£,. His curve in this region is essentially the same 
as ours. 


7 


§7. COMMENTS ON THE N(E) CURVE 


There are three comments to be made on the shape of the N(£) curves in — 
fig. 4. The first is that figs. 4(@) and 4(b) show an infinity at EH-E,=f, or, 
equivalently, H—E,=0-8y. This follows from the fact that the integrals in (21) 
and (22) both diverge when g=1, i.e. B=85= 00 in (21) and (22). Now these 
energies are precisely those for an electron in an isolated double bond, such as 
occurs in ethylene. So it appears as if the most common distribution of bonds 
was one in which the electrons linked themselves in pairs to form independent 
C=C units. In this way we relate our work to that of Pauling, Brockway and 
Beach (1935) who determined the bond order of graphite by considering resonance 
of Heitler-London type among a large number of canonical structures in which 
the bonds were all either single or double. This resonance, it would appear, 
has the effect of spreading out the energies into a band, but the maximum density 
is still located at its value for any one of the unexcited canonical structures. 

Both of the previously published M(Z) curves are wrong because they failed 
to take account of this infinity at H—-E,=+ 8. It is worth while following our 
argument a little further, with a view to seeing whether additional refinements 
would remove the infinity; and if so, which ones would do so. We shall merely 
report the results of our analysis here, without proof. First, on account of (23) 
the introduction of S is irrelevant in this respect. But we might expect that some 
difference would result if, in addition to 8 between nearest neighbours (now 
called 8,) we introduced £, between next-nearest neighbours. Instead of (19) we 
now find that E—E,= +g6,—(g?—3)B,. Thus the shapes of the contours in 
the («, y) plane are entirely unchanged, though their labelling is slightly altered. 
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‘The infinity in N(£) still remains. However, if we introduce third-nearest 
neighbours, with corresponding 8 =f, we discover that 


(EZ — Eo) + (g? — 3) Bo}? = B,20(«, y) + Bg2e?(2«, 2y) 
+ Bi Bsl(g? — 3) +4 cos («+y){1 + cos x+cos y}]. 


It can be shown that the infinity still remains, though it occurs at a slightly 
different value of E than before. We have shown that an infinite value of N(E) 
will always occur, however many neighbours we consider, provided that we 
consider only a two-dimensional layer. However, as soon as we go to a three- 
dimensional structure, the infinity disappears. This is true for several other 
types of crystal structure besides the hexagonal one, though it does not hold 
for all structures, as the particular example of body-centred-cubic shows (Mott 
and Jones 1936, Chapter III). 

Our second comment on the N(£) curves in fig. 4 is that they are obtained, 
as stated in the beginning, on the assumption that o- and z-electrons can be 
treated separately. Now as Altmann (1952) and March (1952) and others 
(Coulson, Higgs and March 1951) have recently shown, the z-electrons merge 
into the distribution of the o-electrons. It seems possible that the lower 
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(occupied) part of the band is not greatly affected by this, but that the upper 
(unoccupied) part may be more changed. Our third comment is to the effect 
that since the method which we have used considers only atomic 2p7 orbitals, 
there will be mixing of the states which we have found with the band of states 
arising from 3p7 atomic orbitals. ‘The matter has been referred to, for molecules, 
by Coulson, Craig and Jacobs (1951). Again the conclusion seems to be that the 
lower half of the band is unlikely to suffer much alteration, but that the upper 
half may be more affected. 


§8. COMPARISON OF BAND SHAPE WITH EXPERIMENT 


It is possible to compare our predicted band shape in graphite with the 
experimental one, as found by emission experiments in which an electron in the 
occupied part of the band falls into an empty K orbital. If we may suppose, 
for lack of any better knowledge, that the transition probability from any part 
of the occupied band is the same, then the observed intensity distribution is a 
direct measure of the band density N(F). ‘The best measurements for graphite 
are those of Chalklin (1948) whose N(£) curve, on this assumption, is shown 
in fig. 6. At first sight this appears to agree excellently with the theoretical shape 

54-2 
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of fig. 4(d). But the agreement is partly deceptive, for the total width of this 
observed band is about 15 volts. ‘The theoretical width should be 1-7y=5 ev, 


‘since the resonance integral y has a value of about 3ev. Our confusion arises 


from having overlooked the o-band which must be superimposed on the 7-band. 
Now o-electrons have lower energy (greater binding energy) than z-electrons. 
So the region on the left of fig. 6 should relate to 7-electrons, that on the right 
to o-electrons. As there are three times as many o-electrons as there are 
z-electrons we shall not expect the z-electrons to show up clearly resolved from 
the o-electrons. The situation is very similar to that found for the charge 
distribution in benzene by March (1952). It may therefore be significant that 
the most noticeable difference between the N(E£) diagrams obtained by Chalklin 
for diamond (all c-bonds) and graphite, is that in the latter there is a substantial 
shoulder at about 281 ev which does not appear in the former. Now this is just 
where the z-band would be expected, as the following figures show. 

The ionization potential (work function) of graphite is about 4-3ev. This — 
corresponds to the removal of the top electron, and shows that for graphite 
E,=43ev. The calculated width of the occupied z-region is about 5ev, as 
we have seen, so that the 7-band would be expected to stretch from 4-3 to about | 
10 ev below ionization level, with an average at about 7ev. Now the ionization 
potential of a o-electron in a large organic molecule is known to be in the region 
of 12ev. So the centre of the o-band, which is a filled band, should lie about 
12ev below ionization. ‘Thus the centres of the occupied o- and z-regions 
should differ by about 12—7=5ev. If we take the shoulder on the Chalklin 
curve at 281 ev to be the z-band, then the rest is the o-band, and the two peaks 
do differ by about 5ev. Further, if we attempt to draw the separate bands, as 
is done somewhat schematically by the dotted lines in fig. 6, we do find that, in 
very rough terms, the area of the o-band is three or four times as large as that of the 
a-band and the width of the z-band is about right. It seems quite likely that 
our interpretation is essentially correct, though the details must not be pressed 
too closely. In the first place the transition probabilities for o- and z-electrons 
are unlikely to be equal, and in the second place, as Altmann (1952) has shown 
for benzene and ethylene, there is no rigorous division into o- and z-bands. © 
Mixing of the two component bands would ‘smudge’ the separation attempted . 
in fig. 6. But it should still be true that the left-hand (small wavelength) end 
of the emission band should. be predominantly polarized with the electric vector 
perpendicular to the layer planes and the right-hand (long wavelength) end 
should be polarized in the layer planes. 

Consideration of figs. 4 and 6 suggests that, although the main conclusions 
are probably correct, it would be very desirable to discuss both the o-band and 
the resulting hybrid more carefully. 
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ABSTRAC1. ‘The space wave functions for the z-electrons in a graphite layer, as found 
in Part I, are transformed into momentum coordinates, and the radial distribution function, 
including o-electrons and z-electrons, is calculated. Using the theory of Jauncey and 
DuMond the shape of the Compton profile in x-ray scattering is predicted. Both this, 
and the momentum distribution, are in excellent agreement with experimentally determined 
values. The introduction of a plausible scale factor in the atomic wave functions would. 
make the agreement almost perfect. 


$1. INTRODUCTION 


HE present authors have shown, in a series of papers, how the momentum 

| distribution of electrons in various atoms and molecules may be calculated 
(a complete set of references is given by Duncanson 1943). In view of 

the difficulty of handling the wave equation in momentum coordinates, we 
found it more convenient to adopt the Dirac transformation theory to derive 
the momentum wave functions from the corresponding space functions. Even 
this procedure is exceedingly complex, if any but the simplest approximate 
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wave functions are used, as the example of Hj showed. Now the momentum 
distribution cannot be determined directly by experiment, but it can be 
deduced from the detailed shape of the Compton line in x-ray scattering, and 
from the energy distribution of inelastically scattered electrons. In order 
that the experiments may be directly comparable with the theoretical work, 
the atoms or molecules must be effectively isolated from each other, as in a gas, ~ 
with the result that the intensity of the scattered x-ray or electron beam is weak. 
In addition, in the x-ray measurements, an extremely high dispersion is necessary. 
It is for this reason that only a limited amount of experimental work has been 
done which is at all suitable for comparison with theory. 

In view of this situation it seems desirable to consider the corresponding 
problem in relation to the solid state. The theoretical aspect is rather more 
forbidding, since at least some of the electrons will now be described in terms 
of non-atomic wave functions (or possibly combinations of more than one atomic 
orbital on different centres). We have thought it wise to discuss the case of 
graphite first, partly because some careful work by Kirkpatrick and DuMond 
(1938) has yielded a more accurate shape for the Compton line than is available 
for other solids, and partly also because the layer structure of graphite is such 
that it is a good approximation to treat it as a two-dimensional problem by 
neglecting the influence of one plane on its neighbours. The analysis in Part I 
(Coulson and Taylor 1952, to be referred to as I) justifies this procedure. The 
effect of adjacent layers can be calculated in a quite straightforward manner, 
but it does not seem to us worth the extra labour involved, until our knowledge 
of the distribution of the bonding electrons in any one plane is more complete 
than it is at present. 

Our procedure in this paper is (i) to discuss the hypothetical case of an 
infinite straight chain of carbon atoms, (ii) to generalize this and compare it with 
the true hexagonal layer, (111) to deduce the shape of the Compton pronle and 
compare with experiment. Our symbols and general model are the same as 
in I. There are three types of electron orbit to be considered. First there are 
the carbon K-shell electrons, in undistorted atomic orbitals; next there are 
the o-electrons of the localized o-bonds, and finally there are the z-electrons 
whose orbitals extend over the whole crystal. ‘The first two types of electron 
have been dealt with in our previous papers, since they have just the same form 
as in ordinary atoms and molecules. The new features occur with the z-electrons, 
with which, therefore, most of our analysis will have to be concerned. As a 
guide in dealing with the two-dimensional problem, we have first studied the 
momentum distribution of the z-electrons in an infinite straight chain of carbon 
atoms, where each atom provides just one z-electron. The directions of all these 


7-electron orbits are supposed to be parallel, and perpendicular to the direction 
of the chain. 


§2, MOMENTUM FUNCTIONS FOR THE 7-ELECTRONS 
IN AN INFINITE CHAIN 
If the chain is long enough, end effects are negligible and we may apply the 
Born cyclic condition. In this particular problem the use of this simplifying 
assumption 1s not necessary, for it may be shown rigorously to make no effective 
difference to the momentum distribution. But it is useful and, as we saw in i 
at is necessary in the two-dimensional problem. 
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Let us suppose that the repeating unit consists of m similar atoms. Then the 
molecular orbitals (Bloch functions) for the z-electrons are 


y Perec 
See ee OF ig BS Pe 1 
k Nice ib; ; ( ) 
where CS Explain) | oy ea (2) 


N, is a normalizing factor and ys, is the normalized atomic orbital of the p, 
electron at the jth atom. The analytical form of %; is taken to be 


= JF RpeXPee Ohya en (3) 


and the exponent « has the ordinary Slater value «=1-625. If we neglect all 
overlap integrals except those between nearest neighbours, then 


NAL 2S COS AnRiN)e i wags wk (4) 
where S is the overlap integral defined by 
S= | Ue are (> abe pee Ae (5) 


A knowledge of the space wave functions (1) enables us without difficulty to 
deduce the corresponding momentum functions. 

Adopting the Dirac transformation theory, and remembering that atomic 
units are being used throughout, the momentum function, corresponding to 
7,18 ¥;, where 


X~=(20)-3? | exp{-ip.n}¥,dr eae (6) 


If we substitute for ‘Y, from (1) and choose appropriate origins for each of the 
separate resulting integrals, this may be written 


vfs he WPA 
Xe = pee expii(— = pp) 0 NN aiciel ese tote (7) 


where p, is the component of the momentum p along the chain, p is the interatomic 
distance and A is the atomic momentum wave function 


A=(2n)-32 i exp{—i(p.r)}Wjdr. —  besase (8) 
The summation in (7) is straightforward, giving us 
<A . {27k 1—exp (—impp,) 
=H OP Pet ep ital ape ”) 


The corresponding momentum density function is 
paras sin” npp,|2 

XwXk —"N sin® (2ahk/n—pp,)|2- 

To find the mean momentum density function xy* we must sum terms such 

as (10) for all z-electrons and divide by the number x of electrons. ‘These occupy 

the 4n lowest energy orbitals, in pairs. These orbitals correspond to the range 

—nl4<k<+n/4. So 

n/4 sin” (npp,/2) 1 (11) 


He Pe eee eet ee 
XX aan Hes n* sin® (2ark/n — pp,)/2 1+25 cos (27k/n) 


(Strictly there are 47+ 1 terms in the sum (11) and we only need 4m. But when 
we let m tend to oo the extra term is irrelevant.) ‘This summation is difficult. 
But as n— oo it simplifies. For the first fraction after the summation sign in (11) 
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which we shall denote by f(pp,,) has the value zero except when pp, =2nk/n + any 
even multiple of 7. Thus 


2AA 
* eS eS 
us = Se, HeP x), 
where, after some rather tedious algebra, which we shall not reproduce, it may be 
shown that 


Kphs) =1 if p, lies between (41+ 1)n/2, / integral 
=0 otherwise. 


If we call these ranges 1 and 2 respectively, then 
xx* =2AA*/(1+2S cos pp «) range l, 
=() range 2." "ite (13) 


In this formula A is the atomic momentum function corresponding to any of 
the %,;. It may be found by use of (8) and straightforward integration. The 
result is 

23a he ps 
For comparison with experiment it is more useful to average over all angles. 
If we call the resulting function /(p), then J(p) is the mean radial momentum 
distribution. After some reduction this may be put in the form 


SE al Desi da ay, 
( eee | 1 +25 cos pp, apis) Ch Alot Vee (15) 


except that all values of p, in range 2 are excluded from the integration. We 
found it easiest to perform the integration graphically. From a knowledge 
of I(p) the Compton profile for scattered x-rays is soon found. We shall return 
to this later. 


AA* = 


§3. MOMENTUM FUNCTION FOR THE TWO-DIMENSIONAL 
GRAPHITE LAYER 
It is not difficult to generalize the preceding discussion to deal with the 
two-dimensional graphite layer. Again, let us concentrate first on the z-electrons. 
Both the nature of the unit cell, and the molecular-orbital wave functions, have 
been described in I. ‘Thus we write for a typical Bloch function 


(aes ia Nea ese ae Cana oe ers Nae (16) 


where ‘’,,,,, and ®,,,, are defined by I, eqn. (1); ¢,,, is a coefficient given by I (8) 
and N,,,, is a normalizing constant defined by I (10). The Dirac transformation 
theory leads to a momentum wave function (cf. eqns. (6) to (8)) which may be 
written 


Xmn = Nieeceuee a Cm Ges a eco (eheheye (17) 
eae En = (2n)-%8 [ Y,,, exp (—ip «r) dr, 


Grn aa =(27)° ial Gh exp ( is ip. r) dr. 


If we use the oa (1) of V’,,,, and the definition II (8) of A, we can write 


+v 


Be h ‘ 
le ea y xX exp 2a a + > exp (—sp [Gl s.} sea ee (19) 


j=—u k=—v 
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where [rl jx is the vector position of carbon atom numbered [jk]. A similar 
formula exists for £,,,, except that (r);, replaces [r],,. This means that (17) can 
be written in the form 


Keen Netnttkek Coon XD (IP. O)}ojas = * eee (20) 


where e is the vector joining [jk] to (jk). This means that p is the C-C bond 
length. The momentum density is : 


eS es =2N uy *{1 + | cosp . e+ cos (> — Pp e e) 


+ cos (= —p. e) |e} eh eae Ay ei RE ee (21) 


So far we have only needed to treat p as a vector. When its components 
are required, it is most convenient to express them relative to the reciprocal 
lattice. The relationship of this lattice to the unit cell is exhibited in fig. 1. 


by 


bj 


Fig. 1. Lattice vectors a and reciprocal vectors b for a graphite layer. 


Here the three fundamental space vectors are a,;, a,, a,, where a; and a, are each 

of length 1/3 directed at an angle of 60°, and a, is a vector, conveniently taken to. 

be of unit length, perpendicular to the layer plane. ‘The reciprocal vectors are 

b,, b,, b, where the magnitudes are 2/3p, 2/3p and 1 respectively: b, is parallel 

to a,; b, and b, are perpendicular to a, and a; and make an angle of 120°. 
Relative to these axes let us write 


Dee D tO pare te 8 8 foes s (22) 
Then [r]j,.=["loo+ja;+ka,, and p-[r];,=P-[Plo tjPj+kPy. 9 ---ee- (23) 


If we substitute from (23) into (19), it becomes possible to complete the summations: 
overjandk. Each separate summation is similar to that involved in passing from 
(7) to (9). The corresponding density function €,,,€,,,,* is similar to (10) except 
that there are two ‘diffraction’ factors instead of one. An analysis entirely 
similar to that which led from (10) to (13) shows that in the summation 

.. 1 


Ths a om mn eee, 
XX UV ahs x 
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which is needed to give the mean momentum density function, when U, V> o, 
the previous range 2 in eqn. (13) disappears, and for all values of the momentum p 
it turns out that we may write 

G( Fj, 
xx* = AA* ( J» i) 


1+ Sg(Pi, Px)” 


In this expression A(p) is the momentum function for a single atomic orbital, 


as in (8) and (14), g(P;, P;,) is our previous g-function written in terms of the 
momentum Si so that 
[eo(P;, P,)?=3 +2 cos P; ae oe Pls prea (25) 


and G(P,, P,,) is a function defined by 
G(P,;, P;,)=!+- _ loos 4(P,+ P,)+cos 4(2P;—P,)+cos#(2P,—P,)]. ...(26) 


If we recall ac P,, P,, P, are components of the momentum in the reciprocal 
space, we shall recognize that the units in which they are measured are 2/3p, 
2/3p and 1 respectively. This means that the magnitudes of the three components, 
in ordinary units, are p,;, p;, p, where 


p; =2P;/3p, Px =2P;/3p, p= P.. ante S (27) 


This analysis gives us the momentum distribution for the z-electrons. Before 
going further, however, there are two comments to make. In the first place 
the third component p, appears only in the term AA* of (24). This is a direct 
result of the single-layer character of our model. In the second place, despite 
appearances, the functions g and G do both possess trigonal symmetry. It may 
be shown algebraically that they are both unaffected by a rotation of the axes 
through 120°. 

Our next task is to average yy* overall angles. Just as in passing from (13) 

o (15) this will give us the mean radial distribution function J(p). In fact 
I(p) =integral of (24) over all angles, keeping the magnitude of p constant. If 
we imagine a series of concentric spheres drawn in p-space, then J(p) dp =integral 
of (24) over all space between the spheres p, p+ dp. Now the element of volume 
may be put in the form 34/3 dp; dp,(p/p,) dp. This shows us that 


G p v3 
eac a 


1(p)=2| | A(p)A*) Ci inks arid Orel (28) 


‘The initial factor of 2 is inserted to take account of both top and bottom halves 
of the sphere p=constant. If we now insert the value (14) for AA* this gives 


128+/3a? G 
I(p) = oS Aae (02 + p68 i De 145g dp; AP rs woke (29) 


where the integration is taken over the whole of the inside of a circle of radius p 
lying in the jk plane. ‘The relation between p, and the other components is 


pf =P" —pP— P24 pp. OS Fee (30) 


Also, when using (25) to obtain the value of g, the positive root is to be taken. 
The nature of the integrand in (29) shows that numerical, and not algebraic, 
methods are simplest for evaluating the integral. 


v= 
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_ This stage of our work is now complete. 'T'o obtain the total momentum 
distribution, we must add to (29) the contributions from the o-bonds and from 
the carbon Is electrons. Since these contributions are purely atomic and 
molecular, and have been discussed in our earlier papers (e.g. Duncanson and 
Coulson 1941) we do not reproduce the analysis here. 


§4. CALCULATIONS 


It is necessary to decide what value shall be taken for the exponent « in (29). 
As we stated in (3) this has been given the ordinary Slater value 1-625. The 
same exponent is used for the basic o-bonds. There is some doubt concerning 
the suitability of this value. We shall return to the matter a little later. The 
internuclear distance is known to be 1-421 A which is almost equal to the rounded-oft 
value 2-7 atomic units. We have therefore put p =2:7 A.U. 
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Fig. 2. Mean radial momentum distribution. Comparison of the radial momentum distribution 
of the p-electrons (or 7-electrons where relevant) in an isolated carbon atom, a carbon pair, 
an infinite chain of carbon atoms, a graphite layer with parameter «=1-625 and the same 
but with parameter «=1-9. 


Figure 2 shows the mean radial momentum distribution /(p) for several 
sets of z-electrons. This includes a graphite layer, an infinite carbon chain, 
a pair of carbon atoms as in ethylene and an isolated carbon atom. In all four 
cases « = 1-625 and (apart from the last one) p=2-7. It will be seen that as we 
move to systems of greater complexity the maximum in the /(p) curve moves to 
a smaller p, and the spread of the curve decreases. In passing from the infinite 
chain to the two-dimensional layer there is only a small change in the peak value, 
though the spread of the curve is reduced. We may reasonably suppose, on 
account of the large interlayer spacing, that the inclusion of other planes than 
the one we have considered would scarcely change the /(p) curve any more. 
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In fig. 3 we have drawn curves for the total momentum distribution, including 
the o-bonds and 1s electrons, weighted according to their number. This curve 
does not differ greatly from that appropriate to an isolated carbon atom. Just 
as in earlier papers we may use the curve of fig. 3 to calculate the shape of the 
Compton profile, using the type of argument advanced first by Jauncey. It has 
been shown (e.g. by DuMond 1933) that all incident wavelengths and all angles 
of scatter may be included in one formula, by use of a suitable ‘reduced’ variable. 
Thus the intensity J of the scattered x-rays is defined in terms of a variable q, 
where in DuMond’s notation g=//2A*, and / is the wavelength displacement 
from the centre of the Compton line, and A* is a quantity, with the dimensions 
of a wavelength, depending only on the primary wavelength and the angle of 
scattering. ‘The formula for J(q) is 


y(q)=h | Pap, Pale (31) 


where & is an arbitrary multiplying constant. If we require the total area of the 
J(q) curve to be unity, then k=}. But for convenience in comparison with 


experiment it is more usual to define k so that the peak value of J(q) is given by 
J(Q)= 1. 


—#— Experimental Compton Profile (Du Mond) 
—-+-- Theoretical Compton Profile 


06L ~-+-—Theoretical Curve for eb aa 06 
—2— Experimental Curve for Graphite 
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Fig. 3. Mean radial momentum distribution, Fig. 4. Shape of one half of the Compton 
a- and 7-electrons together. profile for graphite. - 


Figure 4 shows one half of the Compton profile for a graphite layer calculated 
in this way. We have also shown the experimental curves obtained by DuMond. 
The two curves are very similar, and differ chiefly in the fact that the half intensity 
is calculated to occur at g=0-96, whereas it is found experimentally at g=1-06. 
The - difference is only 10°, a surprisingly good agreement in view of the 
approximations which we have been obliged to make in our model. 

The greater half-width of the J(q) curve is in complete accord with what we 
have found before with molecules. ‘The matter has been discussed in some detail 
by ourselves (Coulson and Duncanson 1942) for methane CH,. It is true that 
there are no z-electrons in this molecule, but there is no reason to doubt that 
similar conclusions would be obtained with them. It was found that if the 
Slater values for the exponents in the atomic orbitals for the 2s and 2p electrons 
of carbon and for the hydrogen 1s electrons were each multiplied by a parameter f, 
then the half-width of the Compton profile was increased by about 10% if B 
was changed from 1-0 to 1-1. ‘The change in half-width was almost linear in B 
at least up tof =1-2. In the case of methane reasons were given for taking B = 1-1. 
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Chiefly these are an acknowledgment of the fact that an electron in a molecule is 
under attraction to two or more centres of force instead of only one, as with an 
atom. And, as specific calculations for several simple molecules have shown 
convincingly, this contracts the charge-cloud. This is precisely what happens 
if 8 is increased. Now it would be rather laborious to repeat all our numerical 
integrations for several values of «; we have, however, made as careful an estimate 
as we could of the shape of the J(p) curve when «=1:9. This is shown in fig. 2. 
Its relationship to the previous curve, which was appropriate to the case of 
ax=1-625, is of the required type to make agreement with experiment almost 
perfect. This may therefore be regarded as a further justification for the increased 
value of « as a result of bonding. 

DuMond also gives a curve for the mean radial momentum distribution [(p) 
in graphite. This he deduces from his experimental J(q) curve by means of the 
inverse relation to (31), viz. 


RI(p)=pldI(q)/@lony ees (32) 


We have reproduced his J(p) curve in fig. 3. Our theoretical curve here is 
normalized in such a way that JI(p) dp=1. Unfortunately it is difficult to do 
this for the experimental curve, since the experiments do not yield sufficiently 
accurate values of J(p) for large p. We have therefore adjusted its scale, for 
purposes of comparison, so that the maximum values of the experimental and 
theoretical curves are identical. Now, as (32) shows, the experimental curve is 
obtained by a differentiation process. It will therefore be less reliable than the 
undifferentiated J(q) curve in fig. 4, and there is little point in comparing the 
detailed shape of the two curves. In general terms the agreement is good. If 
we were prepared to change the theoretical curve by introducing the factor f 
previously discussed, it would shift it to higher values of p, and would lead to 
a quite remarkable measure of agreement with DuMond’s curve. 

There is one final comment that we wish to make. When dealing with atoms 
it is certainly the mean radial momentum distribution function which is important, 
since all directions of orientation of the atoms must occur. But with crystals, it 
may not be necessary to average over angles and we should need to work directly 
from (24) rather than from (29). Until experiments using a single crystal of 
graphite, instead of a powdered crystal, are available, it does not seem worth 
while making these calculations. 
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ABSTRACT. An analysis is made of the distribution of z-electrons in the crystal of 
boron nitride. The insulating character of this substance is shown by the fact that these 
electrons completely occupy one band, of width approximately 2-4 ev, separated from the 
next higher band by a gap of about 4-6 ev. The density of states in the band is calculated, 
accurately for a single layer plane, and approximately for the three-dimensional crystal. 
Until the energies of the o-electrons are better known, further progress in understanding 
this substance will be difficult. 


§1. INTRODUCTION 

ORON nitride bears the same relationship to graphite as borazole (B,N,H.) 
B does to benzene (C,H,). ‘Thus (see for example Bauer 1938) borazole, 

which has sometimes been called ‘inorganic benzene’ has a regular plane 
structure just like that of benzene, and boron nitride (Hassel 1926, Pease 1950) 
consists of parallel layer planes in each of which the boron and nitrogen atoms 
occupy alternate sites on a regular hexagonal lattice. The distance which separates 
the planes is almost the same (3-30 A) as in graphite (3-35 A) and the interatomic 
distances are 1-45A (boron nitride) and 1-424 (graphite). ‘There are, however, 
two essential differences between graphite and boron nitride. The first is that 
on account of the greater electronegativity of nitrogen as compared with boron 
there will be an uneven distribution of the electrons, giving a partly ionic character 
to the lattice. ‘The second is that, as Pease has recently shown, successive layers 
are stacked with their hexagons lying directly above and below each other, in 
such a way (fig. 1) that boron atoms lie between nitrogen atoms, and vice versa. 
This differs from graphite, where there is an additional parallel displacement 
between successive layers. ‘The boron nitride arrangement has the advantage 
that atoms with a net positive charge are completely surrounded by atoms with 
a net negative charge, and vice versa. It is doubtless the increased electrostatic 
attractive energy which makes this configuration the stable one. It is our 
intention, in this paper, to study the electron distribution in boron nitride, in 
much the same way that we studied graphite in Part I (Coulson and Taylor 1952 
to be referred to as I). We shall interpret the nature of the forces that cause 
the bonding, and estimate the band distribution (or, at least, that part of it which 
is due to the 7-electrons). A good deal of the analysis in I may be adapted without 
difficulty to apply to our present problem. 


§2. TYPES OF BONDING ELECTRON 
The hexagonal structure of a single layer suggests strongly that there is sp? 
hybridization at the boron and nitrogen atoms. This would lead to a o-bond 
distribution precisely as in graphite. But it would leave two z-electrons on each 
nitrogen atom and none on each boron atom. If we suppose that ionic structures 
* Now at T.R.E., Malvern, Worcs. + Now at Mathematical Institute, Oxford. 
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exist, in which one of the z-electrons on a nitrogen atom moves into one of the 
empty z-orbitals of an adjacent boron atom, then a double bond may be formed. 
In borazole, for example, we can imagine that the true wave function is a resonance 
hybrid of functions corresponding to structures such as (a), (b) and (c) below. 


B B Bo 
8 he a ALO 
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N N+ N+ 


(a) (b) (c) 


The stability of borazole may be attributed to the possibility of such resonance. 
Presumably the same situation will occur in boron nitride. If so, we may describe 
the bonding as due to a combination of (i) localized o-bonds joining adjacent 
BN atoms and (ii) z-orbitals which, in molecular-orbital language, would cover 
the whole set of boron and nitrogen atoms. Just as with graphite in I, there 
would be a o-band and a z-band. Arguments similar to those used before show 
that the distinction between o and z is somewhat blurred; and also, that the 
main characteristics of the band structure can be obtained by considering a 
single layer plane, since the effect of including adjacent planes will be to widen 
out each separate energy level into a thin band. This widening, just as in I, will 
exert only a small influence on the shape of the complete band. 


§3. THE 7-BAND SYSTEM 


Let us begin by considering the z-electron bands, and fix our attention on 
one plane. The unit cell is identical with that in I, fig. 1, except that the (jk) 
atoms may be taken to be all boron, and the [jk] atoms all nitrogen, as in fig. 2 


ae 
ee 
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Fig. 1. The BN lattice. Fig. 2. Nomenclature of atoms in the unit cell 
of a layer plane. 


of the present paper. It will be convenient again to consider an array of UV cells, 


containing 2UV 7z-electrons, just as in I, fig. 3. Later we shall let both U and V 
tend to infinity. 

Our analysis closely resembles that in I. ‘Thus we form auxiliary functions 
and Y’,,,, exactly as in I(1), except that ®,,,, contains only orbitals of boron 


® 


mr 
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atoms and ‘’,,,, contains only orbitals of nitrogen atoms. ‘The Bloch orbitals 
(¥),,, can be written, entirely as I(2), in the form 
CS \inm =) welt sigs aun asad: Piteurias (1) 

The coefficients c,,,, and the energy £,,,, of this orbital must be found from the 
secular equations. ‘These are quadratics. ‘They involve (i) 8, which is the 
resonance integral for a B—N bond, as in I (3), (ii) S which is the overlap integral, 
as in I (4), and (iii) xy, «pj which are the Coulomb terms for a nitrogen and boron 
atom respectively. These are defined in terms of the one-electron Hamiltonian H 
by the relations 


ae | ad fae Oe oe | Sidr (2) 


We neglect all overlap integrals and resonance integrals except between nearest 
neighbours. 
The secular equation whose roots are the two values of E,,,, is easily written 
down. It is 
ee D (B — ES){1 + exp (2mim]/U) + exp (2min/V)} 
=x),) Saas (3) 


(8 —ES){1 + exp (—2mim/U)+exp(—2min/V)} ay—E 
When expanded, this may be written in the form 
(p> E)(On EA oR ES) 0 ee (4) 


where g is the same function as in I(7) and is given by 


2 Bam alls é A >) 
[e(m, n)]?=3+2cos = +2 cos 7 +2 cos 2 Ta a. (5) 
g” ranges continuously from 0 up to 9. 

It is interesting that eqn. (4) is identical with the equations obtained by 
Roothaan and Mulliken (1948) in their equivalent calculations for the molecule 
borazole, provided that g? is given the values 1 (twice) and 4. Like them, we 
shall find it convenient to change the origin for energy to (xy+«p)/2, which we 
call Fy. Let us put 


E-—E,=Z, }(ap— ay) =8, B—foy+opS=y. 2S (6) 
‘These equations are the generalizations of I(5). ‘The secular eqn. (4) becomes 
Z? — 37 — g2(y— SZ PHO = ae ee (7) 

and the explicit solutions are 
Z =(1—g2S*)1{ — g2yS + yg? +(1—g2S2)84y2 P23}. (8) 


For any values of m/U and n/V (i.e. of g) the positive sign in (8) gives an energy 
value less than £y; the negative sign gives an energy greater than Ey. We may 
call these Z,{g) and Z_(g). As g* increases from 0 up to 9, Z, decreases from 
—6 to some other definite negative value, and Z_ increases from +8 to some other 
definite positive value. In saying this we have implicitly assumed that 8 is positive. 
This follows from the fact that nitrogen is more electronegative than boron, so 
that, as Roothaan and Mulliken (1948) have shown, ax <ag<ap <0. 

The argument above shows that the complete z-electron band is divided into 
two sub-bands, separated by an energy gap of magnitude 28. Now there are 
just sufficient z-electrons completely to fill the lower sub-band. Boron nitride 
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would therefore be an insulator. If we put =0 in (8), we recover the energy 
distribution found for graphite in I, where the two sub-bands just move to 
touch each other. The separation of the two sub-bands depends only on the 
existence of a non-vanishing 5. It would not be affected (except to a very minor 
extent) either by including second-nearest neighbour interactions, or more than 
one layer plane, or the distribution due to the o-bonds. 

We have tried to estimate suitable numerical values in order to draw the 
N(E) curve, as in I. The important parameters are the overlap integral S, the 
electronegativity difference term 5 and the ‘reduced’ resonance integral y. 
The value of S is easily calculated (cf. I) if we accept the usual Slater atomic 
wave functions. We find that S=0-22. This is very much the same as its value 
0-25 in graphite. For 5 and y the most sensible procedure is to adopt the values 
which Roothaan and Mulliken (1948) have found most suitable in their discussion 
of the energy and spectrum of borazole. This would make 5=2-3ev, and 
y=—2-3ev also. A possible alternative set of values suggested by Roothaan 
and Mulliken is 6=2-7ev and y=—1-8ev. The resulting difference in the 
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Fig. 3. Graph of N(E) for a single layer, overlap integrals included. 


N(E) curves is not large, and probably lies within the reliability of the method 
itself. For that reason we select the first set of values. Then the energy gap 
between the two sub-bands is 4-6ev. ‘This would be expected to represent the 
lowest amount of energy that could be absorbed electronically. A gap as large as 
this between the occupied band and the beginning of the conduction band would 
make the crystal a good insulator. 

For reasons given in I the top of the conduction band is most unlikely to be 
adequately predicted by this analysis. But the calculated width of the occupied 
sub-band would be expected to be more reliable. According to (8) this width 
is 2-4 ev, which we may compare with the value of about 5 ev predicted similarly 
for graphite. 

The detailed shape of the N(E) curve is easily found by combining (8) with 
the shape already calculated for graphite in I. ‘This is because the expression 
g(m, n) is the same both in (8) and in I, figs. 4(a) and 4(4), and in each case the 
variables are distributed uniformly in the range —}to +3. ‘Thus the distribution 
function for Z is simply 1/Z’(g) times the distribution function for g. Now 
this distribution function is, apart from a scale factor f, precisely the quantity 
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which is plotted in I, fig. 4(a). But dZ/dg is easily obtained from (8), or (7), and 
by this means the N(E) curve in fig. 3 is drawn. The two regions, which we have _ 
previously called Z, and Z_, give rise to the two sub-bands shown. Normally 
only the Z, band is occupied, and this band is full. The infinity present in the — 
graphite case still persists, and is removed in a perfectly analogous way by taking 
into account other neighbours in all directions. 

It is quite straightforward to include the parallel layer-planes in this analysis. 
The overlap integral between corresponding nitrogen orbitals on adjacent planes, 
calculated according to the prescription in I, is 0-011, and that between 
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Fig. 4. Graph of N(E) for three-dimensional crystal, overlap integrals included. 


corresponding boron orbitals is 0-056. "These may be compared with the value 

0-034 found in I for the carbon atoms in graphite. Following the argument 
that led from I, fig. 4 (8) to I, fig. 4(d) we may expect that the N(£) curve in fig. 3 

becomes the N(£) curve in fig. 4 which is parabolic at the tops and bottoms of 
both sub-bands. In view of our ignorance of the true values of the parameters 

this is probably about as far as it is worth going. It would be more satisfying if 

we could incorporate the o-bond electrons as well as the z-electrons. But the 

situation described at the end of I, where we attempted to interpret the x-ray 

emission spectrum of graphite, warns us that our knowledge of these o-bonds is 

still inadequate to justify any detailed calculations. 
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ABSTRACT. 'The7z-electron diamagnetic susceptibilities of polyacenes of up to ten rings 
are calculated by a method introduced in an earlier paper ; the method, which involves only 
a single numerical integration, is simplified. An alternative method of calculation, which has 
certain attractions and is widely applicable, is sketched in an Appendix. 

The relationship between the diamagnetic anisotropy and the number of rings in the 
molecule is found to be almost linear. 


§1. INTRODUCTION 

N this paper we shall deal with hydrocarbon chains of n rings (fig. 1), a typical 
] member of which is anthracene (n=3). Quite a number of molecules of 

this series have been prepared; consequently there is every prospect of 
assessing the accuracy and value of the theoretical treatment by direct comparison 
with experiment. ‘The polyacene chain is also of particular interest since the 
degree of condensation is intermediate between that of the singly linked polyphenyl 
chain (McWeeny 1951 a, containing Parts I and II, and to be referred to as I, IT) 
and of the highly condensed, hexagonal structures such as coronene (McWeeny 
1951 b, to be referred to as III) and we hope to find ultimately some relationship 
between degree of condensation and magnetic anisotropy. 


2 4 2n 
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Fig. 1. Polyacene of 1 rings. 


The method of calculation employed is essentially that of previous papers, 
involving expansion of the field-dependent secular determinant and a subsequent 
numerical integration (III), but a quite distinct approach (using second-order 
perturbation theory) exists and is briefly sketched in an Appendix. ‘The 
second approach is immediately applicable to any molecule whose zero-field energy 
levels and wave functions are already known and does not involve expansion of the 
secular determinant; but in the present paper, where a whole series of related 
molecules is considered, the expansion may easily be systematized, using recurrence 
relations, and the second method loses its particular attractions. 
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§2. EXPANSION OF THE SECULAR DETERMINANT 


Using the numbering shown in fig. 1 and taking the centre of symmetry of the | 


molecule as the origin of the magnetic vector potential it is possible to factorize 
the 2(2n + 1)th degree secular equation by choosing as basis functions the sets 


Pr = by rr be n+2—r)> Vr = br = Peont2—r) 
for these functions span subspaces providing different representations of the 


symmetry group of the system, C3. 
Thus, using the notation of II, with x=(E—«)/B and k=exp (27t«HS/12), 
S being the area of one benzene ring, we find the secular determinant of the 


n-ring molecule AWE DAD Go ae (1) 
where, for example, | 
SELEY +h" 
Re -— Re R* 
kx tk** 
Ro? =x k 
D+(4)= RY 4:1—8 ho SS SN Lee (2) 
Rs —x  k ; 
+ k4 RE x. R* 
Rito tsie SRP 
+8 R*a = % 


By suitable multiplications of rows and columns the k-factors may be removed 
from all elements except those lying on the “second diagonal’; in what follows we 
shall assume D is in this form. With this reduction it is possible to write 


=% 1 nee hehe 


D+*(n+2)= D*(n) 


+ gente 1 —Xx 
where z=k*6=e" say. There are two series of molecules to be considered 
according as is odd or even, but the same analysis applies in either case. It is 
also easily shown that D-(n; x)= —D*+(n; —x) so that D-(n) may be obtained 
from D*(n) simply by reversing the sign of every even power of x: we need now 
only consider D*t(n) and the superscript + or — may be omitted. 


Let D,(n) be the determinant remaining when 3” and z*” are struck out from 
D(n) and define also the following minors: 


D,(n)—by striking out from D(n) the first row and the last column, 
D,(n)—by striking out from D() the last row and the first column, 
3(”)—by striking out from D(n) the first row and the first column, 
4(n)—by striking out from D(n) the last row and the last column, 

;(7)—by striking out from D(n) the first and last rows and 


oO 


D 
D 
D 


the first and last columns. 


eee 
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From the symmetry of the array it follows that 
D,(n)* = D,(n) =D'(n), say, 
and D,(n) =D,(n)=‘D(n), say, (=‘D,(n)). 
With this notation and with D,(n)=D"(n), (=D "(n)) it is possible to expand 
D(n) in the form 
D(n) =Do(n) + 2"Do'(m) + 2*"Doy'(m)*—Dy"(n). va ee (4) 

In order to construct a series of D’s by recursion it is necessary to express each 

of the terms with suffix 0 as a combination of corresponding quantities associated 


with the preceding member of the series: the following relationships, derived by 
simple expansions, serve this purpose. 


Dy (n+ 2) =Dy'(n) —3*"Dy"(n), 

Do" (n+ 2) =x?D(n) + 2x‘ Dy(m) + Dy'(n), 

‘D(n +2) = —xDo"(n+ 2) +xD(n) +‘D,(n), 

Do(n + 2) = xD" (n + 2) — (2x? — 1) D(n) — 2x" D,(n). 

Thus, from Dy (7), Do’(n), ‘Do(n), Do(m), we may construct Dy'(n+2)....Do(n +2), 
using (4) to get D(m) and D(n+2). The ‘starting values’ are, for the even ring 
molecules, 

Di(0)=f, D,y"(0)=0, *D(0)=1, D,(0)=(t—x) ~ :..0.. (6) 
(being fixed by convention so that D(2) is given correctly) and, for the odd ring 
molecules, 


eee) (lyfe od). Dt) = ax 2) se (7) 


It is possible to make quite rapid progress by this method: thus five recursions 
would yield D(11), the expansion of a determinant of 23 rows and columns. 


§3, EXPRESSION FOR THE SUSCEPTIBILITY 

In (III, §3) it is shown that the z-electron susceptibility may be determined 
from A(x) by a numerical integration (III, (10)). Whilst this method is applicable 
generally, it would be extremely cumbersome in the present case owing to the form 
- of the integrand (the various polynomials involved reach the forty-second degree 
if completely expanded, and expansion would be desirable since various derivatives 
occur). The integration may, however, be thrown into more elegant form and 
the subsequent numerical treatment greatly simplified. 

Suppose the polynomial A(n; x) =/(x), say, is written 

Oras 0] oe teens (8) 

When A refers to the whole molecule (or represents a factor whose roots include 
all members of any initially degenerate set of levels) this form is quite general 
since no terms in @ can occur: this follows easily from I, Appendix If. 

Also ,for alternant hydrocarbons of the type we consider, A is an even function 
of x (Coulson and Rushbrooke 1940). Starting from the result (III, (7)) for 
the sum of the positive roots of f(x), 


ee ed ane d (9) 


+ In this connection an error in the published equations for coronene (III) should be noted: 
the equation labelled F, should differ from I’, only in the sign of the # term and similarly with I, 
and DF F 


842 R. McWeeny 


we put, f being an even function, f(iy) = F(y) (a real function of y) and have 
iyf'(ty)=yF'(y). In anticipation of subsequent manipulations of the integrand 
we then write (9), more cautiously, as 
Luli ee \ 
Sys eal {Pay —h ee (10 
i 7 fs o | Fy) e 
Integrating by parts 


R 
etd 2 yes {Iylog FO) { log F(y) dy —nR, ie (11) 
Ro 0 


and with F(y)=F,(y)+ 62G(y), @2G(y) being always very much less than F)(y), 
we find, on expanding the logarithms and carrying out the limiting process, 


Sea ln + (- dis ah ; dy) f) See (12) 


=x denoting the sum of the roots in the limit +0. 
Thus, in appropriate units, the z-electron susceptibility, and hence the 
diamagnetic anisotropy Ay, is given by 
1° Gy) 

Seat |G 5 dS Pa eee (13) 
Since no derivatives occur in the integrand, there is no need to expand completely 
the polynomials G(y) and F,(y): in the present application, for example, these 
quantities are conveniently expressed in terms of polynomials of lower degree. 


The Effect of Overlap 

The inclusion of non-orthogonality of amount 6 between adjacent atomic 
orbitals leads only to slight modifications of (13). Starting from (III, (7)) it is 

a simple matter to show that 

35 /* 14-45%" 96) 

Ax =(Ax)s—0- a J) C52" Fy) OA Gere 
Here (Ax),s_9 is a purely formal quantity since when overlap is included a change 
is implied in the empirical parameter 8 which determines the units of Ay: the 


16-0 


-G/F, 


0 0-5 1-0 
a 


Fig. 2. G/F, for n=6. 
terms in 6” are therefore not to be regarded as corrections to (Ay)s_9. Nevertheless, 
the form of G(y)/Fy(y) (fig. 2) allows us to conclude that the term involving the 
overlap directly makes only a small* contribution to Ax and will, therefore, have 
little effect on comparisons of Ay between two molecules: in fact, in all the cases 
considered the second term in (14) never exceeded 5°, of the first. 
* With the overlap usually assumed, 82=0-0625. 
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§4. RESULTS 

The numerical integrations involved in (13) and (14) are easily performed 
using simple integration rules, the form of the G(y)/Fo(y) curve being particularly 
satisfactory for this purpose (cf. III, where the integrand exhibits a sharp peak in 
the middle of the range). 

It is interesting to note that the method of calculation commonly employed in 
previous papers (London 1937), where modifications of the individual energy 
levels are separately considered, becomes impracticable when the degree of the 
polynomial A(x) becomes sufficiently high, as in the present case: the quantities 
determining the perturbations of a level (e.g. A’(x)) then occur as the differences 
of very large and nearly equal terms, with the result that seven- or eight-figure 
accuracy of the unperturbed roots may give only two-figure accuracy in Ay. 
This difficulty never arises in the contour integral method and, in addition to the 
fact that no zero-field energy levels need be calculated, the accuracy of the final 
result is limited only by that of the numerical integration; for practical purposes 
only three or four figures need be carried. 


1 2 3 4 5 6 8 10 
‘0 DAG 3-458 ee foe e600" 7:39 eel OO/m al 2 ot 
a0) 222 Ses 531) et OS 020 0-04 10-460 3-33 


n 
Ay (without overlap) 1 
Ay (with 8=0-25) 1 

Our calculations are summarized in the table where Ay is given for a number of 
molecules (with and without overlap). No discussion will be attempted at this 
point, for it is more illuminating to make a wider survey of the results obtained in 
these papers and elsewhere (e.g. Berthier et al. 1951, 1952); this will appear in a 
short concluding paper. 

It is worth noting here, however, that Ay is an almost linear function of , 
the linearity improving as m increases. Over a limited range the curve may be 
accurately fitted by Ay =a+6n+cn?: the coefficient c, which measures the depar- 
ture from linearity, then falls from 0-04 in the range 1 <n <4 to 0-007 in the range 
4<n<10. With overlap 5=0-25, the data in the table may be fitted to within 
less than 1°% by the formula 


Ay = —0-47+1:31n+0-007n2 — (n >4). 


APPENDIX 


A Perturbation Method 

The z-electron energy levels and wave functions in the absence of a magnetic 
field are known for a great many molecules: in some instances they may be very 
easily obtained using symmetries which are destroyed by the field. ‘The molecules 
considered in this paper come within this latter category and it seems worth while, 
in view of other possible applications, to outline a perturbation treatment based 
upon a knowledge of the zero-field solutions. It also seems worth while recording 
the complete solutions (eigenvalues and eigenvectors) for the general polyacene, 
on crude molecular orbital theory. In the absence of the field the secular deter- 
minant of an n-ring polyacene is conveniently factorized using a plane of symmetry 
and, choosing basis functions 


u,=¢,+ by, U,=h,— Pry 
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Coulson (1948) has shown how the factors may be evaluated: but in our case the | 
eigenvectors, U,=,U,d,,, Vi,=%X,0,b,, are required and it is more profitable — 
to examine the difference equations for the coefficients @,,, b,,. By suitable 


manipulations the equations may be thrown into a form satisfied by 
rko 
Gy, =sin=—— | x, r even 
rk 2(n fe 1 ) / k ( i 


Sees w —4)H2 
stingy |e 1) (r odd), 


where x, is a root of x(x, — 1) =4 cos? k7/2(n + 1), with rather similar solutions 6,;. 
These equations give 2x U’s and 2n V’s, two more solutions resulting from the 
special cases k=n +1 and completing the full set of 2(2”+ 1). 

Now the field-dependent determinant factorizes only when we choose the 
basis functions p,, g, (see §2). If H is the matrix of the array remaining when the 
x’s are struck out of D+ (and similarly with D~) then D+ =0 expresses the com- 
patibility of the equations, in matrix form, HC =xC and a perturbation treatment 
must be based on a knowledge of the solutions of Hj)C=x«C, H, being the form 
taken by H when 6-0. It is, therefore, necessary to express the zero-field — 
eigenvectors (15) as combinations of the new basis functions p,, q,: it is easily 
seen that this in fact involves only a re-classification of the eigenvectors U,, V;, 
(i.e. the number sets a,,, b,,) as P’s or Q’s. 

Arranging the zero-field eigenvectors into square matrices 


P = PLE Pp i. Q= |  OPGOE es 


and writing H=H,+70H’+6?H” (expanding the complex exponentials) the 
perturbation of the 7th energy level x, may be written, to the second order in 6, 
Ax, =O pp’ + O72X (Hr) : 
s Xp—Xs 
where %’ = P'H’P, %" = P'H’P, and P? is, in this case, simply the transpose of P. 
In the present case the matrix elements are easily written down: 


KE ys, oer 2[nP, Mies (n—Z)P5,.Psa0 8 ett: Pearl asl 


(the upper or lower sign being taken according as (r odd, s even) or (r even, s odd), 
HA ,., vanishing otherwise) : 


H yy" = +[n®P,,2+(n—2)?Pg2+ ... + Pea? 


(according as r is even or odd). It should be noted that the number of terms in 
these expressions is quite small (e.g. five terms for the 42-atom molecule) so that 
the computation is not unduly cumbersome; this result appears to be fairly 
generally true, i.e. by manipulation of the original determinant the number of 
complex elements may always be enormously reduced, only a few elements then 
remaining in H’ and H” and the calculation being kept within manageable 
dimensions. 

The perturbation of a root of D- =0 is calculated in exactly the same way, 
with a sign reversal in H’ and H” and the substitution of @ for P: actually no new 
matrix elements need be calculated, the resultant matrices being identical with 
those already computed (except for a reversal of order of the rows and columns). 
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Finally, in the same units as in §3, Ay =»,Ay, where the summation extends 
over all the (doubly) occupied energy levels and 
(Hee. eo 


Rien 


Ax, Toe eres 
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ABSTRACT. The various approximations which have been used to calculate the 
probability of excitation of an atom by impact with an electron whose energy is not greatly 
in excess of the threshold are known to fail under certain conditions. This is especially 
true in the important cases in which the excitation proceeds mainly through electron 
exchange. In these cases the exchange coupling is too strong for approximate methods, 
which assume weak coupling, to be valid. By using a schematic model exact solutions 
for strong exchange coupling are obtained and compared: with the results of different 
weak coupling approximations. It is found that the Born—Oppenheimer approximation 
is never reliable at low electron energies but if distortion of the incident and scattered 
waves by the atomic field is allowed for correctly much better results may be obtained. 


§1. INTRODUCTION 

T was first pointed out by Mott (1932) that the contribution to the cross section 

for an inelastic collision between two interacting systems which arises from 

impacts in which their relative angular momentum is {/(/+ 1)}/ can never 
exceed (2/+1)A?/47 where A is the wavelength of the initial relative motion. 
This follows from quite general arguments and is the quantum theory version of 
the obvious classical result that, if the relative angular momentum lies between 
J and J+dJ, the target area is 27-J dJ/ Mv" where v is the relative velocity and 
M the reduced mass. 

Many applications (see Mott and Massey 1949) of this limit theorem have 
been made in the interpretation of nuclear collision phenomena in which the 
interaction is strong and the usual perturbation theory (Born’s approximation) 
breaks down completely. For collisions of electrons with atoms the interaction 
is weaker and comes into play more gradually and more reliance has been placed 
on methods such as that of Born in these cases. Recently, however, it has been 
pointed out by Bates, Fundaminsky, Leech and Massey (1950) that there are many 
instances in which calculations of cross sections for excitation of atoms by electrons, 
based on current approximations, have given results in excess of the allowed 
maximum. In all these cases the large values obtained are associated with the 
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possibility of electron exchange occurring in the collision and arise at electron 
energies near the threshold for the process. It appears in fact that a calculation, 
in which the usual approximations are made, is likely to give more than the 
allowed maximum value for the cross section near the threshold of excitation 
of one level from another if both belong to the same electron configuration and 
may do so in other cases in which there is no change of azimuthal quantum number 
of the atomic electron concerned in the process. 

Unfortunately it is just under these conditions that information about 
excitation cross sections is required for the interpretation of astrophysical and 
geophysical phenomena. ‘Thus of the examples of failure noted by Bates e¢ al. 
(1950) the excitation of O?+, O+ and He (21S and 23S) are of major importance 
in the interpretation of nebular spectra while the deactivation of O('S) and 
O(!D) by electron impact would, if rapid enough, be important in the emission 
of the air glow and aurora. Another excitation process which is important in 
this connection is the dissociation of hydrogen by electron impact due to excitation 
of the *X, and 3X, states. This process provides the most effective means of 
obtaining atomic hydrogen. A knowledge of its probability under different 
conditions would be very valuable, but here again calculations carried out using 
the Born—Oppenheimer approximation give results exceeding the possible 
maximum. 

Before proceeding further it is necessary to distinguish between the various 
approximations which have been used in the calculation of inelastic cross sections 
for electron—atom collisions. In the Born approximation it is assumed that all 
matrix elements of the interaction between the electron and atom, whether 
diagonal or non-diagonal, are small and the possibility of electron exchange is 
ignored. ‘The cross section is then proportional to the square of the matrix 
element of the interaction between initial and final states in which the free electron 
wave functions are plane waves and the overall wave function is simply an 
unsymmetrical product. 

The Born—Oppenheimer approximation includes the possibility of electron 
exchange but still assumes that all matrix elements are small. The overall 
wave functions used to calculate the probability amplitude are properly sym- 
metrized combinations of plane waves and atomic wave functions. 

As a next step it may be assumed that only the non-diagonal matrix elements 
are small. ‘The plane waves which represent the initial and final free electron 
wave functions are then replaced by waves distorted by the mean interaction 
with the atom in the initial and final states respectively. If exchange effects 
are ignored in calculating the distortion we shall call the approximation the 
distorted wave Born (D.W.B.) approximation. This is to be distinguished 
from the distorted wave Born—Oppenheimer (D.W.B.O.) approximation in which 
the effect of exchange in producing distortion is allowed for. 

Of the cases cited by Bates et al. in which the calculations yielded cross sections in 
excess of the possible maxima, H,(4X,+—>*X,,+), He(1!S — 23S) and Hg(6!S + 68P) 
were calculated using the Born—-Oppenheimer and O(23P-+2!D and 21'S), 
O*(24S > 2?D) and O?+(2?P + 2!D and 21S) using the D.W.B. approximations. 
It is natural to enquire whether, in general, there is any advantage gained by 
using the D.W.B. or the D.W.B.O. method and whether, if the cross section as 
calculated by either method comes out to be in excess of the maximum, the true 


cross section is actually close to the maximum. ‘The same question can be put 
in a slightly different way. 


———— ae 
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- We may measure the strength of the interaction by a parameter A. Both the 
Born—Oppenheimer and D.W.B.O. approximations will give a cross section Q 
increasing as A? but the constant of proportionality will in general be different 
with the two methods. Which of these gives the best results for small A? 

In any case the true cross section cannot increase indefinitely as A?. We 
wish to know how Q varies with A when A is too large for the approximate formula 
to be valid. Ifit should turn out that, as A increases, either approximate formula 
‘continues to give a fair approximation until it gives O=Q,,, (the allowed 
maximum) and that fora considerable range of A beyond this the correct Q remains 
close to Q,,,x, then we have a result of immediate application. In most actual 
circumstances it is probable that the effective value of A is not an order of magnitude 
‘greater than that for which the approximate 0 =Q,,,,, i.e. the coupling though 
not weak is not very strong. 

A reasonable indication as to the way these questions may be answered may be 
obtained by studying a schematic ‘model’ of the process of excitation involving 
electron exchange. While the essential features are retained such a model may 
be set up which involves exactly soluble coupled integro-differential equations. 
In this paper an account of an investigation on these lines is given and very useful 
conclusions are arrived at which should make it possible to predict without serious 
error the cross sections for excitation processes in which electron exchange is of 
primary importance. It may be said in anticipation that it appears highly 
probable that the predictions of the D.W.B.O. method are not seriously incorrect 
if O<430Q,,,x while if it predicts a value of O considerably greater than Q,,,, the 
true value is close to Q,,,,. On the other hand, if distortion is not allowed for, 
little use may be made of the results for electron energies near the threshold. 


§2. THE INTEGRO-DIFFERENTIAL EQUATIONS FOR THE PROBLEM 


The problem of the excitation of the 2s level of atomic hydrogen by electrons 
has recently been considered by Erskine and Massey (1952). ‘They find that, 
for electron energies near the threshold, the cross section for the excitation is 


given by 
OzAorrR:| Bika ate SAb We) ala aly (hele os (1) 


k, and k, are the wave numbers of the initial and final motion of the electron 
relative to the atom. |f| is such that the asymptotic form of the solution /; 
of a certain pair of integro-differential equations 1s 


WS sO ieee a or (2) 


Since the maximum possible value for Q is 7/ky” we may call 47k,|f|?/ky the 
probability of the particular inelastic collision concerned. 
, The integro-differential equations expressed in atomic units are of the form 


d? 7 oO 4 ; ; ‘ 

| sa 2Volr) the? | fot [Kaa fale") 
=Welr flr) — | Koulrsr ile’) a 

oe Vale) tht | fit | Kuler Vile a 


Cie 
=r) flr)— | Krolrs fale”) 
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and the solutions must be proper functions satisfying the asymptotic conditions 
(2).and (ae Oe (4) 

In these equations Vy) and V,, are the direct interactions between the electron 
and the atom arranged over the initial and final states of the atom, Ko) and Ky, 
are interaction: kernels which represent the contribution of electron exchange 
to the mean interaction in each case, Vo, is the direct non-diagonal matrix element 
of the interaction and K,,, Ky, are the corresponding contributions from exchange 
effects. If V,,.Ko, and Ky are zero the probability of the transition vanishes. 

If the equations (3) are solved exactly the resulting cross section should be 
accurate, provided that the wave numbers ky and k, are sufficiently small for the 
contribution to the excitation from incident electrons with angular momenta 
greater than zero to be ignored and provided that the influence of other excited 
states is also negligible. As far as the first of these is concerned it is quite clear 
that almost all, if not all, of the violations of the conservation law occur for electron 
energies near the threshold in which case the main contribution comes from 
head-on collisions. The second approximation amounts to a neglect of dynamic 
polarization effects. At the electron energies concerned these are probably 
not very important though it is difficult to estimate their order of magnitude. 

In any case the exact solutions of eqns. (3) will certainly give cross sections 
which obey the conservation laws. It is only when the equations are solved by 
approximate methods that violation of these laws can occur. All approximate 
methods used hitherto proceed by solving eqns. (3) by successive approximation 
on the assumption that V),, Ko, and Ky, are small. The D.W.B.O. method 
makes no further assumption, but the Born—Oppenheimer approximation assumes 
that. Voose9 115 Non. ANG. Keane: aso negheie and the Born approximation 
neglects Kyo and Ko, as well. 

It is to be expected that the D.W.B.O. method will give accurate results if 
Vo,, Ko, and Ky, are indeed small but this is clearly not sufficient to justify the 
Born—Oppenheimer method. In fact, for excitation of atoms by electrons with 
energy near the threshold the distortion introduced by Vo) and V,, at least is 
very marked. ‘The detailed calculations carried out by Erskine and Massey 
(1952) illustrate this for the excitation of the 2s level of atomic hydrogen. 

To examine the usefulness of the various approximations it is very desirable 
to effect comparison with an exact solution of equations of the form (3). If the 
functions Vo9, Viz, Vox, Koo, Ky1, Ko, and Ky, have the actual values appropriate 
to an atom it is quite clear that it would be a very laborious task to obtain an 
exact solution. It is therefore necessary to simplify the problem by substituting 
simplified forms for the various interactions in (3) so that exact solution is possible 


in analytical terms. We now proceed to discuss the choice of such a schematic 
model for the equations. 


§3. SCHEMATIC MODEL ASSUMED 


‘The effect we are particularly concerned to investigate arises from the exchange 
terms Ko, Kj. It is this term rather than the direct non-diagonal matrix 
element Vy, which is responsible for the large cross sections given by the Born— 


Oppenheimer method for electron energies near the threshold. For the 1s—2s 
excitation of atomic hydrogen 


Kox(1, 7’) =4arbo(r) (Ro? — 2yo(7,1')r'h (1) =Kio(r', 7), see (5) 
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where yp and yw, are the respective wave functions for the 1s and 2s levels of 
hydrogen, ky,? =k? — + and y,(7,7’) is such that 


[PetCt oto rer Ir tCOSU). 9 Saas es (6) 


‘To represent the various features of these kernels we take 


ho(7) = (479?) 1*, r<79; Yo(7) =0, r>70, ] a 
v(t) =(4n79)*, r<r,; y,(r) =0, r>r. I Mernicts 


We further ignore ky,” and multiply the whole by a constant $4 which is equivalent 
to replacing the effective interaction e?/|r—r’| between the two electrons by 
A e?/|r—r'|. A can be taken then as the parameter which measures the strength 
ofthe coupling. Apart from A, our assumed kernel includes two other adjustable 
parameters 7) and7, which measure the spread of the initial and final wave functions. 

The direct interaction V5, is put equal to zero as the main effect with which 
we are concerned arises from exchange. 

The combined effect of Vo) and Koy in producing distortion of the incident 
electron wave is represented in the schematic model by a single function z(r) 
which is such that 


Bier Oy ry, OT) SO, Try. es, (8) 
A similar function v;(7) given by 
OA) ee iat 0, try See (9) 


replaces V,, and K,,. Over the comparatively narrow range of electron energies 
concerned the constants C and D may be chosen so that the distortion produced 
by the potentials vy and v, approximates quite closely to the actual distortion 
due to Vo, and Ky, and to V,, and K,, respectively. In any case we are concerned 
here with obtaining results of general interest and application rather than studying 
excitation of a particular atom. It is only important that the range covered in 
the investigation should be that in which particular cases are likely to fall. 
The equations which we use as a model representation of (3) are therefore 


lo yBfy = st 2 fir’)r' dr’ +r[ flr) ar eed Scere (10a) 


@ 
oP + byt fy=0, ety eee ere (10d) 
af, : ON. / = f f - > 
aa Yo f= | at ites Pah) Or a FSF, kt (10c) 
dr 0 
= Sll2 | : folr’)r' dr’, mor Sa Pgs ox a (10d) 
J0 
d* : 
oh + af, =0, ap ae te (10 e) 
As yp, refers to the excited state we have taken 7,>7). In these equations 
(pekptel, vee Re 2D, MY =3AKr or) 7, eens (11) 


Similar equations have been discussed previously by Massey and Mohr (1937) 
in connection with nuclear collisions. ‘They involve an extension to inelastic 
collisions of a model used by Morse and Allis (1933) to discuss the importance 
of exchange effects in elastic scattering of electrons by atoms. 
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§4, SOLUTION OF THE EQUATIONS 
For r9<r<r, we have 


fre=Bosinkr+Ceoskrrl,  ~— 9 9 “aqerns (12) 
1 =(0/v,2) i EG ew ana Oh ee (13) 
and, for r>7,, f, = exp (k,7). eeautnnity of f, and df,/dr at r=r, gives 
“B=if—Isink,r,  C=B—Icosky,  — .... (14) 
so for 7, =r 7, fi, =Bexp(ikir)+1{l—cosk(7,—7r)}. we ee (15) 
Bone hcrane fp a int de ley oe eee (16) 
s=1 
he ; CGC game hie ak AR as Pg (17) 
so that to satisfy (10a) Bi. 
Da,(vy2— x2) sin yr =A? { —r=b, = as +55, ae £ rs} Sort. (18) 


where 
j= [ [1{1 — cos k,(7, —7)} + B exp (2k,r)] dr 


=I fr, —r9—ky* sin ky(7, — 79)} — (¢8/Ry) {exp (thy) — exp (thyr)} (19) 
and, from (13) and (16), 


[=0"%a, (=e fers e820) Reb Eee Os ae (20) 
To satisfy (10c) ' 
2b, (v,2— x,2) sin y,7 = 1? { _rDa, en +a, a} egetint (21) 
s s 


Equating to zero the coefficients of ry on the right-hand sides of (18) and of (21), 


Xb, COs o = AUASy, = Yo ae at sin Ri(r, Py ro) }La, (Ae af ro mee xe) 
Ss Ss s 
= —iBk, {exp (tk,r,) — exp (tky7)}. - (22) 
a, oo ee (23) 
aXe 


Once these equations are satisfied the complete equations will also be satisfied if 
A,( V9" — Xo") = A?b,/x,°, 
bv? a Xe) = Sh a./yx.?, 

fors=1to4. This will be so if the y,? are the roots of the biquadratic in x?, 

Me — vy AOe— yy) =e. ee eee (25) 
Having determined the x, from (25) the ten unknowns a,, 6,, « and f are obtained 
from the four equations (24), the two equations (22) and (23) and the four contin- 
uity conditions at the boundary r=7y. Since for r>7o, fo must have the form 
(4) these four conditions are 


La, SIN X79 =SiN Roy + % Exp (éRoTo), 
La, X, COS X%9 = Ry CoS Rory + tkyx exp (tRoo); 

Xb, Sin x79 = 8 exp (tk) + 1{1 —cosk,(r, —79)}; 
Ub,V5 COS X19 =1k,B exp (tkyr9) —k, I sin k,(7, — 79). 
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Eliminating the b, by the relations (24) we may obtain in the form A,/Ay where 
A, and Ay are two six-row determinants. For given values of the constants 
A, ky, ky, vo, vy, 7 and r,, the procedure is to solve eqn. (25) numerically to obtain 
the x, and then proceed to the numerical evaluation of A, and Ap. 
When .9 is small we find for 8 the D.W.B.O. approximation 
fot: +: 
B= af 2p, - Ee {= (vo—%) _ sin eet x = {cos * (cos ne sin 4) 


He S “0 _ cos ts) (cos u, — cos “| st BFsfaat-s (27) 


0 


where 
K = {sin® uy +(v92/Ry2) cos? up}, L ={sin® 0, + (v,2/k,2) cos® “22. 
Ug = Volo» Uy == V119 and Uy = V41y. 


The Born—Oppenheimer approximation follows by putting vy =», v,; =, in (27). 


§5. NUMERICAL RESULTS AND DISCUSSION 

Detailed numerical calculations were carried out according to the following 
scheme. ‘The constants 7) and r, were first chosen to correspond roughly to 
atomic magnitudes. ‘I'wo separate cases were investigated. In the first of these, 
referred to as the ‘hydrogen’ model, the values chosen were 7) =1, 7, =2 and in 
the second, referred to as the ‘ helium’ model, 7) =0-6,7, =2. Detailed calculation 
of || as a function of ky was first carried out for the ‘distortionless’ cases in 
which C and D in (11) were put equal to zero. In this case the relevant small 
coupling formula is the Born—Oppenheimer approximation. The calculations 
were carried out for a number of different values of YW including particularly 
several values of such magnitude that the small coupling approximation gives 
results for the inelastic collision probability 4,|8|?/ky exceeding unity, the 
maximum possible value. For each value of .~ the maximum value of the 
inelastic collision probability considered as a function of ky was obtained and the 
curves Ia of figs. 1 and 2 illustrate the variation of this maximum value with 7% 
for the ‘hydrogen’ and ‘helium’ cases respectively. It will be seen that in each 
case the probability rises to a maximum of nearly unity at a particular value of .o/ 
and then falls gradually at greater values. Probabilities calculated by the Born— 
Oppenheimer approximation are illustrated by the curves Id of figs. 1 and 2. 
They agree quite closely with the exact calculations provided the probability is 
less than 0-6. 

In both figs. 1 and 2 the value of the coupling constant which corresponds 
roughly to atomic conditions is indicated. ‘This is given by putting A =2 in 
(11) so that, in absolute units, the effective interaction between the incident and 
atomic electrons has the true value e?/|r—r’|. 

The results obtained in the ‘ distortionless’ cases are encouraging in that they 
suggest that when the weak coupling approximation gives an excessively large result 
the correct probability is probably fairly close to unity. This satisfactory feature 
remains when distortion is taken into account, the relevant weak coupling formula 
being now that given by the D.W.B.O. approximation (27). Detailed calculations 
were carried out again for a ‘hydrogen’ and a ‘helium’ model differing from the 
‘distortionless’ cases in that C and D were no longer put equal to zero. For the 
‘hydrogen’ case C and D were so chosen that vg? =k? + 2:2, vy =k,? + 3-5, and for 
the ‘helium’ case so that vp2=hk,?+ 8-36, v,2=k,?+2-75. These choices were 
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based on the following considerations. ‘The values of vp are such that, for Re 
in the neighbourhood of the threshold, the phase shift produced in an incident 
wave by the potential well assumed agrees approximately with that already 
calculated for the actual atomic field (Morse and Allis 1933). v, was determined 
in the same way the phase shifts now referring to waves scattered by the respective 
atoms in the 2s state. For hydrogen the phase shifts are available from the 
calculations of Erskine and Massey (1952) but for helium they were estimated 
specially. . 
The results of the calculations, in which distortion occurs, are illustrated in 
figs. 1and2. The curves IIa represent the exact variation of collision probability 
with the coupling constant .°%, and the corresponding curves IIb the values given 
by the D.W.B.O. approximation (27). It is clear that this approximation is a 
good one provided the probability is less than 0-5. As the correct probability 
has a very flat maximum for a coupling constant which is much greater than that 
which would be expected to arise in actual atomic problems it also follows that 
when the D.W.B.O. approximation gives a result greater than unity the exact 
value is between 0-5 and 1:0. Knowledge that the probability in any actual 
case lies between these limits would be of considerable value in many applications. 
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Probability of an Inelastic Collision 


Provability of an Inelastic Collision 


A, 
ee 72 oe ee lin eee 
Coupling Constant .4 (atomic units) : Coupling Constant .4 (atomic units) 

Fig. 1. Probability of an inelastic collision as Fig. 2. Probability of an inelastic collision as 
a function of the coupling constant .o7 a function of the coupling constant .of 
(see (11))—‘ hydrogen’ case (estimated (see (11))—‘ helium’ case (estimated 
actual value of .of for hydrogen is indicated actual value of .o/ for helium is indicated 
at Go). at Go). 


Curve Ia. Exact calculation—no distortion present. 
Ib. Born—Oppenheimer approximation. 
Ila. Exact calculation—distortion present. 
IIb. D.W.B.O. approximation. 


Comparison of the curves Ila, IIb and Ib shows that the Born—Oppenheimer 
approximation (curves Ib) is completely unsatisfactory when distortion typical 
of actual atomic fields is present. ‘Thus the curves Id do not give a good approxi- 
mation to the actual curves IIa for any coupling constant. Entirely misleading 
conclusions would be drawn from them when distortion is present. We can 
therefore state quite definitely that the Born—-Oppenheimer approximation is 
wholly unreliable if appreciable distortion is present, no matter how small the 
coupling may be. In the cases considered it gives greater probabilities than the 
D.W.B.O. method but this is not a general result. Cases were easily found for 
which the D.W.B.O. method gives greater probabilities but remains more 
accurate. ‘I'he following procedure may be recommended for the estimation of 
the probability of an inelastic collision in which electron exchange predominates : 
Evaluate the probability using the D.W.B.O. method. If the value obtained is 
less than 0-5 it is probably fairly reliable. If it gives a probability much greater 
than 0:5 the true value is likely to be between 0-5 and 1:0. In carrying out this 
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prescription it is important that distortion should be fully allowed for. Thus, 
referring to the eqns. (3) for the excitation of atomic hydrogen, it is important 
to include the effect of the exchange interaction kernels Ky) and K,, as well as 
that of the direct interactions Vp) and V,,. 

The way the energy variation of the inelastic cross section varies with the 
coupling constant .¥ is also of interest. Figure 3 illustrates results obtained for 


° ) ° 
n w - 


Inelastic Cross Section in Units of tra, 
J 


1-5 2-0 25 3-0 35 
&? (atomic units) 


Fig. 3. Variation of the inelastic cross section with electron energy for different values of the 
coupling constant .ef—‘ helium’ case (exact calculation with distortion present). 


the ‘helium’ model. There is very little change of shape until .W exceeds 100 after 
which the maximum begins to move out to higher energies. ‘This effect occurs for 
values of .Y greater than that for which the maximum probability of an inelastic 
collision is attained (see curve Ila of fig. 2). For most atomic problems it is 
probable that the D.W.B.O. approximation gives a fairly good representation of 
the actual energy variation of the cross section. 
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Note on the Fluctuation Problem in Cascade Theory 
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ABSTRACT. A general relation between two fundamental functions appearing in 
fluctuation problems in cascade theory is derived. This leads to a direct method of solving 
for the 7th moments of the distribution function. It is also shown that the Janossy 
G-equations are superfluous for the solution of the fluctuation problem. 


Si INTRODUCTION 
W E have given the complete analytical solution of the fluctuation problem 


arising in nucleon cascade theory both for homogeneous nuclear matter 
and a finite absorber (Messel 1952, Messel and Potts 1952, to be 


——— 
o 


cee dle 


referred to as I, II). Since doing so we have established a number of general ~ 


relations which greatly simplify the previous work and which at the same time 
make the so-called Janossy G-equations superfluous for the solution of the 
fluctuation problem (Janossy 1950). It will be recalled that Janossy introduced 
his G-equations in an attempt to solve the fluctuation problem in nucleon and 
electron—photon cascade theory. ‘These equations, though appearing very 
simple in form, have proved to be of no value whatever in obtaining the actual 
distribution functions themselves (I, II). However, until the present, it appeared 
that differential equations for the nth moments of the distribution functions could 
only be obtained from either Janossy’s G-equations or Bhabha’s (1950) method 
of utilizing the energy as a continuously varying parameter. ‘The partial integro- 
differential equations for the moments so obtained have not been solved directly, 
though a recurrence solution may be given. On the other hand we have given 
in I and II the complete solution for the moments, by an indirect method with the 
Janossy G-equations as the starting point. 

It is the purpose of this paper to show that the solutions for the moments may 
be obtained by an alternative general method, which is straightforward and 
direct, using a simple relation between two fundamental functions appearing in 
cascade theory. No use is made of the Janossy G-equations. Though through- 
out this paper we shall for convenience consider the distribution functions which 
were given in II, the arguments are of general applicability and may be applied 
in particular to the fluctuation problem in electron—photon cascade theory. 


§2. FUNDAMENTAL RELATIONS 

As in II- we define H,(E,;£,,...,£,;0)dE,...dE,, as the differents 
probability that after depth @ a primary nucleon of energy E, has given rise to 
n nucleons with energies in the ranges E,, E; + dE,,j=1,...,m in any order, and 
d(£y; £,N;0) as the probability of finding N particles with energies greater 
than £ and an arbitrary number of particles with energies less than EF at a depth 0. 
Furthermore J,,(Eo; Ey,...,£,;0)dE,...dE,, is defined as the differential pro- 
bability of finding n particles with energies in the ranges E,, E;+dE,,j=1,..., 
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at a depth @ and any number of particles with arbitrary energies, due to a primary 
of energy Ey. The nth factorial moment of the distribution function ¢ is defined 
as T(E); E; 0). 

Then by definition 


A, oo 1 es) ce) E EH 
$(E,;E, N UD eeaaran aie aB...| dBy | dBy ss. | HIN 


x Hy +(Eo; Fy,... » Ey +539), ser g8\272 (1) 
T,(Ey; B36) = x a EOE Bde O seek pce vee k (2) 
eg hi. ees nis 9)= Ea aE, pis. sh FB 4 Wee At $d ree sled DR) 
c= 


We now establish the general relation we require. Substituting (1) and (2) 
we get 


oO © co E E 
T (Eo; E; @) = x ie dE, . ahs dE nia aE eee dE +a+b 
a=0 b= eS p 0 a) 
rl ea een ea 0). 2 |} + es SS Sok (4) 
Setting b=c—a (4) becomes 
co (oe) o fe c! 
T’,(Ey; E38) i By... | (Dee ica, 


* fe dE 41 os i. dEnta is dE n+a+1 $5 a iy dE cH n+( Eo ) E,, sions iDiea', ) 4) 


>) dE,...| a, 3 =| dP ce.) qlee, 
E E Pues /0 


MTR ER ea NO) e wee (5) 


Note that eqn. (5) gives a much simpler relation between H and T than (4). 
Now using definition (3) we find the rather surprising relation 


(Ey; B38) =|_ dB...) DRE Wem Feiler .%,-(6) 


-The same relation holds between the S*(e,x) and C,(e,...,€,;*) in I. It 
should be pointed out that the relation (6) may be readily eablehed from 
Ramakrishnan’s (1950) product density function. 

Hence equation (6) gives the nth moment expressed in terms of the distribution 
function J, and is obtained from J ,, by an elementary integration over the energy 
variables. Scott (1951) derived a similar relation for the special case of E=0 
corresponding to the physically unimportant quantity, the th moment of the 
distribution function for the total number of particles of all energies. 


§3. FIRST-COLLISION AND LAST-COLLISION DIFFUSION EQUATIONS 

It is well known that by considering just the first collision in a cascade a 

diffusion equation for the distribution function H may be obtained. This 

‘first-collision’ equation is a non-linear partial integro-differential equation and 

has not been solved directly. By considering all possible last collisions a further 

diffusion equation satisfied by H,, is obtained. The ‘last-collision’ equation is 
56-2 
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also a partial integro-differential equation but is linear, and we have solved 


it directly in II. : ; 
Just as for the case of the H,,, both a first-collision and a last-collision equation 


may be derived for the function J,,. The last-collision equation may be solved 
directly and hence using (6) the general solution for T,, is obtained immediately. 
In II the solution for J, was given as 


J (Ey; Exy.++; B38) 
—n Uy +70 +io rA+i0 
Ey | sees fe ds, | dp e?® 


~ (Qnijn 


Uy—40 Up—-to fico 
E ol E,.\snti 
‘ (z) ot (3) Y,(S1 se Sn :P), GA eat (7) 


where Y,, is a known function. 
From (6) therefore, it follows that (cf. II, eqn. (36)) 


Uy +i Upn+ico Atio 
T,,(Eo; E;) | see So dp e® 


(2772) " Sq, —600. 51 Un—t0 Sy J Aico 
i Sit... +8, 
x (3) DOP PE ey iy ue (8) 


It is now apparent that all that was required to solve the fluctuation problem 
was the solution of the last-collision diffusion equation for H, and J,. This 
method has also been used by us to solve the similar problem in electron—photon 
cascade theory. 

It should be noted that in the above work, no use is made of the Janossy 
G-equation. Since the only value of these equations originally was in obtaining 
the differential equation for the th moments, it now appears that there is little 
reason for their introduction into the theory. 
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LETTERS TO THE EDITOR 


Diffuse Reflection of Neutrons from a Single Crystal 


Diffuse streaks have been found in the Laue diffraction pattern of an iron single crystal 
placed in a beam of neutrons from the Harwell pile. 

The crystal, which had an effective irradiated volume of 0:025 cm*, was mounted with 
the [110] zone axis vertical, and surveys of the background scattering were made at various. 
fixed crystal settings. The intensity near the 110 reflection is illustrated in the figures. 
Figure 1 shows the counting rates along the equatorial layer line, on both sides of the 
Laue position, fig. 2 the counting rates on a line at right angles through A. The Laue 
reflection is seen to be accompanied by substantial wings, which cannot be attributed to: 
instrumental broadening. By subtracting the disordered scattering and the Laue peak, 
whose outline is known from the later studies, the excess counting rate resolves into the 
long hump of intensity shown. This is the (110) diffuse streak: it contributes up to 
140 counts per minute from vibrational and magnetic inelastic scattering processes. 
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Fig. 1. The equatorial distribution of background Fig. 2. A transverse survey through A of fig. 1. 
intensity around the Laue reflection 110 of Taken at 0p (110)=26° 50’. 


iron. Experimental points are indicated by 
circles. The residual intensity above the 
level of the disordered scattering is shown 
resolved into the Laue contour and the 
associated diffuse streak. The small peak 
at G is a Laue reflection from a disoriented 
inclusion in the iron crystal. @p(110)=19° 41’. 


As the Bragg angle of (110) is increased from zero, the diffuse streak first appears 
displaced towards the side of greater scattering angle and, in accordance with theory, 
gradually passes through the Laue reflection position. 

The intensity of the extra scattering is, as with x-rays, a few thousandths that of the 
elastic reflection. Between 292° k and 585° k it is observed to be closely proportional to 
temperature, as required by the equations of Weinstock (1944) and Moorhouse (1951). 

The wavelength changes which take place on scattering have been detected by 
absorption measurements at the settings A and B (with 0;(110)=32° 41’). It may be 
deduced from the geometry of scattering theory that neutrons undergoing first-order 
inelastic scattering associated with a given set of lattice planes gain energy if they are too 
slow for Bragg reflection, and lose energy if they are too fast. "This is strictly true for the 
vibrational scattering of slower-than-sound neutrons, and true for the predominant type 
of scattering in other processes. It follows that the equatorial survey of fig. 1 should 
almost completely separate out the two cases, accelerated neutrons going to the right of 
the Laue peak and decelerated neutrons to the left. The readings at A and B have 
verified this theorem. 

Weinstock’s theory of the vibrational scattering component predicts counting rates 
at C, D, E and F of about 40, 75, 105 and 55 per minute respectively. The observed 
streak is slightly more intense, and also of greater transyerse breadth, than thus estimated. 
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‘The difference corresponds to an inelastic scattering cross section of about ten millibarns. 
It is of the magnitude required by the theory of magnetic inelastic scattering, as discussed 
by Halpern and Johnson (1939) and by Moorhouse. 

Full details will be published in due course. 


Atomic Energy Research Establishment, R. D. Lowbe. 
Harwell. 
19th June 1952. 


Hacrern, O., and Jounson, M. H., 1939, Phys. Rev., 55, 898. 
Mooruowuss, R. G., 1951, Proc. Phys. Soc. A, 64, 1097. 
WEINSTOCK, R., 1944, Phys. Rev., 65, 1. 


Intensities in Band-Systems of O, and O,* 


In a note published recently (Pillow 1951) an improved method was suggested for 
distorting wave functions for a simple harmonic oscillator in order to give approximations 
to the vibrational wave functions for a diatomic molecule. ‘This method is here applied 
to the calculation of transition probabilities, and hence of relative intensities, for bands of 
two oxygen systems—the Schumann—Runge system of O,, and the second negative system 
of O.*. : 

In both systems Feast (1950 and private communication) has listed very fully the bands 
that have been observed, and which are presumably the strongest; and in the case of O,* he 
has (following Johnson) given visual estimates of comparative intensities, though without 
claiming high accuracy. In the tables below are set out the values for the intensities 
calculated on the assumption of uniform distribution among the excited levels—that is, of 
‘infinite vibrational temperature’, and the bands noted by Feast are indicated. ‘The 
intensity scale is arbitrary. 


Schumann—Runge Bands O, 


Excited State (constants calculated from data of Knauss and Ballard (1935)).* 
We’ =801°7, 0 a= 19°27, Te = 1-604 x 10— cm. 

Ground State (constants from Herzberg (1950)). 
®e = 1580, N12 ON, te =ai207 X10, com 


Calculated intensities at 7’ 00 


Die 5 6 7 8 9 10 11 12 13 14 15 16 


3 17 4S 67S eo O mn asORS, 5:6 dies 923 | 3:5 2:4 
F F F F F 


F=bands listed by Feast as having been observed. 


* ‘The method of distortion requires a potential function of the Morse form, and hence a 
quadratic expression for the energy levels. In most cases the measured values give such 
an expression to a very close approximation. For the excited state of the Schumann—Runge 
system the levels given by Knauss and Ballard need small terms in v+4 of higher order 
than the second to represent them exactly, but a quadratic expression has been fitted to these 
levels by the method of least squares and the errors are quite small. 
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Second Negative System O,+ 


Excited State (constants from Herzberg (1950)). 
@,’ =900, We == 118.4) Pe 1-409 x 10-* cm. 


Ground State (constants from Herzberg (1950)). 
We” =1876-4, to, == 16553; Fo I01227)< 10=2-cm: 


Calculated intensities at T’— 00 


0 1 2 4:5 6°5 8 a5 6 4 
(8) (7 (8), HIS) OO). .0).-3) 
1 0:5 2 5°5 9 10 9 5 1:5 1 
(3) (5 (6) (7) (8) (1) 
2 0:5 2) 6 11 13 7 4 it 3 3 
(1) (2) (7) (8) (8) (4) (2) 
3 2 6 11 15 2: + D 4 y) 
(1) (2) (8) (7) (6) (a) (a) (a) 
r 1 4 11 A7, 15 6 
(3) (6) (7) (7) (3) 
5 1°5 8 13 oN 12 2 ; Numbers in brackets represent 
(2) (8) (7) visual estimates made by Feast. 
(a)=band observed but not 
Po 3 13 23 19 6 5 estimated. 
(a) (2) (6) (4) (4) (5) 
7 6 20 28 15 3 3 7 
(a) (3) (6) (1) (3) (2) 
8 11 28 29 8 
(a) (1) (1) (1) 
Imperial College and M. E. PILLow. 


Northern Polytechnic, 
Holloway, London. 
- 18th July 1952. 


Herzperc, G., 1950, Molecular Spectra and Molecular Structure. I, Diatomic Molecules (New 
York: van Nostrand), 2nd edn. : 

Feast, M. W., 1950, Proc. Phys. Soc. A, 63, 557. 

Knauss, H. P., and BALLarpD, S. S., 1935, Phys. Rev., 48, 796. 

Pittow, M. E., 1951, Proc. Phys. Soc. A, 64, 772. 


Energy per Ion Pair for Argon with Small Admixture of other 
Gases 


Measurements of the value of the average energy w required to form one ion pair in a 
range of gases may be made by observation of the ionization current produced in a gas by 
the complete absorption of, say, alpha-particles of known energy emitted at a measured 
rate from a source. 

In the course of measurements of this kind on argon, air, Nz, Oz, CH, and CO, it has 
been observed that the ionization current in argon may be increased by about 1% (corre- 
sponding to a decrease in w) by the admixture of 2 to 5% of CO, or Og (other gases have 
yet to be tried). Since the value of w for O, and CO, is larger than that for argon (see table), 
and since the particle energy is completely absorbed in the gas, the effect must be due to a 
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more efficient conversion of the particle energy into ionization, or to a more efficient collection 
of the current, such as might be caused by an increase in electron mobility offsetting 
recombination. 

The latter hypothesis was rejected by observing that saturation was complete at the 
working value of the collecting field, and that the ratio w(A)/w(A+5%CO,) was constant 
at 1:01 over a range of collecting fields, and for alpha-particles from plutonium (thin source), 
uranium (thick source), and for beta-particles from *°S slowed down by foils to a maximum 
energy of approximately 40 kev. 

It would therefore appear that the contaminant increases the ionization current by the 
fruitful destruction of excited states, and that these account for about 1% of the energy 
abstracted from a charged particle in argon. 

In the measurements of w for A—O, mixtures by Haerberli, Huber and Baldinger (1950) 
there are no experimental points in the effective region, and the effect is not shown. ‘Their 
values were obtained from ion pulse measurements, and it would be of interest to determine 
whether the collection time in such apparatus is sufficiently long to allow for a high proba- 
bility of energy transfer. 

The absolute values of w obtained for Pu alpha-particles for the gases mentioned above 
are, to within 1%, 


Gas A Air O. N, CH, CO, A F5%CO; 
w (ev) 26-3 35-6 32-9 36-4 29-1 34-2 26-1 


This note is published by permission of the Director of this Establishment. 


Ministry of Supply, J. SHARPE. 
Atomic Energy Research Establishment, 
Harwell, Didcot, Berks. 


21st July 1952. 
HaerBerii, W., Huser, P., and BALDINGER, E., 1950, Helv. Phys. Acta, 23, 481. 


Paramagnetic Resonance in Potassium Chromicyanide 


Resolution in the paramagnetic resonance spectra of most compounds is limited by the 
line width arising from the magnetic interaction of the paramagnetic ion with the nuclear 
magnétic moments of the neighbouring diamagnetic atoms. In hydrated compounds the 
protons give rise to a half width at half intensity of about 6 gauss, and using heavy water 
this can be reduced to 2 gauss (Bleaney, Bowers and Ingram 1951). In a cyanide such as 
K3;Cr(CN)., where the neighbouring nuclear magnets are small and further from the 
paramagnetic ion, the limit should be reduced to less than 1 gauss. 

Single crystals were examined of potassium cobalticyanide, which is diamagnetic, 

containing between 1% and 0-1°% of chromicyanide. Since no paramagnetic resonance 
spectrum could be observed of the undiluted manganicyanide in the temperature region 
12°K to 290°K, manganicyanide was later also used as dilutant. In each case the spectrum 
oe to two ions in the unit cell, and can be fitted to the following spin Hamiltonian, 
with S=3/2: yo ghHS + D{S,2—4S(S+1)}+-E(S,2—S,) + ASI. 
Using cobalticyanide as dilutant, the experimentally determined parameters are g=1-:992, 
D=-+0-083 em~!, E=0-011 cm~, and with manganicyanide as dilutant the parameters are 
g=1:993, D=+0:054 cm, E=0:012 cem-!. The experimental errors are + 0-002 in g 
and + 0-001 em~'in Dand FE. The sign of D was found from the unequal intensities of the 
(3 <> 4) and (—} <> —3/2) transitions at low temperatures. The minimum line width, 
observed for the (<> —}) transition in a crystal of cobalticyanide containing probably 
rather less than 0-1% of chromicyanide, was 1:0+ 0-3 gauss. 

In the determination of the nuclear spin of *8Cr as 3/2 (Bleaney and Bowers 1951) an 
alum was used in which the spectrum was not fully resolved. Ina cyanide crystal containing 
chromium enriched in ®*Cr a hyperfine structure of four lines of equal intensity (see figure) 
was clearly resolved at 90°k and 20°x, confirming the value of 3/2 for the nuclear spin. 
From the separation of these lines |A|=0-00147+ 0-00005 cm-'; A is isotropic. 
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To obtain an estimate of the nuclear magnetic moment of **Cr it was assumed previously, 
on empirical evidence (see Bleaney and Bowers 1951), that the magnetic field acting on the 
nucleus is the same as for the isoelectronic ion V2+. Comparison of hydrated double 
sulphates of Cr°* and V** then gave a value of 0-46 n.m. A second value can be obtained 
by comparison of the above value of A with that of diluted K,V(CN),.3H,O, for which 
J. M. Baker (unpublished) finds |4|=0-0056 cm-. This gives 1 =0-58 n.m., and suggests 
that 0:5+0-1 n.m. is the probable value. 

Single crystals of undiluted potassium chromicyanide gave a single broad isotropic line; 
the value of g is 1:99+0-01. The line width was anisotropic, and its variation had roughly 


H 
oe 


Hyperfine structure of °*Cr in the electronic transition ($< —4) of diluted potassium chromi- 
cyanide. ‘The four lines are clearly resolved; the central line is due to even isotopes with 
zero nuclear spin. ‘The small line on the right (low field) is due to an organic free radical, 
used as a field marker. 

a 

the same angular dependence as the separation of the outer two lines found in the diluted 

chromicyanide. The absence of resolved fine structure, the shape of the line, and the angular 

variation of the line width suggest that there are considerable exchange interactions between 
neighbouring ions in the concentrated salt. 

The author wishes to record his gratitude to Dr. B. Bleaney for generous help and 
discussion. His thanks are also due to the Director and the Electromagnetic Separator 
Group of the Atomic Energy Research Establishment, Harwell, who supplied the enhanced 
*8Cr, and to Dr. P. F. D. Shaw for the preparation of various chemical compounds. 
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A Search for the Nuclear Reaction *He(a,n) 7Be 


According to the atomic masses listed by Li et al. (1951), the nuclear reaction 
4He(a,n)*Be should be endothermic to the extent of 18:98+0-04 Mev. ‘This means 
that, if an w-particle beam is used to bombard a helium target, the threshold energy for the 
reaction in the laboratory system of coordinates is 38:0+0:1 Mev. At the end of 1951, 
using the Birmingham 60-inch cyclotron, we made an effort to detect this reaction at an 
energy a little above the calculated threshold. This attempt was not successful, and we hope 
to repeat the experiment with some refinements at a slightly higher cyclotron energy when 
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an opportunity occurs. However, as this may not be for some time, it seems worth while 
to report now the upper limit which we established for the cross section of the proposed 
reaction. 

The layer of helium comprising the target was compressed to 3 atmospheres and was 
9-2 cm thick. The external w-particle beam from the cyclotron was admitted to the target 
through a half-inch diameter hole covered by a mica window 3 mg/cm? thick. A measure- 
ment was made of the energy of the a-particle beam, before it entered the target chamber, 
by absorption in aluminium. This gave 39-9 Mev, with an estimated r.m.s. spread of beam 
energy (after allowing for straggling in the aluminium) of 0:2 Mev. This estimate of energy 
is based on the theoretical range—energy relation of Aron et al. (1949), which is believed to be 
accurate to within }°%, at this energy (cf. Aron 1951). The energy agreed with an indepen- 
dent measurement of the range of the a-particles in a photographie plate. Owing to the 
loss of energy in traversing the window and the helium, the mean energy of the «-particle 
beam in the helium target was actually 39-0 Mev. This corresponds to an energy, in the 
centre-of-mass system, of 0:5 Mev above the threshold of the reaction *He(a,n) ‘Be. 
Beryllium nuclei produced by this reaction will therefore be strongly collimated in the 
forward direction in the laboratory system. Aluminium foils were placed behind the helium 
target and a search was made in these foils for the characteristic radioactivity of “Be. No 
such activity was detected, and from the conditions of the experiment an upper limit for the 
cross section of the proposed reaction can be calculated. This upper limit is 7 x 10~*5 cm?, 
which is about 10~? of the nuclear geometric cross section for the collision of two «-particles. 
(The upper limit was in part set by long-lived radioactivity produced directly in the foils by 
a-particle bombardment, and could probably be reduced appreciably by the use of foils of 
very pure aluminium.) 

A system of two a-particles can exist only in states of even angular momentum. If the 
reaction *He(a,n)*Be were to compete favourably (at energies near the threshold) with 
the elastic scattering of the two a-particles, the neutron would have to be emitted with 
zero angular momentum, and this would require the ground state of ‘Be to have even parity. 
The present result indicates otherwise, and agrees with what one would expect from the 
nuclear shell model for the ground state-of a nucleus of 4 protons and 3 neutrons, namely, odd 
parity (see for example Klinkenberg 1952). 

The reaction *He (a, n)*Be is the mirror of *He(a,p)’Li. The latter reaction is the 
inverse of the well-known reaction “Li(p,a«)*He, of which the excitation function 
(Rumbaugh et al. 1938*, Heydenberg et al. 1948) exhibits a maximum at a proton energy of 
3 Mev. From general considerations based on this excitation function it seems possible 
that the reaction “He (a, n) 7Be should be observable, in a more refined experiment with the 
present method, at an a-particle bombarding energy of about 41 Mev. 

We should like to thank Professor W. E. Burcham and Professor R. E. Peierls for their 
interest in this work. 


Physics Department, D. WALKER. 
University of Birmingham. W. T. LINK. 
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*'The uncertainty in the absolute cross section for the reaction *Li(p, «) ‘He may amount to 
as much as a factor of 3. 
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REVIEWS OF BOOKS 


Low Temperature Physics: Four Lectures, by F. E. Stwon, N. Kurtt, J. F. ALLEN, 
and K. MENDELSSOHN. Pp. vi+132. (London: Pergamon Press, 1952.). 
21s. 


In recent years, and particularly since the war, there has been a considerable growth of 
interest in research at very low temperatures. On the one hand it is a relatively unexplored 
field where there are still many phenomena which are only imperfectly understood, so that 
comparatively simple experiments can yield results of fundamental interest, and on the other 
hand the development of the means of providing the very low temperature—particularly 
the development of the Collins helium liquefier as a commercial product—has made it 
possible to carry out such experiments at an ever increasing number of research centres. 
This little book collects together the texts of four lectures on low temperature physics 
delivered at the Royal Institution in 1950 by leading authorities in this country; the lectures 
have been revised to bring them up to date to about the end of 1951, and provide an 
interesting popular account of the present state of the subject. 

In the first lecture Professor Simon surveys the whole field and emphasizes the special 
interest of low temperatures in making apparent quantum phenomena which at higher 
temperatures are masked by thermal agitation. His elegant account of the basic principles 
involved in the lowering of temperatures and the limitations imposed by Nernst’s theorem 
serves also as an introduction to the second lecture by Dr. Kurti who deals with the method 
of adiabatic demagnetization and the new knowledge of the magnetic behaviour of para- 
magnetic salts gained with its aid. The third and fourth lectures, by Professor Allen and 
Dr. Mendelssohn respectively, deal with the two topics which are perhaps the most fascin- 
ating—the properties of liquid helium and superconductivity. It is probable that both 
involve manifestations of quantum phenomena on a macroscopic scale, but the basic theory 
is still very incomplete. It is this and the intrinsic beauty of many of the experimental 
phenomena, well brought out in both lectures, which make these topics so fascinating. 

It is perhaps too much to ask for completeness in such a short account of a large and 
rapidly growing subject, but it should be noted that there is practically no mention of recent 
studies of the electrical and magnetic behaviour of non-superconducting metals at low 
temperatures—for instance of the high frequency behaviour of metals when the mean free 
path is very long, and of the peculiar field dependence of the magnetic susceptibility 
discovered by de Haas and van Alphen—studies which promise to throw new light on the 
fundamentals of the electron theory of metals. Otherwise the book may be confidently 
recommended as an attractive and lucid introduction which should be valuable not only to 
-anyone embarking on research in the subject, but also to the general student of physics who 
is looking for a broad picture unburdened by too many technical details. D. SHOENBERG. 


Verstandliche Elemente der Wellenmechanik, Vol. 1, by KARL JELLINEK. Pp. xii+ 
304. (Basel: Wepf, 1950.) 34 S.fr. 


This book, the first of two volumes, addresses itself primarily to experimental physicists 
-and sets out to provide yet another introduction to wave mechanics. ‘The usual topics are 
treated : a lengthy discussion of the wave and corpuscular aspects of light and of material 
‘particles leads to the introduction of Schrédinger’s equation, which is then used to discuss 
the simpler one-electron problems, including the harmonic oscillator, the rigid rotator, the 
hydrogen atom and various types of potential barrier. 

The author presupposes a minimum of mathematical knowledge on the part of his 
readers ; his preoccupation with elementary mathematical details is a little wearisome 
_at times, but no doubt many readers will regard this as a fault in the right direction. And it 
is good to find that the mathematical properties of the various types of wave function are 
fully and clearly discussed and illustrated by copious diagrams. Other welcome features of 
the book include the author’s habit of illustrating his arguments by means of simple order- 
of-magnitude calculations, and his comparison on several occasions of the predictions of © 
wave mechanics with those of classical theory and of the old quantum theory. 

The value of the book is unfortunately greatly diminished by the author’s eccentric 
ideas about the physical interpretation of wave mechanics. Not content with the accepted 
-probability interpretation of Schrédinger’s wave function, he attempts to picture the motion 
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of a particle in a pseudo-classical way which is quite out of keeping with the philosophy of 
the subject. To this end he is at great pains to postulate and justify the existence of a 
‘world ether’, basing himself on quotations from writers on the theory of relativity ; 
wave functions are apparently to be regarded as describing the ‘ ability ’ of the “ ether’ to 
transport particles or photons. One can imagine the bewilderment of the innocent reader, 
who on one page has learned for the first time that e#—=cos x-+i sin x, and who is asked a 
few pages later to consider (out of their context) some of Einstein’s more abstruse remarks 
on the philosophical basis of general relativity. The author’s attitude leads him to discuss 
whether or not particles ‘ cause’ wave functions and to ask ‘ how’ a particle can possibly 
move through points where its wave function is zero, or ‘ how’ a particle can ever be in a 
region where its kinetic energy is negative; to answer such questions he feels compelled to 
assume that the energy of a particle must be fluctuating rapidly owing to an exchange of 
energy between the particle and the ‘ ether’. 

The reader must be warned to beware of the author’s hobby-horse; if he can successfully 
avoid the evil eye of this strange animal, he will find much valuable information in the rest of 
the book. E. H. SONDHEIMER. 


Philosophic Problems of Nuclear Science, by W. HEISENBERG (translated by 
F.C. Hayes). Pp..126. (London:, Faber, 1952.) 16s, 


Under an inept title, this book is a somewhat halting English translation of a collection 
of Heisenberg’s essays on general aspects of physics. These essays cover a wide range of. 
subjects and stretch over a long period of time, the first dating from 1932 and the last 
from 1948, but they show a remarkable unity and continuity of thought. They deal mainly 
with historical and epistemological aspects of atomic physics, but they also include some 
consideration of the more directly human and social relations of science. One may take 
exception here and there to some interpretation of historical situations, as for instance to: 
the overestimation, fashionable in German circles, of Goetheé’s scientific dilettantism. 
One may doubt the author’s political wisdom when (p. 77) he describes Hitler’s savage 
onslaught as a “‘ reshaping’ of the world, which (this was written in 1941) “‘ absorbs all 


our powers’. But whether or not he agrees with the author’s views, the reader is sure to 
derive inspiration from an intercourse with one of the most powerful minds in the realm 
of physics. L. ROSENFELD. 


Introduction to the Transfer of Heat and Mass, by E.R. G. ECKERT. Pp. xiii+284. 
(New York, London: McGraw-Hill, 1950.) $4.00, 34s. 


There are plenty of books on heat transfer, either theoretical or practical in outlook, but 
I know of none which presents the best of both in the way this does. ‘The theory is quite 
straightforward in pure conduction and in pure radiation problems, and the numerical 
data are adequately known, but when convection plays a part the whole story is different. 
A few years ago the best one could do was to go to either of the two good books, one British 
and one American, in which experimental data were critically examined and collated in the 
light of similitude, and to take from them the nearest analogous case to the one in hand. 
One could also look at the mathematical attempts being made by men like Karman and the 
theoretical workers in aerodynamics, but they seldom contributed to a complex practical 
problem. ‘The third possibility, favoured by some chemical engineering teachers, was to- 
pretend that everything could be reduced to an equivalent problem in pure conduction, 
taking the numerical constants from experiment. 

‘The present book, designed primarily for chemical engineering students, treats the whole 
matter rationally on the basis of boundary-layer theory, thus providing logical reasons for 
the form of the basic equations. It has the additional merit that, when the form of the 
equation has been established, good values of the constants are given based on experiment. 
The student trained in this method will certainly see heat transfer as a field where a rational. 
and not a rule-of-thumb approach is possible and desirable, and he will almost certainly 
design better apparatus for having studied in this way. 

Naturally there are plenty of problems left where further study is needed, notably the 
problems of heat transfer between solid surfaces and fluids which are evaporating from or 
condensing on them. Here, the author gives a very good summary of the present position. 

taste. 
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ABSTRACTS FOR SECTION B 


A New Photographic Method of Measuring the Dispersion of the Optical 
Constants of Metals, by J. Bor. 


ABSTRACT. The information from which the optical constants are determined is 
produced on a photographic plate as a vertical band of variable intensity for each wavelength. 
This is compared microphotometrically with an adjacent band of graded intensity produced 
by the same wavelength. 


Photoelectric Analysis of Elliptically Polarized Light, by J. F. Arcuarp, P. L. 
Ciecc and A. M. ‘TAYLor. 


ABSTRACT. ‘Two photoelectric methods of analysis of elliptically polarized light have 
been devised. -One is a development of that of Kent and Lawson in which a rotating 
analyser yields a zero signal only when the light is circularly polarized. We use a double ~ 
image prism after a rotating analyser and the two resulting modulated beams are in anti- 
phase only when the axes of the double image prism coincide with the axes of the vibration 
ellipse. In our second method, the light is modulated by chopping, and a double image 
prism selects two components in phase. In principle these are equal only when the axes of 
the double image prism are at 45° to the axes of the vibration ellipse. ‘Thus the directions 
of the major and minor axes are found, and if the outputs from the two photocells which 
receive the two components are balanced on an electrical bridge the axial ratio may be 
obtained. ‘The double image prism plays here the part of the half shadow device in visual 
observations, but as its half shadow angle is 90° instead of only a few degrees, the photocells 
are used under optimum conditions, i.e. the signal-to-noise ratio is maintained as large as 
possible and the potential accuracy of the method greatly exceeds the limitations imposed by 
the rest of the apparatus. By slight modifications the method may be made partially 
self-recording. 

Some examples are given of measurements made on metallic films undergoing surface 
changes. 


The Absorption of Light by Noble Metals and its Relation to the van der Waals 
Contribution to the Cohesive Energy, by J. FRIEDEL. 


ABSTRACT. Noble metals show absorption in the blue or ultra-violet; to this is due the 
colour of copper and gold. As was pointed out by Mott, the main absorption seems due 
to internal photoelectric effect, with excitation of a d electron. ‘This absorption could be 
analysed—as for the alkalis—as a transition between two bands, described by Bloch orbitals, 
here d ands. But we must take into account the correlation between the hole created in the 
nd shell and the Fermi electrons. We can assume that the hole is screened by a (x+1)s 
electron bound to it, so that absorption should begin for the frequency of formation of the 
excited nd°(n+-1)s configuration of the monovalent ion. Agreement with experiment 
obtained using this model is satisfactory both in the pure metals and their alloys. The same 
model may be used for x-rays, and optical absorption energy is related to the energy separ- 
ation between the s and d parts of K and Ly, emission bands. An order of magnitude 
for the van der Waals interaction energy between d shells could perhaps be deduced from 
the structure on the short wavelength side of the absorption peak. . 


The Optical Constants of Thin Metallic Films Deposited by Evaporation, by 
P. L. CLEGG. 


ABSTRACT. The optical properties of evaporated films of silver, gold, tin and indium 
have been determined by a photoelectric method. 'The influence of factors concerned with 
the evaporation procedure has been investigated in the case of the silver films. 'The 
experimental results are accounted for by the agglomerated state of the films, and the 7 
application of Maxwell Garnett’s theory is shown to lead to qualitative IATA with the 
experimental values. 
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Intensity Limitations in a Point Source of X-Rays, by V. E. CossLetr. 


ABSTRACT. When the diameter of the focal spot of electrons on an x-ray target is reduced 
to the order of a micron, the conditions of heat dissipation become so favourable that a greatly 
increased specific loading is possible. 'The use of an electron lens to obtain such a small 
spot introduces a limitation on the rate of supply of energy, as its high spherical aberration 
requires insertion of a very small aperture, the size of which must be reduced pari passu with 
the size of spot desired. For magnetic lenses of usual design the rate of energy supply is 
then proportional to the radius of the spot. Comparison with the rate of heat dissipation 
shows that, for spots of the order of a micron, the limitation on x-ray intensity is not the 
thermal properties of the target, as in tubes with large focus, but the emission obtainable 
from the cathode. Specific loadings of the order of 10’ wattscm~? can be tolerated in a 
spot of 1 micron; the corresponding emission current density would be of the order of 
100 amp cm~-?. 


The Optical Properties of Liquid Selenium, by E. W. Saker. 


ABSTRACT. ‘The optical properties of liquid selenium in the near infra-red have been 
measured for a temperature range up to about 400° c. The absorption edge moves out to 
longer wavelengths as the temperature increases by an amount corresponding to a change 
of energy gap of 1-410 ev/c°. Absorption is very small at low frequencies. For 
long wavelengths the refractive index decreases as the temperature increases by an amount 
corresponding to the thermal expansion. 


The Contamination in fy ees Films by the Material of the Source, by O. S. 
HEAVENS. 


ABSTRACT. ‘Two methods have been used to detect contamination in evaporated films 
of silver and germanium by the boats or filaments used as sources. It is shown that a 
limit of a few parts in ten thousand of tungsten or molybdenum is attainable by microchemical 
methods. ‘The much lower limit of one or two parts in 10’ is attained by the radioactive 
tracer method. ‘The extent of contamination of films of silver and germanium by tungsten 
and tantalum boats has been studied using the tracer technique and the minimum contam- 
ination attainable found to be of the order of a few parts per million. The possibilities of 
transfer of the boat material by direct evaporation or by reaction with the residual gas in the 
system have been minimized; the results suggest that contamination arises by the solution 
of the boat materials by the molten silver or germanium. 


Self- Absorption in Arc Sources in Thermal Equilibrium, by H. EpELS. 


ABSTRACT. Ananalytical determination of the self-absorption of a spectral line emitted 
by a source is not normally possible if account is taken of the spatial variations in the source 
of the excited atom concentrations and spectral line profiles. "The problem is simpler for 
uniformly excited sources, for example sources in thermal equilibrium, but not for high- 
pressure arc sources normally assumed to be in thermal equilibrium, because of the existence 
of temperature gradients. In the paper an analysis is developed which shows that if such 
arcs are treated as uniformly excited, an error of less than 5°% is introduced into the deter- 
mination of the self-absorption of spectral lines, emitted by transitions between an upper 
state s and a lower absorbing state t, which satisfy the limits hvgq/RT, > 10 and vgt/rpq<0°5, 
where q is the ground state and J, the axial source temperature. ‘The theory is then applied 
to the experimental data for a high-pressure mercury arc discharge and to new measurements 
of the line widths at 5461 and 4358 A. It is shown that the calculated value of the 
self-absorption at 5461 A agrees with the experimental value only when the line has a 
resonance profile. The existence of this form of profile is also indicated by considerations 
of the broadening processes in the discharge. 


Ring Shims for Coned Magnet Polecaps, by E. R. ANDREW and F. A. RusHworTu. 


ABSTRACT. The calculation of ring shims of the kind suggested by Rose has been 
extended to include polecaps in the form of truncated cones. Design data are presented in 
graphical form. Actual and predicted performances are compared for one magnet. 
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Observations on Growth and Etch Phenomena on Haematite (Fe,O3) Crystals, — 
by A. R. VERMA. +} 


ABSTRACT. New observations of growth and etch features on haematite (Fe,O3) — 
crystals are reported. Of the many growth patterns observed, the simplest are the “i 
elementary ‘ molecular growth spirals’ in accordance with the theory of Burton, Cabrera 
and Frank. The step height between the successive arms of an elementary spiral has been — 
measured by the application of multiple beam interference fringes and is shown to be equal 
to the x-ray value of the repeat distance. More complex growth patterns are illustrated and 
discussed. The etch figures reported are also of molecular thicknesses and consist of a very 
large number (~107/cm? of the crystal surface) of similarly oriented small triangular figures, 
and some more complex shapes. 


On the Theory of Kinking, by F. C. FRANK and A. N. STROH. 


ABSTRACT. The concept of a kink-band is introduced. It is a thin plate of sheared 
material in a crystal, transverse to a slip direction, bounded by opposite ‘ tilt walls’ of 
dislocations. It is shown that if the angle of kinking shear y exceeds a critical value 
(estimated to be 34°) the stress concentration at the edges of the band is sufficient for the 
creation of new dislocation pairs or loops, and that if the applied stress exceeds a certain 
value, inversely proportional to the square root of the diameter of the kink-band, the kink- 
band will grow by this process. The stress field of a kink-band is similar to that of a slip- 
band, a plate of deformation twin, or a crack under tangential stress. “In suitable ~ 
circumstances these may act as sources for each other. ‘The successive generation of 
kink-bands, inside each other, from the same source, will produce a kink of the kind 
experimentally observed. 


